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Abstract 


The  use  of  wind  power  is  in  a  period  of  rapid  growth  worldwide  and  wind  energy  systems  have 
emerged  as  a  promising  technology  for  utilizing  offshore  wind  resources  for  the  large  scale 
generation  of  electricity.  Drawing  upon  the  maturity  of  wind  turbine  and  floater  technologies 
developed  by  the  wind  energy  and  oil  and  gas  industries,  respectively,  large  offshore  wind  energy 
systems  have  been  developed  and  are  being  proposed  for  operation  in  offshore  areas  where 
environmental  restrictions  are  ^less  restrictive,  large  wind  resources  exist,  and  open  sea  areas  are 
available  for  wind  farm  development. 

I  ' 

A  fully  coupled  dynamic  analysis  technique  was  developed  to  predict  the  response  of  a  floating  wind 
turbine  system  in  a  stochastic  wind  and  wave  environment.  This  technique  incorporated  both  non¬ 
linear  wave  loading  on  the  submerged  floater  and  the  aerodynamic  loading  on  the  wind  turbine.  A 
tension  leg  spar  buoy  was  designed  to  support  the  wind  turbine.  This  design  was  chosen  due  to  its 
relatively  small  size  and  hence  lower  potential  cost  per  wind  turbine.  The  system’s  tethers  were 
attached  to  the  ends  of  spokes  which  radiated  out  from  the  spar  cylinder.  This  arrangement  of  lines 
and  spokes  promised  to  be  very  stiff  in  the  roll  and  pitch  modes  of  motion. 

The  fully  coupled  analysis  technique  was  used  to  evaluate  the  feasibility  of  the  chosen  floater  design; 
Damping  properties  of  the  combined  floater  /  wind  turbine  system  were  determined  by  conducting 
simulated  free  decay  tests  for  the  different  modes  of  motion  and  wind  turbine  operating  conditions. 
Numerical  simulations  for  operational  conditions  were  also  carried  out.  The  response  of  the  floating 
wind  turbine  to  three  different  sets  of  environmental  conditions  was  determined  and  compared  to  a 
fixed  base  system.  Additional  simulations  were  conducted  to  determine  extreme  wind  and  wave 
event  response.  Stiffness  of  the  floating  system  in  roll  and  pitch  was  found  to  be  a  desirable  attribute. 
The  results  of  the  analysis  demonstrated  that  the  tension  leg  spar  buoy  has  the  potential  to  support  a 
wind  turbine  in  an  offshore  environment  without  adversely  affecting  the  loading  on  the  system 
components. 
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Barm-  length  of  moment  arm  (m) 
Ltether-  length  of  tether  (m) 

Espoke-  length  of  spoke  (m) 
m:  mass  (kg) 

M:  moment  (Nm) 

N:  generator  rotational  speed  (rpm) 
Niines'  number  of  lines  per  spoke 
P:  power  (kW) 


P:  number  of  generator  poles 
Q:  torque  (N-m) 
r:  radius  (m) 
t:  time  (s) 

S^co):  ambient  spectral  density 
T:  draft  (m) 

Tm  period  of  oscillation  (sec) 

U:  wind  velocity  (m/s) 

Ureh  relative  wind  velocity  (m/s) 

w:  spoke  characteristic  dimension  (m) 

z:  vertical  position  (m) 

a:  angle  of  attack 

a:  power  law  exponent 

p\  angle  between  nacelle  and  waves 

V :  displaced  volume  of  the  floater  (m^) 

0:  wave  potential 

0:  floater  excitation  force 

rj:  wave  elevation  (m) 

A:  tip  speed  ratio 
(o:  frequency  (rad/sec) 

(Owind:  frequency  (rad/sec) 
p:  density  (kg/m^) 

Q:  rotor  frequency  (rad/sec) 
a:  stress  (Pa) 

F:  angular  position  of  blade  (rotational) 

9t'  twist  angle 

dp:  section  pitch  angle 

Owind:  wind  incidence  angle 

dyaw'.  nacelle  yaw  angle 

S:  complex  response  amplitude 

(p:  angle  relative  to  wind 

(f:  floater  motion  (six  degrees  of  freedom) 
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1  Introduction 


f.f  Motivation 

Use  of  wind  turbines  for  electrical  power  generation  has  become  a  major  means  of  generating  clean 
renewable  energy  around  the  world.  Over  the  last  five  years  the  average  annual  growth  rate  for  the 
wind  energy  sector  has  been  35%  [12].  The  main  factors  driving  this  growth  have  been  a  desire  for 
some  countries  to  reduce  their  dependence  on  fossil  fuels  and  an  improvement  in  wind  power 
technology  that  has  made  it  more  economically  feasible.  The  European  Wind  Energy  Association 
(EWEA),  a  group  representing  European  wind  energy  companies,  has  set  a  target  of  75,000  MW 
installed  wind  power  by  2010. 

Wind  energy  systems  have  emerged  as  a  promising  technology  for  the  utilization  of  offshore  wind 
resources  for  the  large  scale  generation  of  electricity.  The  relatively  low  surface  roughness  of  the 
ocean  results  in  higher  wind  speeds.  Some  of  the  problems  associated  with  wind  power,  like  taking 
up  too  much  space  on  land,  are  lessened  when  the  machines  are  moved  offshore.  Current  projects  for 
offshore  wind  power,  like  Homs  Rev  [11],  are  all  situated  in  relatively  shallow  water  using  wind 
tiurbines  fixed  to  the  sea  floor.  The  use  of  bottom  mounted  wind  turbines  may  not  be  the  best  option 
for  deeper  offshore  locations. 

Ideally  an  offshore  wind  farm  would  be  placed  as  close  to  the  load  it  would  be  supplying  and  where 
the  winds  are  of  sufficient  velocities.  The  water  depths  of  these  locations  are  sometimes  too  deep  for 
bottom  mounted  wind  turbines  to  be  economically  feasible  (example:  Hawaii).  A  wind  farm  should 
not  interfere  with  coastal  fishing,  shipping,  and  recreational  boating.  Therefore,  mounting  wind 
turbines  on  floating  supports  may  prove  to  be  the  best  solution  to  utilize  the  wind  resources  in  areas 
with  water  depths  in  excess  of  50  meters. 

The  major  goal  of  this  research  was  to  develop  a  fully  coupled  dynamic  analysis  technique  which 
could  be  used  for  floating  wind  turbines.  The  analysis  method  had  to  incorporate  wind  loading  on  the 
rotor  blades  and  wave  loading  on  the  floater.  Also,  a  design  for  the  floating  support  was  developed  to 
achieve  the  desired  performance  characteristics.  The  fully  coupled  analysis  method  was  used  to 
evaluate  a  wind  turbine  system  mounted  on  a  tension  leg  spar  buoy  operating  in  realistic  wind  and 
wave  enviromnents. 


The  fully  coupled  dynamic  analysis  was  performed  by  combining  wind  turbine  loading  and  control 
subroutines  and  a  floater  hydrodynamic  loading  subroutine  using  a  commercially  available 
mechanical  dynamics  program.  A  floating  wind  turbine  system  was  analyzed  in  the  time  domain 
using  a  computer  model.  Damping  properties  of  the  system  were  evaluated  by  simulating  free  decay 
tests  in  the  six  modes  of  motion.  Simulations  included  normal  operations  in  various  wind  and  wave 
conditions  and  extreme  wind  and  wave  event  cases.  The  loading  responses  on  the  system 
components  for  test  cases  were  evaluated  to  detennine  the  feasibility  of  the  chosen  means  of  floating 
support. 

1.2  Previous  Studies 

Floating  wind  turbine  analysis  and  design  research  has  been  limited  in  the  past.  European  institutions 
have  been  the  leaders  in  this  field  of  research.  Most  of  the  work  has  pertained  to  possible  designs  for 
the  floating  support  system  and  economic  analyses  of  various  wind  farms  using  these  designs.  The 
following  list  is  a  partial  account  of  past  research  relating  to  floating  wind  turbines: 

•  In  the  United  Kingdom,  the  FLOAT  study  [35]  was  conducted.  This  research  looked  at 
various  designs  for  floating  wind  farms.  Cost  estimates  were  performed  on  the  designs  and 
found  to  be  approximately  twiee  that  of  bottom  mounted  wind  turbines. 

•  Research  at  University  College  London  investigated  a  design  which  mounted  several  turbines 
on  a  single  structure  with  the  potential  advantage  of  reduced  motion  response  and  shared 
moorings  [15],  [16].  The  results  of  this  work  were  that  the  design  was  prohibitively 
expensive  and  might  be  subjected  to  excessive  wave  loading. 

•  The  use  of  a  toroidal  shaped  floater  was  evaluated  by  a  group  of  Italian  researchers  [2].  This 
floater  supported  a  single  wind  turbine  in  water  depths  between  30  and  100  meters.  The 
support  was  held  in  place  by  mooring  lines  under  tension.  The  shape  of  the  floater  was 
chosen  to  minimize  wave  interaction,  but  the  design  was  deemed  too  expensive  to  build. 

•  Recently  researchers  at  Delft  University  of  Technology  in  the  Netherlands  performed  a  study 
on  using  a  floating  support  for  wind  power  applications  in  50  meters  of  water.  This  research 
was  presented  at  the  Offshore  Wind  Energy  in  Mediterranean  and  Other  European  Seas 
Conference  in  April  2003  [6]. 
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1.3  Chosen  Floater  Design 

In  the  past  many  means  for  supporting  a  floating  wind  turbine  have  been  proposed.  Some  of  the 
possible  designs  are  the  pontoon,  spar  buoy  (with  catenary  lines),  and  tension  leg  platform.  The 
requirements  for  a  floating  support  are: 

•  Must  support  its  weight  and  that  of  the  wind  turbine 

•  Must  remain  upright  for  all  design  wind  and  wave  conditions 

•  Must  prevent  excessive  loading  on  the  wind  turbine  due  to  waves  acting  on  floater 

•  Must  be  economically  feasible 

Every  design  has  advantages  and  disadvantages.  Mounting  multiple  wind  turbines  on  a  single 
pontoon  reduces  the  amount  of  power  line  which  must  be  run  along  the  ocean  floor.  However,  the 
disadvantages  are  high  structural  cost  for  the  extremely  large  pontoon  which  would  be  required  for 
adequate  spacing,  hence  high  cost,  and  laige  wave  loads  on  the  structure.  The  mounting  of  a  single 
wind  turbine  on  a  spar  buoy  is  the  simplest  and  easiest  design  to  build.  However,  with  catenary 
mooring  lines  the  system  is  not  very  stiff  in  the  pitch  and  roll  modes  of  motion  and  extremely  large 
loads  result  on  wind  turbine  components.  The  tension  leg  platform  has  very  good  performance 
characteristics  (stiff  in  pitch  and  roll).  However,  the  structure  is  relatively  complex  and  hence  the 
cost  is  likely  to  be  higher. 

A  design  used  in  this  research  combined  the  best  features  of  both  the  spar  buoy  and  the  tension  leg 
platform.  It  was  theorized  that  a  spar  buoy  which  had  been  stiffened  in  the  roll  and  pitch  modes  by 
using  a  combination  of  spokes  and  lines  under  tension  could  meet  the  requirements  for  a  floating 
support.  The  proposed  design  was  called  the  tension  leg  spar  buoy  (see  Figure  1.1),  The  system 
would  be  set  up  such  that  operators  could  board  the  floating  wind  turbine  via  a  walkway  which  circles 
the  base  of  the  tower.  The  inset  in  the  figure  shows  the  system  attached  to  gravity  anchors  on  the 
seafloor. 


Figure  1.1  Floater  design. 


Figure  1.2  is  a  cross  sectional  view  of  the  floater  structure.  The  design  consists  of  a  tower  extension 
that  connects  the  wind  turbine  tower  to  the  spar  cylinder.  The  orange  portions  represent  the  part  of 
the  structure  that  is  normally  under  water.  The  relatively  small  diameter  of  this  structure  near  the  free 
surface  reduces  wave  loading  on  the  floater.  A  majority  of  the  buoyancy  comes  from  the  spar 
cylinder.  The  figure  shows  the  ballast  chamber  which  has  fill  and  drain  pipes  which  run  up  the 
structure  to  the  catwalk  where  they  can  be  connected  to  pumps  onboard  a  support  ship.  Installation  of 
the  system  is  accomplished  by  filling  the  ballast  chamber  and  attaching  the  tethers.  Then  the  ballast 
is  pumped  off  and  tension  put  on  the  lines. 
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Figure  1.2  Cutaway  of  floater. 

The  wind  turbine  tower  is  attached  to  a  vertical  cylinder  which  provides  a  majority  of  die  buoyancy 
for  the  system.  The  spokes  are  attached  to  tethers  which  connect  the  system  to  either  anchors  driven 
into  the  seafloor  or  gravity  anchore.  The  design  has  relatively  low  restoring  coefficients  in  the 
horizontal  translational  modes  of  motion. 

The  computer  model  used  in  this  research  models  the  spokes  as  rectangular  members  with  square 
cross-sections.  These  spokes  are  reinforced  to  increase  their  stiffness  in  bending.  The  spokes  could 
be  made  buoyant  to  reduce  the  required  size  of  the  spar  cylinder.  However,  the  use  of  a  truss 
structure  might  be  the  most  economical  means  to  fulfill  the  function  of  the  spokes  which  is  to  increase 
the  stiffness  of  the  system  in  all  of  the  rotational  modes  of  motion.  One  of  the  benefits  of  using  a 
tension  leg  spar  buoy  is  that  the  pitch  and  roll  motions  are  lower  than  other  types  of  floating  support. 
As  shown  in  Figure  1.3,  the  structure  pitches  very  little  when  the  tension  leg  spar  buoy  undergoes 
surge  motion.  The  pitch  motion  is  resisted  by  the  tethers.  Having  pre-tension  in  the  lines  allows  the 
system  to  absorb  pitch  /  roll  moments  caused  by  the  wind  and  waves.  This  rotational  stiffness  is 
determined  by  the  spring  coefficient  provided  by  the  lines  and  the  length  of  the  spokes. 


Figure  1.3  Diagram  of  system  motion  with  horizontal  excitation. 

The  floater  is  assumed  to  be  a  rigid  body  connected  to  the  seafloor  by  flexible  lines.  The  tethers  only 
provide  rotational  stiffness  so  long  as  they  maintain  tension.  Therefore,  the  amount  of  reserve 
buoyancy  is  a  key  design  parameter  which  must  be  high  enough  to  ensure  that  the  lines  remain  under 
tension  to  prevent  excessive  amounts  of  pitch  and  roll  which  have  a  large  effect  on  wind  turbine 
loading. 


1.4  Outline  of  Thesis 

Many  equations  used  in  the  modeling  of  the  fully  coupled  analysis  technique  are  described  in  this 
report.  However,  not  all  of  the  basic  theory  is  covered  in  depth.  Where  necessary  the  applicable 
sources  for  this  theory  are  given.  The  following  is  an  outline  of  this  thesis. 

•  Chapter  2  provides  a  basic  description  of  the  wind  turbine.  The  primary  components  of  the' 
machine  are  described.  This  chapter  also  includes  the  wind  turbine  outputs  and  a  general 
description  of  the  wind  environment. 

•  Chapter  3  describes  the  problem  of  modeling  the  fully  coupled  floating  wind  turbine  system. 
The  different  forces  and  moments  acting  on  the  floater  and  turbine  rotor  are  described. 
Finally,  the  different  loading  conditions  used  in  evaluating  the  floater  and  analysis  technique 
are  presented. 

•  Chapter  4  is  a  brief  description  of  the  hydrostatic  properties  of  the  system.  The  effect  of 
static  loading  on  system  deflection  and  the  restoring  coefficients  are  discussed. 

•  Chapter  5  details  the  hydrodynamic  analysis  of  the  system.  The  six  modes  of  motion  and 
their  associated  equations  are  detailed. 

•  Chapter  6  deals  with  hydrodynamic  loading  on  the  floater.  Both  the  inertial  and  drag  loading 
equations  are  described.  The  production  of  random  wave  records  is  discussed. 

•  Chapter  7  describes  the  means  by  which  the  floating  wind  turbine  system  model  is  generated 
and  analyzed.  The  various  computer  codes  which  are  used  to  conduct  the  fully  coupled 
dynamic  simulations  are  discussed.  The  wind  turbine  model  is  also  detailed. 

•  Chapter  8  presents  the  results  of  the  damping  analysis  of  the  system.  Each  of  the  individual 
modes  of  motion  is  discussed.  The  final  portion  of  the  chapter  gives  the  results  of  the 
multiple  mode  damping  tests. 

•  Chapter  9  details  the  results  of  the  numerical  simulations  of  the  system  model  for  normal  and 
extreme  operating  environments.  The  loading  responses  of  the  wind  turbine  components  and 
the  lines  are  discussed. 


2  Wind  Turbine 


A  wind  turbine  is  a  device  that  coverts  the  kinetic  energy  in  wind  into  useful  work  output.  The  output 
can  be  in  the  form  of  electrical  power.  In  the  last  decade  there  has  been  an  increased  interest  in  using 
wind  turbines  to  produce  electrical  power.  This  chapter  is  included  to  give  the  reader  a  basic 
understanding  of  wind  turbines.  A  full  description  of  how  a  wind  turbine  operates  can  be  found  in 
“Wind  Energy  Explained”  by  Manwell,  McGowan,  and  Rodgers  [24]. 

2.1  Components 

The  most  common  wind  turbine  design  is  the  Horizontal  Axis  Wind  Turbines  (HAWT).  The  output 
shaft  of  the  HAWT  is  parallel  to  the  ground.  A  wind  turbine  consists  of  seven  major  subsystems. 
These  subsystems  are  the:  rotor,  drive  train,  nacelle,  generator,  tower  and  foundation,  control  system, 
and  electrical  output  system. 

2.1.1  Rotor 

The  rotor  consists  of  normally  two  or  three  blades  attached  to  a  hub.  The  system  performance  of  the 
wind  turbine  is  based  on  the  selection  of  blade  number,  shape,  and  length.  The  rotor  can  be  either 
upwind  or  downwind  design.  Most  wind  turbines  currently  in  use  are  three  bladed  upwind  designs. 
The  blades  can  be  designed  with  either  constant  pitch  or  controllable  pitch.  The  power  of  the  rotor  is 
controlled  via  aerodynamic  control  (stall  regulation)  or  variable  pitch  blades. 

2.1.2  Drive  Train 

Wind  passing  across  the  rotor  blades  produces  lift  and  drag.  The  lift  force  generates  a  torque  about 
the  rotational  axis  of  the  rotor.  The  output  shaft  of  the  wind  turbine  rotates  at  a  relatively  slow  speed. 
For  most  electric  generators,  the  rotor  output  must  be  stepped  to  a  higher  speed.  A  gearbox  is  used  to 
step  up  the  rotor  output  to  spin  the  generator.  The  rest  of  the  drive  train  consists  of  journal  and  thrust 
bearings  which  support  the  shafts  and  minimize  movement  of  the  drive  drain  in  the  radial  and  axial 
directions.  The  gearbox  in  wind  turbines  is  one  of  the  major  sources  of  failure  for  the  system.  For 
this  reason  some  wind  turbines  are  designed  to  have  generators  that  operate  at  the  speed  of  the  rotor. 


2.1.3  Nacelle 


The  nacelle  is  the  structure  that  houses  the  generator  and  drive  train.  This  structure  protects  the  wind 
turbine  equipment  from  the  elements.  Offshore  wind  turbines  are  different  from  land  based  wind 
turbines  in  that  the  nacelle  is  sealed  to  prevent  entry  of  sea  air.  A  cooling  system  is  installed  to 
transfer  the  heat  generated  by  the  wind  turbine  machinery  to  the  ocean  air  passing  over  the  nacelle. 

The  alignment  of  the  nacelle  is  controlled  by  the  yaw  system.  This  system  connects  the  nacelle  to  the 
top  of  the  tower.  The  yaw  direction  of  all  upwind  wind  turbines  must  be  actively  controlled  so  that 
the  rotor  faces  into  the  wind.  Some  downwind  turbines  use  free  yaw  systems,  which  align  themselves 
with  the  wind  automatically.  Yaw  of  a  wind  turbine  generates  a  gyroscopic  moment  proportional  to 
the  rate  of  yaw  and  the  rotational  inertia  of  the  rotor.  Therefore,  rate  of  yaw  of  the  nacelle  must  be 
limited.  Active  yaw  control  is  usually  accomplished  using  a  pinion  and  bull  gear  arrangement. 

2.1.4  Generator 

The  generator  converts  the  mechanical  work  input  of  the  wind  turbine  into  useful  electrical  output. 
Most  wind  turbines  use  either  induction  or  synchronous  generators  to  accomplish  this  tesk.  The 
majority  of  wind  turbines  use  induction  generators.  The  rotor  of  an  induction  generator  spins  slightly 
faster  than  the  rotating  magnetic  field  on  the  stator.  The  difference  between  the  generator  rotor 
frequency  and  the  rotating  magnetic  field  on  the  generator  is  called  slip.  The  slip  of  an  induction 
generator  determines  flie  amount  of  power  produced  by  the  machine, 

2.1. 5  Tower  and  Foundation 

The  wind  turbine  luu^elle  and  rotor  must  be  supported  at  their  operational  height.  The  height  of  the 
supporting  tower  can  be  over  100  meters  high.  The  height  of  the  tower  is  determined  by  the  rotor 
diameter  and  the  wind  conditions  at  the  site.  The  speed  of  wind  increases  with  height  above  ground 
due  to  wind  shear.  Higher  wind  speeds  result  in  increased  power  output.  It  is  therefore  desirable  to 
raise  the  wind  turbine  nacelle  up  off  the  ground  to  take  advantage  of  this  effect. 

The  stiffness  of  the  support  tower  must  be  chosen  such  that  there  is  ho  dynamic  coupling  between  the 
rotor  and  the  tower.  Coupled  vibrations  lead  to  reduced  fatigue  life  of  the  machine.  Downwind 
turbines  also  have  to  deal  with  vortices  shed  off  of  the  tower,  called  tower  shadow.  As  the  rotor 
blades  pass  through  the  tower  shadow  noise  and  oscillations  are  produced.  This  results  in  increased 
noise  and  vibrations  of  the  rotor  blades. 


2.1.6  Control  System 


The  control  system  changes  the  blade  pitch,  nacelle  yaw,  and  generator  loading  of  a  wind  turbine. 
The  yaw  is  changed  to  keep  the  nacelle  aligned  with  the  wind  in  order  to  maximize  power  output.  If 
torque  on  the  rotor  becomes  too  high  the  control  system  can  also  change  the  yaw  such  that  the  rotor 
presents  a  smaller  frontal  area  to  the  wind,  and  thereby  reduces  the  amount  of  torque  produced.  The 
control  system  can  also  change  the  pitch  of  the  blades  to  alter  the  amount  of  torque  produced  by  the 
rotor.  The  purpose  of  the  control  system  is  to  maximize  power  output  and  fatigue  life.  Manufactures 
use  various  means  to  accomplish  this  purpose. 

2.1.7  Electrical  Output  System 

Electrical  output  of  the  generator  must  be  transferred  to  the  electric  grid  or  electric  load.  The  output 
passes  though  other  electrical  components  like:  cables,  transformers,  power  factor  correction 
capacitors,  solid  state  power  converts,  and  switchgear.  Some  wind  turbine  designs  use  solid  state 
power  converts  for  increased  flexibility.  Normally  the  output  voltage  is  stepped  up  to  34.5  kilovolts. 
Offshore  wind  turbines  typically  send  power  through  cables  buried  on  the  seafloor. 

2.2  Wind  Turbine  Output 

Wind  passing  over  the  rotor  blades  generates  torque  and  thrust.  The  torque  is  then  converted  into 
electrical  power  output.  Rotor  torque  increases  with  wind  speed.  In  order  for  the  generator  to 
produce  useful  power  output  the  wind  speed  must  exceed  the  cut-in  speed  of  the  generator.  Wind 
speed  increases  to  the  rated  power  output  of  the  wind  turbine.  Above  the  rated  speed  wind  turbine 
power  is  controlled  by  changing  the  blade  pitch  or  by  stalling  a  portion  of  the  blade.  The  speed  above 
which  the  wind  turbine  can  no  longer  operate  is  called  the  cutout  speed.  Figure  2.1  is  a  graph 
showing  the  power  versus  wind  speed  for  a  500  kilowatt  wind  turbine.  This  graph  was  generated 
using  a  program  called  PROPID  which  is  available  on  the  National  Wind  Technology  Center  website. 


2.3  Wind  Environment 

One  of  the  benefits  of  building  wind  turbines  offshore  is  the  low  surface  roughness  Of  the  sea.  Wind 
passing  of  obstacles  such  as  trees,  houses  and  hills  acts  to  lower  the  wind  speeds  closer  to  the  surface. 
Wind  speeds  increase  with  height  above  the  surface.  The  wind  profile  is  normally  modeled  using 
either  the  power  law  profile  or  the  logarithmic  profile  [29]. 


Power  Law  Profile:  F(z)  =  F(z  J  • 


Logarithmic  Profile:  F(z)  =  F(z 

ln(z^/zj 


(2.1) 


where: 

F(z)  is  the  wind  speed  at  height  z 
V(zr)  is  the  reference  wind  speed 
z  is  the  height  above  ground 

Zr  is  the  reference  height  above  ground  used  for  fitting  the  profile 
Zo  is  the  roughness  length 
a  is  the  power  law  exponent 
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A  draft  of  lEC  61400-3  Safety  Requirements  for  Offshore  Wind  Turbines  [20]  states  that  for  offshore 
environments  the  power  law  exponent  should  be  assumed  to  be  0.14.  A  value  of  0.2  is  typically  used 
for  normal  wind  environments  on  land.  When  the  logarithmic  profile  is  used  for  calm  seas,  the 
surface  roughness  length  (z„)  is  0.50  millimeters  (Table  2.2  of  [24]).  The  wind  profile  out  at  sea  is 
such  that  wind  speeds  are  typically  higher  than  on  land  for  the  same  height  above  the  surface.  The 
difference  in  wind  profiles  improves  the  operation  of  the  wind  turbine.  An  offshore  wind  turbine 
rotor  can  be  mounted  on  a  shorter  tower  and  still  produce  the  same  amount  of  power  output  as  a  land- 
based  counterpart. 

The  wind  speed  increases  with  height  above  the  surface  of  the  water.  As  the  rotor  spins,  the  blades 
pass  through  the  wind  profile.  This  results  in  a  changing  inflow  speed  to  each  blade  that  varies  with 
time  in  a  periodic  manner.  Oscillating  thrusts  and  moments  are  produced  by  the  blades  with  a 
frequency  corresponding  to  the  number  of  blades  (N)  times  the  rotational  speed  of  the  rotor. 
Modeling  the  wind  in  a  turbulent  manner  is  important  since  actual  winds  are  never  constant.  Figure 
2.2  is  a  wind  record  for  a  mean  wind  speed  of  20  meters  /  second. 


3  Problem  Statement 


Several  problems  must  be  solved  in  order  to  determine  the  feasibility  of  a  tension  leg  spar  buoy 
support  for  an  offshore  wind  turbine.  One  part  of  the  analysis  is  to  determine  the  effect  of  floater 
motion  on  wind  turbine  performance  and  fatigue  life.  Another  important  part  of  the  analysis  is  to 
determine  the  loading  on  the  tethers  that  attach  the  floater  to  the  ocean  floor. 

The  wind  turbine  /  floater  system  are  acted  upon  by  both  waves  and  wind.  The  combination  of  the 
various  components  leads  to  the  analysis  of  a  coupled  dynamic  system.  When  wave  loads  cause  the 
floater  to  move,  the  nacelle  moves  with  it.  The  motion  of  the  nacelle  results  in  changes  to  the  inflow 
wind  to  die  turbine  rotor.  Changes  in  the  inflow  to  the  rotor  can  result  in  changes  in  rotor  thrust 
loading.  This  change  in  thrust  is  felt  by  the  floating  support.  Angular  motion  of  the  floater  can  also 
result  in  g5n'oscopic  loads  on  the  rotor. 

The  wind  turbine  itself  is  a  very  complex  system  to  model  accurately.  Several  wind  turbine 
simulation  codes  can  be  used  to  accomplish  this  task.  For  this  research  the  wind  turbine  is  modeled 
using  the  FAST  wind  turbine  code  [4]  in  conjunction  with  ADAMS  [1].  FAST  (Fatigue, 
Aerodynamics,  Structures,  and  Turbulence)  code  is  produced  by  the  National  Renewable  Energy 
Laboratory.  ADAMS  (Automatic  Dynamic  Analysis  of  Mechanical  Systems)  is  commercially 
available  software  that  is  widely  used  to  model  mechanical  systems.  Codes  modeling  the  floater 
loading  and  mechanical  properties  are  coupled  to  the  ADAMS  wind  turbine  model  to  determine  the 
effect  placing  it  on  a  floating  platform  had  on  the  wind  turbine  components.  The  system  is  modeled 
as  a  group  of  71  rigid  bodies  (called  parts)  connected  by  springs  and  dampers.  For  example  each 
rotor  blade  consists  of  15  different  parts,  one  for  each  blade  element. 

Loads  due  to  wind  and  waves  on  the  individual  parts  are  calculated  relative  to  each  part’s  body  based 
coordinate  system.  All  of  the  forces  and  moments  from  the  individual  components  of  the  system  are 
incorporated  into  a  complete  dynamics  analysis  relative  to  the  global  coordinate  system.  The  motion 
of  every  part  is  determined  by  evaluating  of  the  combined  dynamic  equations  for  the  system.  These 
calculations  are  carried  out  in  the  time  domain  using  the  ADAMS  Solver.  The  analysis  tools  and  the 
system  model  are  discussed  in  Chapter  7. 


3. 1  Modes  of  Floater  Motion 

The  combined  rigid  body  of  the  wind  turbine  and  floating  support  undergoes  oscillatory  translation 
and  rotation  due  to  wind  and  wave  loadings  (see  Figure  3.1).  Surge  and  sway  are  horizontal  motions 
along  the  x  and  y  axes,  respectively.  The  horizontal  axes  of  the  system  are  aligned  with  the  four 
spokes  to  which  the  floaters  tethers  are  connected.  Heave  is  the  vertical  motion  of  the  floater.  The 
oscillatory  angular  motions  are  referred  to  as  roll,  pitch  and  yaw.  Surge,  sway,  and  yaw  are  the  three 
rigid  body  motions.  The  tether  lengths  are  assumed  to  remain  constant  during  rigid  body  motion. 
Pitch,  roll  and  heave  occur  when  the  tethers  undergo  elastic  elongation  /  contraction.  Heave  also 
occurs  when  the  structure  undergoes  surge  or  sway  motion  (like  an  inverted  pendulum).  The  origins 
of  the  three  axes  are  located  at  the  center  of  the  floater  on  the  calm  waterline.  The  reference  for  the 
global  coordinate  system  is  a  point  at  the  center  of  the  tower  where  it  connects  to  the  floater  (on  the 
free  surface).  The  orange  portions  of  the  structure  are  submerged  when  the  system  is  at  rest. 
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Figure  3.1  Floater  modes  of  motion. 
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3.2  Forces  and  Moments 

Forces  and  moments  due  to  gravity,  wind,  and  waves  are  exerted  on  both  the  wind  turbine  assembly 
and  the  floater.  Wind  forces  and  moments  are  primarily  exerted  on  the  rotor  of  the  wind  turbine, 
while  wave  forces  ^t  only  on  the  floater.  The  forces  and  moments  produced  on  the  wind  turbine 
blades  are  transferred  to  the  floater.  The  following  is  a  discussion  on  the  forces  and  moments  acting 
on  the  combined  system. 

3.2.1  Wind  Turbine 

Wind  passing  across  the  rotor  blades  produces  lift  and  drag  forces.  The  two  forces  result  in  a  torque 
on  the  rotor,  a  thrust  along  the  rotor  axis,  and  a  pitching  moment.  The  thrust  on  the  blades  is  also  felt 
by  the  tower  and  the  floater.  Figure  3.2  is  a  diagram  of  forces  and  moment  on  a  stationary  airfoil  in  a 
wind  stream.  For  ideal  non-separated  flow,  these  forces  and  moment  act  along  the  cord  at  a  distance 
of  approximately  c/4  from  the  leading  edge  [24],  where  c  is  the  cord  length. 

Lift  Force 


Figure  3.2  Forces  and  moments  on  a  stationary  airfoil. 


The  lift  force  is  defined  as  the  force  component  that  is  perpendicular  to  the  incoming  air  flow.  The 
drag  force  is  defined  to  be  parallel  to  the  direction  of  the  incoming  airflow.  Drag  arises  fi'om  vortex 
shedding,  viscous  friction,  and  separated  wakes.  The  pitching  moment  is  the  summation  of  the 
moments  created  by  the  individual  forces  acting  on  the  airfoil.  This  moment  acts  about  an  axis  that  is 
perpendicular  to  the  cross-section  of  the  airfoil.  On  a  wind  turbine  the  rotation  of  tire  rotor  increases 
the  inflow  into  the  airfoil,  while  at  the  same  time  changing  the  effective  angle  of  attack. 


The  thrust  on  individual  blades  causes  a  moment  at  the  blade  root  which  is  perpendicular  to  the 
pitching  moment.  When  the  nacelle  is  aligned  as  in  Figure  3.1  these  moments  act  about  the  y  and  z 
axes  (rotor  torque  acts  around  the  x-axis.)  The  orientations  of  the  blades  change  as  they  rotate. 
Therefore,  cyclic  moments  are  produced  with  frequencies  equal  to  the  blade  frequency.  These 
moments  are  excitation  for  the  floater  in  both  yaw  and  pitch.  It  should  be  noted  that  the  effects  of 
these  moments  will  be  small  relative  to  the  rotor  torque  and  moment  due  to  rotor  thrust  at  the  base  of 
the  tower.  The  reasons  for  this  will  be  shown  in  the  analysis  section. 

Wind  speed  increases  with  elevation.  Oscillating  forces  and  moments  are  produced  as  the  rotor 
blades  pass  through  this  varying  wind  field.  These  cyclic  loads  can  affect  the  fatigue  life  of  large 
diameter  wind  turbines.  The  magnitudes  of  these  oscillating  loads  on  the  floater  are  small  compared 
to  the  forces  and  moments  generated  by  the  rotor  torque  and  the  rotor  thrust  acting  with  the  moment 
arm  of  the  tower.  The  oscillating  forces  may  excite  resonant  motion  of  the  floater  and  wind  turbine 
and  must  therefore  be  considered. 

Gyroscopic  motion  of  a  wind  turbine  is  a  factor  in  turbine  fatigue.  In  land  based  turbines  this  results 
primarily  from  yawing  of  the  wind  turbine.  When  the  wind  turbine  is  mounted  on  a  floating  structure, 
both  yaw  and  pitch  of  the  floater  cause  gyroscopic  moments.  When  the  plane  of  rotation  of  the  wind 
turbine  rotor  is  changed,  a  precession  torque  is  developed.  This  precession  torque  acts  about  an  axis 
perpendicular  to  both  the  rotor  axis  and  the  axis  about  which  the  plane  of  rotation  is  being  changed. 
Pitch  motion  will  generate  a  gyroscopic  moment  while  roll  motion  will  not  produce  a  gyroscopic 
moment  since  the  plane  of  rotor  rotation  does  not  change  during  roll  when  the  wind  turbine  is  aligned 
as  in  Figure  3.1.  The  gyroscopic  moment  produced  is  the  cross  product  of  the  angular 

momentum  vector  (  and  the  angular  velocity  vector  (<^).  Where  the  angular  velocity  is  the  first 
derivative  of  the  pitch,  roll,  or  yaw  motions  [24,  p.  148]: 

(3.1) 

where  J  is  the  polar  moment  of  inertia  of  the  rotor  and  Q  is  the  angular  velocity  vector  of  the  rotor. 

3.2.2  Floater  Wave  Loading 

As  shown  in  Figure  3.1,  the  tethered  spar  buoy  consists  of  a  vertical  cylinder  that  provides  buoyancy 
and  four  spokes  to  which  the  tethers  are  attached.  The  floater  is  symmetrical  about  the  z-axis. 
Therefore,  only  heave,  surge,  and  pitch  are  considered  in  the  hydrodynamic  analysis.  There  is  veiy 
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little  wave  excitation  in  the  yaw  mode,  since  the  major  portions  of  the  floater  are  vertical  cylinders. 
However,  forces  and  moments  generated  on  the  wind  turbine  may  excite  the  other  modes  of  floater 
motion.  The  following  discussion  deals  with  surge,  pitch,  and  heave  excitation.  Sway  and  roll  are 
similar  to  surge  and  pitch  since  all  these  modes  are  caused  by  horizontal  forces  on  the  submerged 
structure.  Therefore,  the  discussion  for  surge  and  pitch  are  applicable  to  sway  and  roll. 

For  a  typical  spar  buoy,  the  pitch  and  surge  motion  equations  are  coupled.  The  system  is  not  very 
stiff  in  the  surge  direction.  The  stiffness  of  the  system  in  both  heave  and  pitch  will  result  in  very 
small  motions  of  the  floater  in  these  two  modes.  Since  the  system  acts  like  an  inverted  pendulum, 
some  of  the  motion  in  flie  heave  mode  is  due  to  the  set-down  caused  by  surge.  As  the  floater  moves 
away  from  the  rest  position  its  draft  increases.  This  coupling  between  surge  /  sway  and  heave  is  non¬ 
linear  and  determined  by  the  length  of  the  tethers.  Shorter  tethers  have  more  coupling,  since  set- 
down  is  large  for  the  same  amount  of  horizontal  motion. 

The  motions  of  the  floater  may  have  an  effect  on  the  wind  turbine.  Wave  loading  is  assumed  to  be 
produced  by  regular  sinusoidal  waves.  The  motion  of  a  floating  structure  in  irregular  waves  can  be 
determined  by  linearly  superimposing  the  results  from  multiple  regular  waves  to  produce  realistic 
seastates.  This  method  will  discussed  later  in  this  thesis  and  is  fully  described  in  Faltinsen,  1990  [13] 
and  Sclavounos,  2003  [31]. 

The  analysis  of  the  tethered  spar  buoy  uses  the  long  wavelength  assumption  that  applies  when  the 
wavelength  is  much  longer  than  the  spar  diameter  (X  >  5  d).  This  assumption  results  in  no  significant 
wave  disturbance  being  generated  by  the  spar  buoy.  The  excitation  forces  are  caused  by  both 
incoming  wave  and  diffraction  potential  flow  effects.  The  excitation  forces  and  moments  are 
obtained  from  the  incoming  wave  potential  and  using  analytical  expressions  for  the  added  mass  terms. 
The  force  caused  by  the  incoming  wave  is  approximated  using  the  long-wavelength  GI  Taylor 
equations  [31].  The  equations  of  motion  for  the  spar  buoy  are  presented  in  the  analysis  section. 

The  behavior  of  a  floating  wind  turbine  is  complex.  Modes  of  motion  interact  to  produce  complex 
behaviors.  For  example,  normal  operation  of  a  floating  wind  turbine  will  produce  a  yaw  moment  on 
the  floater  two  ways.  The  first  yaw  moment  is  due  to  the  oscillating  yaw  produced  by  the  normal 
operation  of  the  rotor  blades  spiiming  which  results  in  the  yaw  moment  changing  with  time.  The 
second  yaw  moment  is  produced  due  to  the  fact  that  a  floating  wind  turbine  tilts  slightly  while  under 
load.  With  both  wind  and  waves  aligned  with  the  x-axis,  a  roll  moment  is  exerted  on  the  floater  by 
the  generator  torque.  The  tethers  which  lie  on  the  y-axis  (lines  2  and  4)  resist  this  motion.  The 


moment  causes  a  small  roll  deflection  of  the  floater  due  to  the  elongation  of  one  tether  and  the 
contraction  of  another  tether  on  the  opposite  side  of  the  floater.  This  roll  deflection  is  small  due  to 
the  relatively  high  modulus  of  elasticity  of  the  tethers.  Roll  of  the  floater  results  in  the  nacelle 
moving  off  axis.  This  movement  produces  a  moment  arm  along  the  y-axis  that  is  acted  on  by  the 
rotor  thrust.  A  yaw  moment  is  produced.  It  should  be  stressed  that  in  most  cases  this  moment  is 
small,  since  the  stiffness  of  the  system  in  roll  is  high.  Figure  3.3  shows  how  roll  is  generated  by 
torque  on  the  generator  (amount  of  roll  is  exaggerated).  With  wind  into  the  rotor,  you  can  see  how  a 
moment  is  produced  about  the  z-axis. 


3.3  Loading  Conditions  used  for  Analysis 

For  most  of  the  following  discussion  wind  is  assumed  to  blow  in  the  positive  x-direction  and  the  rotor 
is  aligned  along  the  x-axis.  This  assumption  is  only  for  the  purposes  of  this  discussion.  The  actual 
calculations  take  into  account  variations  in  the  wind  and  wave  directions.  However,  any  angular 
difference  between  wind  and  waves  can  be  accounted  for  using  the  wave  incidence  angle,  p.  Five 
different  loading  conditions  are  applicable  to  the  analysis  presented  in  this  thesis. 
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•  Stetic  loading  with  no  wind  or  waves 

•  Static  loading  with  constant  wind. 

•  Static  and  dynamic  loading  with  waves  acting  on  the  floater  and  no  wind. 

•  Static  and  dynamic  loading  for  random  wind  and/or  waves  on  the  system. 

•  Extreme  wind  or  wave  loading  events. 

The  static  loading  condition  accounts  for  forces  and  moments  which  arise  from  gravity  and  buoyancy. 
The  floater  supports  the  combined  weight  of  the  system.  The  floater  has  excess  buoyancy  (buoyant 
force  >  system  weight),  which  puts  tension  on  the  four  tethers.  The  buoyancy  is  dependent  on  the 
submerged  volume  of  the  floater.  For  a  simple  spar  buoy,  the  total  volume  of  seawater  displaced  by 
the  floater  (  V )  is  a  function  of  the  floater  diameter  (d)  and  floater  dmft  (7): 

V  =  -d^T  (3.2) 

4 

This  displacement  results  in  an  associated  buoyant  force  that  is  greater  than  the  total  weight  being 
supported.  The  excess  buoyancy  produces  tension  in  the  four  tethers.  The  following  equation 
presents  how  to  calculate  the  static  tether  force: 

^tethers  ‘  Psw  ^total )  S  (3-3) 

where  is  the  total  mass  of  the  system. 

The  second  loading  condition  is  steady  state  winds  with  no  waves.  This  loading  condition  is  very 
useful  for  the  damping  coefficients  for  the  system.  Steady  state  winds  result  in  nearly  constant  thrust 
and  torque  being  exerted  on  the  rotor.  Movement  of  the  floater  results  in  changes  to  flie  wind  inflow 
to  the  rotor.  This  change  in  wind  inflow  alters  the  thrust  and  torque  loading  on  the  rotor  results.  The 
wind  turbine  acts  as  a  damper  for  some  modes  of  system  motion. 

The  third  loading  condition  is  applicable  when  the  turbine  is  not  operating  and  waves  are  acting  on 
the  floater.  Waves  acting  on  the  floater  produce  excitation  forces  and  moments  in  the  surge,  sway, 
roll,  pitch  and  heave  modes  of  motion.  Since  the  floater  is  a  right  circular  cylinder,  yaw  excitation  is 
generated  only  by  waves  acting  on  the  spokes.  This  condition  is  used  to  model  shutdown  opemtions. 
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The  next  loading  condition  has  both  random  wind  and  waves  acting  on  the  system.  This  condition  is 
useful  for  determining  expected  system  response  during  normal  operations  for  different  reference 
wind  speeds  and  sea  states. 

The  response  of  the  system  to  extreme  wind  and  wave  events  is  a  key  design  feature  which  has  to  be 
understood  if  a  particular  system  is  being  analyzed  for  suitability.  Wind  and  wave  inputs  modeling 
extreme  events  are  used  in  to  the  numerical  simulation.  The  fully  coupled  dynamic  response  is 
evaluated  and  the  resulting  component  loading  values  are  determined. 


4  Static  Analysis  of  Floater 


The  physical  characteristics  and  loading  conditions  of  the  combined  floater  /  wind  turbine  system 
must  be  underctood  in  order  to  develop  a  computer  model  that  accurately  predicts  the  physical 
response.  The  purpose  of  this  chapter  is  to  describe  the  analytical  methods  used  to  solve  the  problems 
presented  in  the  previous  chapter.  First,  the  static  analysis  of  the  system  is  presented.  The  means  by 
which  static  deflection  can  be  obtained  by  calculating  the  floater  restoring  coefficients  are  described. 

An  important  design  parameter  is  tether  line  stress.  If  line  stress  exceeds  the  ultimate  stress  for  the 
tethers  catestrophic  failure  of  the  system  could  result.  Excess  buoyancy  of  the  spar  is  used  to  put 
tension  on  the  tethem.  The  tension  present  in  the  tethers  while  the  system  is  at  rest  is  referred  to  as 
the  pre-tension.  Excitation  forces  and  moments  by  waves  and  wind  on  the  system  cause  the  line 
stress  to  vary  about  this  tether  pre-tension.  The  following  equation  is  used  to  calculate  the  pre-tension 
stress  of  tethers  with  diameter  {d,e,her)  for  a  given  number  of  tethers  per  spoke  (Ntethen)'. 

fy  _  ‘  ~  ^’otot )  S  _ ^tethers _  tA  It 

4'  N  L  •  —  d  L  ^  ^tethers  ’  ^  '  ^tether 

^  ^  tethers  tether 

External  forces  on  the  system  cause  increases  and  decreases  in  the  tether  line  stresses.  As  d^cribed 
above  a  moment  will  cause  lines  to  increase  and  decrease  in  tension.  Changes  in  vertical  forces  are 
shared  by  the  four  tethers.  As  the  floater  moves  line  stresses  change.  It  is  important  to  know  how  the 
line  stress  in  each  of  the  four  tethers  will  change  during  operations. 

In  steady  state  operation,  static  loading  of  the  wind  turbine  produces  a  horizontal  force  along  the  x- 
axis  called  thrust.  This  force  causes  a  horizontal  displacement  of  the  floater  in  the  x-direction  (see 
Figure  1.3).  This  thrust  is  resisted  by  a  horizontal  restoring  force  produced  by  the  tethers  as  their 
fairlead  position  is  moved  from  the  rest  position  which  is  directly  above  the  anchor  point  for  each 
tether.  As  the  floater  surges  away  from  the  rest  position  its  draft  increases,  due  to  the  fact  that  the 
tethers  act  as  rods  in  tension.  This  increase  in  draft  would  cause  an  increase  in  the  line  tension. 
However,  for  small  surge  displacements  (^/  «  water  depth,  D)  a  linearization  assumption  can  be 
made  that  the  line  tension  remains  constant.  This  assumption  is  based  on  the  fact  that  the  tethers  will 
have  a  very  large  tension  stress  while  at  rest.  The  draft  of  the  buoy  is  also  assumed  not  to  increase  for 
small  suige  motions.  To  the  leading  order,  the  surge  and  sway  restoring  coefficients  are  identical  due 
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to  symmetry.  This  restoring  coefficient  can  also  be  obtained  by  setting  the  moment  about  the  tethers 
anchors  on  the  sea  bottom  to  zero  and  solving  for  the  force  which  produces  an  opposite  moment  to 
the  moment  due  to  excess  buoyancy  for  a  given  displacement.  For  small  displacements  (small  angle 
approximation),  the  surge  and  sway  restoring  coefficients  are  given  by  the  following  equation: 

C„=C^=^  (4.2) 

^tether 

where  L,e,her  is  the  water  depth  minus  the  depth  at  which  the  lines  connect  to  the  buoy.  Horizontal 
displacement  for  a  given  loading  can  then  be  found  by  dividing  the  force  applied  to  the  floater  by  the 
restoration  coefficient  as  described  by  expression  (4.3).  These  restoration  coefficients  are  also  used 
when  estimating  the  natural  periods  in  surge  and  sway. 


(4.3) 


Heave  restoring  forces  results  from  both  tether  and  hydrostatic  effects.  If  the  structure  moves  only  in 
the  heave  direction  the  tethers  extend  and  contract.  For  a  system  using  stiff  lines  for  tethers,  this 
restoring  mechanism  is  much  greater  than  the  hydrostatic  restoring  mechanism.  The  heave  restoring 
coefficient  is  given  by  the  following  equation: 
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(4.4) 


where  A,e,hers  is  the  total  cross  sectional  area  of  the  tethers  attached  to  a  spoke.  The  vertical 
displacement  from  a  vertical  load  can  be  found  by: 

^3=7^  (4.5) 

'^33 

Wave  loading  on  the  submerged  stmcture  creates  pitch  and  roll  moments  on  the  system.  Thrust  on 
the  rotor  acting  along  with  the  moment  arm  of  the  tower  creates  a  moment  about  the  y-axis,  which 
also  results  in  pitch  of  the  floater.  Pitch  and  roll  restoring  coefficients  are  similar  due  to  symmetry. 
Pitch  and  roll  restoring  coefficients  are  calculated  by  taking  into  account  the  moments  produced  by 
the  lines  due  to  extension  and  contraction,  the  moment  due  to  line  tension  {Fathers),  the  moment  due  to 
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system  mass,  and  the  moment  due  to  system  buoyancy  (see  Figure  4.1).  For  stiff  lines  the  largest 
restoring  moment  comes  for  the  extension  and  contraction  of  the  lines  on  the  opposite  sides  of  the 
axis  about  which  the  moment  is  being  applied.  The  restoring  coefficients  can  be  obtained  by 
assuming  a  small  rotation  (6  or  5)  at  any  point  along  the  center  of  the  tower  /  spar  and  summing  the 
moments  produced.  The  following  is  a  derivation  of  the  roll  restoring  coefficient. 
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where  d  is  the  spar  diameter  and  T  is  the  draft  of  the  spar.  The  center  of  buoyancy  is  at  z  =  -  J72.  The 
center  of  gmvity  for  the  system  at  rest  is  z  =  The  length  of  the  spokes  extending  beyond  the  spar  is 
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Therefore,  the  roll  and  pitch  restoring  coefficients  are  found  to  be 
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Figure  4. 1  Diagram  used  to  calculate  pitch  restoring  coefficient. 

The  restoring  moments  due  to  the  change  in  water  plane  during  pitch  and  roll  are  small  for  this 
system  and  are  therefore  ignored.  It  is  important  to  understand  that  this  restoring  coefficient  is  only 
an  approximation.  The  actual  model  used  in  the  fiilly-coupled  simulation  has  each  of  the  individual 
forces  (lines,  buoyancy,  gravity,  etc.)  that  produce  the  restoring  effect.  The  approximation  of  this 
restoring  coefficient  is  useful  in  a  linear  analysis  of  the  system. 


For  any  applied  moment,  the  amount  of  pitch  or  roll  motion  can  be  found  by  the  following  equation; 
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(4.7) 


Although  wave  forces  will  not  generate  significant  yaw  motion  of  the  floater,  wind  turbine  loading 
can  cause  yaw  motion.  The  restoring  coefficient  for  yaw  is  similar  to  the  surge  restoring  coefficient. 
Movement  of  the  tether  fairleads  produces  horizontal  force  components  which  act  along  the  spokes 
and  buoy  diameter  to  generate  moments  that  counteract  the  yaw  motion.  The  yaw  restoring 
coefficient  and  yaW  displacement  is  found  using  the  following  equations; 


33 


+  L. 


-'66 


-IL 


spoke 


) 


■F 

tethers 


^tethers 


(4.8) 


The  preceding  equations  defined  the  response  of  the  floating  wind  turbine  system  to  external  loads. 
The  restoring  coefficients  are  estimates  of  the  system  behavior.  The  fully-coupled  floating  wind 
turbine  model  does  not  use  these  coefficiente.  Instead,  the  various  physical  properties  (lines,  masses, 
and  buoyancy)  are  modeled  individually  and  their  combined  behavior  determined  by  the  mechanical 
dynamics  solver.  This  will  be  discussed  more  fully  in  Chapter  7.  The  estimates  of  the  restoring 
coefficients  are  useful  in  predicting  the  natural  periods  of  the  system  and  verification  of  the  time 
domain  hydrodynamic  analyses  using  linear  theory. 
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5  Dynamic  Analysis  of  Floater 


The  linear  frequency  domain  method  can  be  used  to  obtain  an  initial  understanding  of  the  importance 
of  the  various  modes  of  motion  of  the  floater.  This  method  allows  the  response  of  a  system  to  be 
determined  analytically  instead  of  empirically.  This  method  allows  for  fast  calculations  that  predict 
motion  in  a  hydrodynamic  environment.  Linear  analysis  is  a  useful  tool  for  determining  the  motion 
of  the  floater  for  small  wave  amplitudes.  However,  the  combined  wind  turbine  /  floater  system  has 
non-linear  excitation  and  quadratic  damping  properties  that  are  better  handled  in  the  time  domain. 
Large  amplitude  waves  result  in  large  horizontal  displacements.  Many  of  the  restoring  coefficients 
for  the  system  are  non-linear  for  large  displaeements.  Also,  the  non-linear  cross-coupling  between 
surge/sway  and  heave  (due  to  set-down)  becomes  more  pronounced  as  the  horizontal  motions 
increase.  The  non-linear  cross-couplings  are  not  aecounted  for  in  a  frequency  domain  analysis. 

The  frequency  domain  analysis  is  primarily  used  to  verify  the  proper  operation  of  hydrodynamic 
computer  code  coupled  with  the  ADAMS/Solver  system  model  and  estimate  natural  periods  for  the 
system.  Frequency  domain  analysis  is  fully  described  in  the  following  references:  Newman  [28], 
Faltinsen  [13],  and  Sclavounos  [31].  Excitation  forces  for  hydrodynamic  motion  arise  due  to  waves 
interacting  with  the  submerged  strueture.  For  infinite  water  depth,  the  following  equations  give  the 
velocity  potential  (0)  and  elevation  (rj)  for  a  wave  of  amplitude  A  and  direetion  p-. 

(I>{x,  y,  z,  0  =  Re  i — I 

I  J  (5.1) 

ri{x,y,t)  =  Re{^  j 

where  k  is  the  wave  number  and  co  is  frequency.  In  deep  water  k  =  o?l%. 

As  stated  earlier  the  long  wave  length  assumption  applies  to  the  linear  hydrodynamic  analysis.  Also, 
the  damping  is  assumed  to  be  linear.  This  assumption  allows  the  motion  equations  for  surge,  heave, 
pitch,  and  yaw  to  be  solved  in  the  frequeney  domain.  Sway  and  roll  are  identieal  to  surge  and  pitch 
and  will  therefore  not  be  discussed  further.  The  wind  turbine  could  excite  the  structure  in  yaw.  It  is 
assumed  that  no  waves  are  generated  by  the  floater,  due  to  the  fact  that  the  wavelengths  of  the 
ambient  waves  are  assumed  to  be  much  greater  that  the  floater  diameter.  The  wave-floater  interaction 
problem  can  therefore  be  solved  in  a  simplified  manner.  The  excitation  forces  are  obtained  from  the 


incoming  wave  potential  and  analytical  expressions  for  the  buoy  and  spokes  added  mass.  For  a  spar 
buoy  the  pitch  and  surge  motions  are  coupled.  The  coupled  surge/pitch  motion  equations  are: 
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and  the  uncoupled  heave  and  yaw  equations  of  motion  are: 

+4,)i; + = F,{i) 

(i^«+A<)A +{S6.)4 +«,  =F,(0 

Mij.  Inertia  coefficient 
Ay:  Added  mass  coefficient 

Bij’.  Damping  coefficients  which  include  damping  due  to  viscous  drag  and  the  wind  turbine. 

Cyi  Restoring  coefficient 

F,:  Excitation  force  or  moment 


In  order  to  solve  the  equations  of  motion  all  of  the  above  coefficients  must  be  solved. 


The  following  analysis  is  taken  from  a  paper  by  Kim  and  Sclavounos  [21],  In  order  for  the  system  to 
be  analyzed  using  linear  frequency  domain  analysis  all  of  the  linear  quantities  must  be  analyzed 
relative  to  the  wave  elevation: 

=  (5.4) 


where  A  is  flie  wave  amplitude.  The  linear  wave  excitation  forces  fj  acting  on  the  platform  are  defined 
relative  to  the  local  coordinate  system  (x,  y,  z)  and  by  virtue  of  linearity  adopt  the  complex 
representation: 

y;  =  Re(/l.X,en  /  =  !,.. .,6  (5.5) 

where  Xi  is  the  complex  amplitude  of  the  exciting  force  or  moment  in  the  indirection.  The  platform 
motions  4;  in  the  six  degrees  of  freedom  also  adopt  the  complex  representation: 

4,  =  Re(A-E,e‘n,  /  =  !,.. .,6  (5.6) 
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where  Hi  is  the  complex  response  amplitude  for  the  platform  in  the  j-direction.  They  are  called  the 
Response  Amplitude  Operators  (RAO’s).  The  classical  forms  of  the  matrix  equations  for  the  RAO’s 
are: 


+  +  (5.7) 

Where  the  damping  has  three  main  contributors:  damping  due  to  the  wind  turbine  {b^),  damping  due 
to  viscous  drag  on  the  vertical  cylinders  (b),  and  the  viscous  drag  on  the  spokes  (bj).  These  damping 
terms  are  combined  into  the  total  damping  term  (5).  Linear  analysis  is  used  to  verify  the  output  of  the 
hydrodynamics  loading  subroutine  in  this  research.  Therefore,  the  wind  turbine  is  assumed  to  be 
secure  during  linear  analysis  of  the  floater.  When  the  wind  turbine  is  secured  all  of  the  damping 
terms  associated  with  it  are  zero.  The  other  assumption  made  during  the  linear  analysis  is  that  the 
lines  provide  negligible  damping.  This  is  a  valid  assumption  since  the  floaters  will  require  relatively 
short  lines  for  the  water  depths  being  considered.  By  making  these  assumptions  the  frequency 
domain  equations  of  motion  are  obtained  when  the  equations  (5.2)  and  (5.7)  are  combined: 
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Knowledge  of  the  platform  resonant  frequencies  for  the  various  floater  modes  of  motion  is  essential 
for  the  design  of  the  platform.  Resonant  vibrations  in  the  floater  can  excite  similar  motion  in  the 
tower  and  wind  turbine  if  natural  frequencies  are  not  taken  into  account.  The  resonant  frequencies  are 
found  from  the  solution  of  the  eigenvalue  problem: 

[-ty^(M^.  +  4)  +  (c^+cJ]H,  =  0,  /  =  !,.. .,6  (5.10) 

The  eigenvalues  of  this  equation  give  the  natural  frequencies  for  the  modes  of  motion  that  are 
applicable  to  the  floater. 


5.1  Coupled  Surge  and  Pitch  Analysis 

The  purpose  of  this  section  is  to  detail  how  the  various  coefficients  necessary  for  the  linear  analysis 
are  derived.  Horizontal  forces  acting  on  the  floater  and  wind  turbine  excite  motion  in  the  surge 
direction.  These  forces  produce  moments  about  the  y-axis  of  the  system.  These  moments  cause 
pitching  motion  of  the  floater.  The  combined  surge  and  pitch  motion  of  the  floater  results  in  the 
movement  of  the  nacelle.  This  motion  can  result  in  changes  in  thrust  and  moment  loading  on  the 
floater  from  the  wind  turbine.  Equations  (5.2)  and  (5.8)  are  used  in  the  analysis  of  surge  and  pitch 
motion. 

5.1.1  Surge  and  Pitch  Added  Mass 

The  inertia  portion  of  the  surge  equation  deals  with  the  calculation  of  the  force  that  is  necessary  to 
accelerate  the  mass  associated  with  the  structure  in  the  surge  direction.  This  inertia  consists  of  mass 
and  added  mass.  The  mass  Mn  is  the  combined  mass  of  the  wind  turbine,  floater,  and  lines.  The 
inertia  term  M55  in  pitch  is  the  mass  moment  of  inertia  of  the  structure  about  the  point  where  rotation 
is  being  assumed. 

The  added  mass  coefficients  for  the  coupled  surge/pitch  equations  are  determined  by  evaluating  the 
forced  oscillations  of  the  floater  in  surge  and  pitch.  Strip  method  is  used  in  this  evaluation.  The 
spokes  at  the  bottom  of  the  structure  must  be  treated  separate  from  the  cylindrical  spar  (see  Figure 
5.1).  Each  strip  along  the  spar  has  the  acceleration: 

^slrip  ~  ^sliip  '  ^5  (5-11) 

where  Zs,rip  is  the  distance  of  the  strip  being  evaluated  from  the  rotation  point  at  the  base  of  the  spar. 
The  added  mass  creates  an  inertia  force  on  the  surface  of  the  floater  due  to  the  accelerations.  Two  of 
the  spokes  contribute  to  the  added  mass  of  the  structure  in  surge  since  they  lie  perpendicular  to  the  x- 
axis.  The  added  mass  of  the  other  two  spokes  has  a  negligible  effect  on  the  structure  in  pitch.  The 
added  mass  of  the  structure  is  determined  by  assuming  that  the  structure  surges  and  pitches  about  the 
point  where  the  added  mass  coefficients  are  being  evaluated  about. 

The  reaction  forces  due  to  forced  oscillations  are  obtained  by  integrating  the  reaction  force  on  each 
strip.  The  added  masses  of  the  cylinder  are  found  using  the  definition  of  added  mass: 

Keac>ion=-Ajj  (5.12) 
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Figure  5.1  Derivation  of  surge  and  pitch  added  mass. 


where  Ostrip  is  the  acceleration  of  the  strip  being  analyzed  due  to  both  translation  and  rotation  at  the 
point  around  which  the  added  mass  is  being  determined.  The  two  dimensional  added  mass  term  for  a 
cylindrical  strip  is  equal  to  the  volume  displaced.  Next  the  contribution  of  the  spokes  is  added  to  the 
surge  added  mass. 
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5.1.2  Surge  and  Pitch  Restoring  Coefficients 


(5.14) 


The  restoring  coefficients  for  pitch  and  surge  are  derived  in  the  previous  chapter.  The  spar  buoy 
which  consists  mainly  of  vertical  cylinders  is  strongly  coupled  in  the  surge  and  pitch  modes  of 
motion.  The  restoring  coefficients  for  a  tethered  spar  cylinder  are; 
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5.2  Heave  Analysis 

Heave  of  the  floater  is  motion  in  the  vertical  direction.  Motion  in  the  heave  mode  is  minimized  due 
to  the  configuration  of  the  lines  bn  the  tension  leg  support  system.  However,  as  the  system  moves  in 
surge  and  sway  modes  it  settles  lower  in  the  water  since  the  legs  remain  nearly  constant  in  length. 
This  results  in  a  coupling  that  is  non-linear.  This  coupling  is  readily  handled  using  a  fully  coupled 
ADAMS  model,  which  will  be  discussed  later.  Linear  analysis  can  be  used  to  evaluate  small 
amplitude  heave  oscillations  and  determination  of  the  heave  natural  frequency. 


5.2.1  Heave  Inertia  and  Added  Mass  Coefficients 


The  inertia  term  for  heave  motion  is  the  total  mass  of  the  combined  structure.  Added  mass  in  heave  is 
due  to  the  mass  of  water  under  the  spar  that  is  accelerated  when  the  structure  moves  in  the  vertical 
direction.  This  added  mass  can  be  approximated  by  the  displacement  of  a  half  sphere  of  water  with  a 
radius  equal  to  that  of  the  base  of  the  spar  cylinder  [31],  [15].  The  spokes  also  contribute  to  the  heave 
added  mass.  Heave  added  mass  is  given  by: 
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where  r  is  the  radius  of  the  bottom  of  the  spar 

5.2.2  Heave  Restoring  Coefficient 

The  heave  restoring  force  is  generated  by  two  mechanisms.  First,  as  the  buoy  heaves,  the  line  tension 
changes  and  very  large  forces  are  generated  due  to  the  relatively  high  modulus  of  elasticity  for  the 
tethers.  The  other  part  of  the  restoring  force  is  the  change  in  buoyancy  of  the  floater  as  it  moves  in 
the  vertical  direction.  This  second  factor  is  roughly  two  orders  of  magnitude  smaller  than  the  line 
restoring  force  for  water  depths  in  the  range  of  100  to  300  meters  for  relatively  stiff  lines.  The  heave 
restoring  coefficient  is  given  by  equation  (5.17). 
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where  Ay^  is  the  area  of  the  spar  buoy  waterplane. 

5.3  Yaw  Motion 

Yaw  motion  in  the  support  system  is  exited  by  yaw  moments  generated  on  the  wind  turbine  rotor  and 
wave  loading  of  the  submerged  spokes.  These  excitations  are  counteracted  by  the  restoring  force  of 
the  tethers  which  act  over  the  moment  ann  of  the  spokes.  As  the  floater  rotates  the  fairlead  position  is 
moved  away  from  the  rest  position  and  a  horizontal  force  component  is  produced.  The  force  results 
in  a  counter  moment  which  opposes  the  yaw  displacement. 


5.3.1  Yaw  Inertia  and  Added  Mass  Coefficients 


The  yaw  moment  of  inertia  and  added  mass  are  found  by  combining  the  contributions  of  the 
individual  components.  Cylindrical  components  centered  on  the  yaw  axis  do  not  contribute  to  the 
yaw  added  mass.  Therefore  only  the  spokes  make  a  contribution  to  the  yaw  added  mass.  The  spokes 
are  modeled  as  rigid  box  structures  that  have  a  square  cross  section. 

The  yaw  moment  of  inertia  is  the  moment  of  inertia  about  the  z-axis.  The  yaw  added  mass  of  the 
spokes  is  found  using  the  strip  method  along  the  spokes.  The  total  yaw  added  mass  for  all  four 
spokes  with  length  and  width  (w)  is  given  by  [28]: 
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5.3.2  Yaw  Restoring  Coefficient 

For  small  yaw  motions,  a  linear  restoring  coefficient  can  be  assumed.  This  coefficient  is  similar  to 
the  linear  heave  restoring  coefficient,  since  both  coefficients  are  generated  by  movement  of  the 
fairlead  position.  As  the  system  yaws  the  tension  in  all  of  the  tethers  produces  a  moment  in  the 
opposite  direction  of  the  displacement.  The  linear  yaw  restoring  coefficient  is  given  by: 
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6  Hydrodynamic  Loading 


Waves  act  on  the  submerged  portion  of  the  floater  and  cause  motion  of  the  system.  This  chapter  deals 
with  the  determination  of  wave  excitation  on  the  submerged  portion  of  the  tower,  the  spar  buoy,  and 
the  spokes.  The  forces  and  moments  due  to  the  waves  are  found  using  potential  theory.  The  velocity 
potential  of  a  planar  progressive  wave  with  frequency  (co),  wavenumber  (A),  and  direction  (P)  is  given 
by  the  following  equation: 

=  (6.1) 

where  k  is  the  wavenumber 

6.1  Interface  of  Hydrodynamic  Loading  Module  with  ADAMS 

Prior  to  beginning  the  discussion  of  the  derivation  of  hydrodynamic  loading  it  is  important  to 
understand  how  the  calculations  are  used  in  the  fully  coupled  dynamic  analysis.  Every  time  step  the 
ADAMS  solver  sends  the  hydrodynamic  loading  subroutine  the  position,  velocity,  and  acceleration  of 
a  point  on  the  part  for  which  loads  are  being  determined.  The  rotational  components  are  relative  to 
that  parts  body  based  coordinates  (roll,  pitch,  and  yaw),  not  the  rotations  about  the  global  coordinate 
system.  The  local  rotations  change  the  loading  on  the  system.  Loads  are  then  calculated  relative  to 
the  body  based  Coordinate  system  and  sent  back  to  the  solver.  Every  part  in  an  ADAMS  model  has 
its  own  body  based  coordinate  system,  which  moves  with  the  part.  Roll,  pitch,  and  yaw  are  the 
angular  differences  between  the  body  coordinate  system  axes  and  the  global  coordinate  axes.  The 
solver  combines  the  loads  on  the  six  different  floater  components  and  combines  them  with  the 
aerodynamic  loads  from  the  45  blades  elements  in  the  rotor  while  performing  the  totol  system 
dynamic  analysis. 

6.2  Gl  Taylor’s  Loading  Approximation 

The  diameter  of  the  floater  is  small  compared  to  the  wavelength  of  most  of  the  waves  of  interest.  The 
diameter  of  the  spar  will  fall  into  a  range  of  6  m  to  30  m.  A  wave  with  period  (7)  of  10  seconds  will 
have  a  wavelength  of  approximately  (wave  period  rule  of  thumb):  X  ^  +  TV  2  =  150  m.  Therefore 

the  wavelength  is  much  greater  than  the  diameter  of  the  spar.  The  long  wavelength  approximation 
can  therefore  be  used  when  determining  the  excitation  forces  and  mornents  on  the  structure  by  the 


incoming  waves.  This  approximation  allows  calculations  to  be  made  at  the  center  of  the  body  and 
assume  that  the  wave  potential  is  constant  across  the  body.  G1  Taylor’s  equation  for  wave  excitation 
is  given  by  equation  (6.2). 


=  (pV+4) 


du.  du. 
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V  dt  di  j 


dt^ 


(6.2) 


Where  An  is  the  added  mass  of  the  structure  in  the  direction  of  the  force  being  calculated  and  x,  is  the 
displacement  of  the  structure  in  the  direction  of  interest.  The  velocity  of  the  undisturbed  wave  field 
in  the  direction  of  the  force  at  the  center  of  the  structure  is  denoted  by  m,.  The  above  equation  is 
called  G1  Taylor  formula.  It  is  applicable  to  structures  whose  principal  dimension  is  much  less  than 
the  wavelength  of  the  incident  waves. 

6.2.1  Excitation  Loads  on  Submerged  Cylinders 

The  G1  Taylor  formula  is  used  in  a  strip-wise  integration  to  evaluate  the  total  hydrodynamic  force  on 
the  wind  turbine  floater.  The  motion  of  the  floater  is  accounted  for  in  these  integrations.  The 
horizontal  and  vertical  forces  on  the  floater  are  phase  shifted  by  90°.  The  following  is  a  derivation  of 
the  forces  exerted  on  a  vertical  cylinder.  For  the  following  discussion,  refer  to  Figure  5.1.  The  strip 
method  for  excitation  forces  are 
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,  A :  wave  amplitude 


where  u. 


0  to  the  leading  order.  Which  results  in  the  above  excitation  equation  becoming: 


dt 


ICO 


- ik  cos  6 

dx 


(6.3) 


The  velocity  components  of  the  undisturbed  water  at  position  are  found  by  taking  the  x,  y,  and  z 
derivatives  of  the  wave  potential  equation. 
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The  added  mass  terms  in  surge  /  sway  for  a  strip  of  thickness  dz  and  the  heave  added  mass  due  to  the 
bottom  of  the  spar  cylinder  are  given  by  the  following  equations: 


a, ,  =  aj2  =  Area^,  ■  p^-  dz  =  zrj  •  dz 

A  2*;r  3 
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where  Ass  is  the  added  m^s  of  a  half  sphere  of  water  with  a  radius  equal  to  that  of  the  spar  cylinder. 
The  final  part  of  the  equation  that  is  need  is  the  acceleration  terms: 
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Now,  by  using  these  terms  in  (6.2),  the  forces  can  be  found.  The  strip  method  is  used  to  sum  up  the 
contributing  horizontal  forces  along  the  length  of  the  vertical  cylinder.  For  heave  the  excitation  force 
is  assumed  to  the  vertical  force  acting  on  a  displaced  volume  of  water  at  the  depth  of  the  floater 
bottom.  This  volume  is  assumed  to  be  that  of  a  half  sphere  with  a  radius  equal  to  the  bottom  of  the 
cylinder.  The  forces  on  the  individual  strips  are  given  by  the  following  equations. 
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where  Area  is  the  cross-sectional  area  of  the  strip  being  analyzed  (Area  =  7fr^y,i„dJ).  The  first  two 
equations  in  (6.4)  apply  for  both  the  spar  and  the  submerged  portion  of  the  tower.  The  forces  on  the 
bottom  of  the  spar  in  the  z-direction  are  treated  separately  since  only  the  bottom  of  the  structure  has 
vertical  wave  excitation  forces  acting  on  it.  The  derivation  of  the  pitch  and  roll  moments  are  similar 
to  those  of  the  forces  in  the  x  and  y  directions.  Forces  acting  on  the  cylindrical  strips  are  multiplied 
by  the  moment  arm  to  the  bottom  of  the  cylinder  (the  point  at  which  forces  and  moments  are 
summed).  The  line  configuration  of  the  tension  leg  platform  results  in  very  small  motion  in  the  pitch 
and  roll  directions  so  long  as  the  all  the  lines  remain  in  tension  and  stiff  lines  are  used.  Since  the  spar 
cylinder  may  not  extend  all  the  way  to  the  surface,  the  top  of  the  cylinder  is  assumed  to  be  at  depth  — 
di  and  the  bottom  is  at  -dj. 

dM^=2- pArea  ■  Aco^  cos p ■  | +  2) •  e'^dz -  A^s_,,„- 
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6.2.2  Spoke  Excitation  Forces 

The  spokes  of  the  floater  are  used  to  give  the  lines  longer  moment  arms  to  act  over.  Longer  spokes 
result  in  lower  changes  in  stress  when  pitch  and  roll  moments  are  applied  to  the  system.  The  actual 
shape  of  the  spokes  is  not  known.  They  could  be  open  trusses  or  closed  box  structures  with 
buoyancy.  For  modeling  purposes  the  spokes  are  treated  as  boxes  with  a  characteristic  dimension 
equal  to  the  length  of  one  side  (w).  The  wave  excitation  on  the  spokes  can  be  changed  by  varying  the 
characteristic  dimension.  It  should  be  noted  that  the  forces  and  moments  on  the  spokes  are  low 
compared  to  that  of  the  other  portions  of  the  floater.  The  hydrodynamic  forces  on  the  spokes  are 
calculated  at  their  center  of  mass.  Assuming  the  spokes  are  attached  to  the  bottom  of  the  spar  (z  =  - 
^2),  the  loading  equations  for  the  spokes  are: 

Fy  =  +  ^22  )A(o^  sin  p  ■  me  j  _ 

F.  =  ipVoK^,e  +  A,,)Aco^  ■  j  _ 

where  x^,x^,  and  x^  are  the  x,  y,  and  z  accelerations  of  the  center  of  mass  of  the  spoke.  The 

contribution  of  spokes  to  the  floater  loading  is  small  compared  to  the  vertical  cylinders  that  make  up 
the  submerged  portion  of  the  tower  and  the  spar.  This  is  due  to  the  fact  that  the  spokes  have  a 
relatively  small  volume  and  are  mounted  near  the  bottom  of  the  floater.  These  two  characteristics 


45 


reduce  their  influence  on  excitation  forces.  Moments  due  to  spoke  excitation  forces  are  generated  by 
applying  ftese  forces  at  the  center  of  the  spokes. 

6.3  Viscous  Drag 

Most  of  the  underwater  structure  of  the  floater  /  wind  turbine  system  consists  of  a  buoyant  cylinder 
and  the  extension  of  the  tower  into  the  water.  Viscous  drag  on  a  vertical  cylinder  can  be  calculated 
using  Morison’s  equation.  The  horizontal  differential  force  in  the  x-direction  (dF)  on  a  strip  of  length 
dz  of  a  cylinder  using  Morison’s  equation  is: 

dF  =  p^^^dzC^a^  +^C^Ddz\u-x^\{u-jc^)  (6.7) 

Where  the  first  term  in  this  equation  is  the  inertial  forces  and  has  been  handled  using  the  GI  Taylor 
equation.  The  second  term  is  the  viscous  damping  term.  It  is  used  to  calculate  the  viscous  drag  on 
the  submerged  portions  of  the  floater.  There  are  three  non-dimensional  factors  that  influence  the 
viscous  flow  around  a  submerged  body: 

Keulegan-Carpenter  number  KC  = 

Reynolds  number  Rn=I/Z>/K  (6.8) 

Rn/KC  P  =  — 

vT 

where  a  is  the  amplitude  of  oscillation,  D  is  the  characteristic  length  of  the  body,  T  is  the  period  of 
oscillation,  and  v  is  the  kinematic  viscosity  (15°  C  salt  water:  1.19E-6  mVs).  As  described  in  the 
notes  by  Molin  [26],  the  determination  of  viscous  effects  on  submerged  structures  is  very 
problematic.  For  this  reason  simplifying  assumptions  are  made  to  handle  the  drag  forces  on  the 
cylinders  and  spokes.  Drag  coefficients  are  assumed  to  be  constant  values.  Testing  has  shown  that 
drag  coefficients  for  oscillating  bodies  change  with  amplitude  and  frequency  along  with  floater 
geometry. 

6.3.1  Drag  on  Submerged  Cylinders 

For  the  floating  wind  turbine  system,  the  spar  cylinders  used  in  the  floaters  have  diameters  of  fall  in 
the  range  of  7  to  12  meters.  The  amplitude  of  motion  in  most  sea  states  is  in  the  range  of  1  to  12 
meters  for  surge  and  sway.  For  cylinders  with  this  range  of  diameters  and  amplitudes  of  motion  the 


range  for  KC  is  0.5  to  10.  According  to  Faltinsen  (p.  239,[13]),  for  KC<10  the  drag  coefficient  (Co) 
for  a  cylinder  can  be  estimated  using  the  KC:  =  0.2KC .  A  large  range  of  drag  coefficients  are 

therefore  present  on  the  structure  (0.1  to  2.0). 

A  constant  drag  coefficient  of  0.7  is  used  in  equation  (6.7)  for  determining  drag  on  the  cylinders.  The 
relative  velocity  is  the  difference  between  the  horizontal  velocity  (u)  of  the  water  at  the  depth  z  and 
the  velocity  of  the  structure  in  the  surge  direction  (x, ).  By  substituting  the  fluid  velocity  terms  from 
potential  flow  theory  the  second  term  in  Morison’s  equation  in  monochromatic  waves  becomes: 

dFjrag  =  Y  C^Ddz  ^Acocos  P  •  91e  |  |  cos  |  j  (59) 

where  D  is  the  diameter.  This  equation  cannot  be  directly  integrated,  therefore  the  total  drag  force  on 
the  cylinder  at  each  time  step  is  found  by  dividing  the  submerged  cylinders  into  strips  and  summing 
the  forces  for  these  strips.  The  contribution  this  force  has  on  the  pitch  and  roll  moments  is  found  by 
multiplying  each  dFd„g  by  the  moment  arm  of  that  strip  (distance  from  the  strip  to  the  point  about 
which  the  moments  are  being  calculated).  The  differential  moment  (about  the  base  of  the  cylinder,  z 
-  d:)  produced  by  the  drag  force  on  a  cylindrical  strip  is  given  by  the  following  equation: 


(6.10) 


6.3.2  Drag  on  Spokes 

Morison’s  equation  is  also  used  to  determine  the  viscous  drag  on  the  four  spokes  of  the  structure.  The 
drag  force  in  the  surge,  sway,  and  heave  modes  is  found  by  using  the  drag  equation  for  a  single  strip 
with  the  area  term  being  the  frontal  area  the  spoke  presents  to  the  wave.  The  spokes  have  a  square 
cross-section  so  a  drag  coefficient  of  1 .0  is  used.  The  drag  forces  are  given  by  the  equation: 


=  YCo(w  Z,cos(4  - >0))|^rycosy5-  - jc, 
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where  w-Zcos(^g  - P)  is  the  effective  frontal  area  of  the  spoke  to  the  wave,  w  is  the  spoke  width, 
and  L  is  the  spoke  length.  The  x,  y,  and  z  positions  and  velocity  of  the  spoke  correspond  to  the  point 


at  the  center  of  mass  for  the  spoke.  ^  is  the  yaw  angle  of  the  spoke.  The  equation  is  evaluated  at  a 
depth  z  equal  to  the  depth  of  the  center  of  the  spoke.  Damping  in  the  yaw  mode  of  motion  is  due  to 

viscous  dmg  on  the  spokes. 

6.4  Representation  of  Random  Waves 

Hydrodynamic  analysis  of  a  floating  structure  in  a  realistic  sea  environment  requires  that  statistical 
estimates  be  made  to  represent  different  sea  states.  Irregular  seas  are  generated  by  linear 
superposition  of  waves  of  different  frequencies.  Since  the  individual  waves  are  independent  of  each 
other  they  can  be  treated  separately  when  calculating  the  excitation  forces  on  the  floater.  In  depth 
discussions  on  this  procedure  can  be  found  in  Sclavounos,  2003[31]  and  Faltinsen,  1990  [13]. 

A  sea  state  can  be  represented  by  the  ambient  wave  spectral  density  S(m).  llie  spectrum  is  a 
statistical  representation  of  an  ocean  environment.  As  the  sea  state  of  the  ocean  environment 
increases  its  spectral  density  increases.  Figure  6.1  is  a  graph  of  the  spectral  density  versus  frequency. 
The  spectrum  is  broken  up  into  frequency  bands  of  width  Aco. 
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Figure  6.1  Ambient  Wave  Spectral  Density 

The  spectral  density  is  used  to  determine  the  amplitude  for  a  frequency  band  of  width  Am  by  using  the 
equation  (6.12).  This  equation  allows  a  group  of  waves  of  amplitudes  Aj  to  be  generated  for  use  in 
the  excitation  analysis  of  the  floater.  As  the  spectral  density  of  a  band  increases  so  does  its 
amplitude. 
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(6.12) 


(6.13) 


Where  H,/}  is  the  significant  wave  height  which  is  defined  as  the  mean  of  the  highest  1/3  of  waves. 
Tj  is  the  mean  wave  period  for  the  wave  spectrum  being  represented  and  is  calculated  using  (6.14). 


7;  =2;r^ 

m, 


(6.14) 


The  moments  are  defined  as: 


=  j[<y*  •5(<y)]c?ry 


(6.15) 


The  standard  deviation  {&)  of  the  spectrum  is  calculated  by  setting  k  equal  to  zero. 


=  js(^a))d{V  =  a^ 


(6.16) 


The  standard  deviation  can  be  used  to  calculate  the  significant  wave  height  (Sclavounos,  2003): 


//  w  s  =  4(7 


(6.17) 


Another  spectrum  that  was  used  during  this  research  was  the  Joint  North  Sea  Wave  Project 
(JONSWAP)  type  spectrum  which  is  given  in  Faltinsen,  1 990  [  1 3]. 
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where 
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and 
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The  frequencies  of  the  individual  wave  components  determine  if  the  combined  wave  output  will 
appear  to  be  random  or  periodic.  If  the  frequency  bands  are  equally  spaced  Mong  the  frequency  axis 
and  the  center  of  the  bands  are  used  for  the  associated  wave  frequency,  the  output  will  be  periodic. 
Therefore,  it  is  important  that  the  bands  frequencies  not  be  multiples  of  each  other.  Another 
important  consideration  is  computational  time.  The  processing  time  for  a  simulation  is  directly 
proportional  to  the  number  of  bands  being  used  in  the  analysis.  It  is  desirable  to  only  use  frequency 
components  which  have  an  spectral  energy  associated  with  them  to  describe  a  given  sea  state. 

The  technique  used  in  this  thesis  involved  a  two  step  process.  First,  the  wave  spectrum  was  generated 
using  1000  bands  as  described  previously.  Then  the  desired  number  of  bands  was  generated  by 
choosing  the  band  widths  such  that  the  area  under  the  wave  spectrum  curve  was  approximately  equal 
for  each  of  the  frequency  ranges.  The  area  of  the  new  frequency  bands  was  found  using  the  following 
equation: 

KKK) 

-  (619) 

The  spectrum  was  broken  up  into  bands  which  had  equal  areas  (see .  This  resulted  in  the  bands  being 
more  closely  spaced  near  the  peak  and  wider  spaced  at  the  tails  of  the  spectrum.  If  the  bands  were  of 
equal  width,  the  tails  (high  frequency  and  low  frequency)  would  consist  of  wave  components  with 
very  low  (sometimes  zero)  amplitudes.  Processor  time  would  be  wasted  analyzing  these  extraneous 
bands.  Instead,  by  using  bands  with  equal  spectral  energy,  every  wave  component  contributed 
significantly  to  the  floater  loading  and  processor  usage  was  optimized.  The  frequency  of  each  band 
was  found  by  taking  the  moment  of  area  for  the  band  and  dividing  by  the  area  instead  of  the  center  of 
the  band. 
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Figure  6.2  Example  of  sea  state  6  wave  spectrum  divided  into  bands  with  equal  spectral  energy. 

The  two  wave  spectra  given  by  (6.13)  and  (6.18)  are  used  to  generate  random  sea  states  for  use  in  the 
hydrodynamics  analysis  of  the  floater  portion  of  the  floating  wind  turbine  system.  The  amplitudes 
associated  with  each  frequency  band  is  then  calculated  using  a' (6. 12)  with  the  area  contained  by  the 
band  substituted  for  S((Oi)Aco  term.  The  superposition  of  the  different  waves  results  in  an  irregular  sea 
elevation.  The  wave  elevation  (q)  for  long  crested  seas  is  given  by: 

j 

=  +  +  (6.20) 

N^\ 

where  ej  is  the  random  phase  angle  for  the  frequency  band  drawn  from  a  uniform  distribution. 

The  time  varying  wave  elevation  is  very  important  in  the  floater  load  analyses.  Nonlinear  drift  forces 
result  due  to  the  fact  that  more  of  the  structure  is  in  the  wave  fteld  when  a  crest  is  passing  than  when 
.  the  trough  is  passing.  These  forces  will  be  discussed  in  the  following  section. 

6.5  Nonlinear  Hydrodynamic  Forces 

The  time  varying  wave  elevation  generates  second  order  drift  forces.  These  forces  cause  the  floater  to 
be  displaced  from  the  rest  position  and  oscillate  about  a  mean  drift  position.  Second  order  forces  are 
not  reflected  in  linear  analysis.  This  is  due  to  the  fact  that  linear  theory  assumes  that  all  analysis  can 
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be  carried  out  using  a  constant  z  equal  to  zero  for  the  upper  portions  of  the  floater.  The  potential  flow 
equations  give  identical  magnitudes  for  forces  acting  in  both  the  direction  of  wave  propagation  and 
the  reverse  direction.  This  results  in  linear  theory  producing  oscillating  forcing  function  with  a  mean 
of  zero.  The  mean  horizontal  drift  force  due  to  wave  diffraction  and  radiation  effects  is  assumed  to  be 
negligible  because  of  the  small  cylinder  diameter  relative  to  the  ambient  wavelength. 

In  this  research,  nonlinear  hydrodynamic  forces  are  modeled  using  the  approximation  method 
described  in  Faltinsen,  1990  [13].  This  method  involves  accounting  for  changes  in  wave  elevation  by 
performing  strip  theory  calculations  over  the  actual  volumes  instead  of  using  the  mean  free  surface  as 
one  endpoint  in  the  computations.  Wave  elevations  above  the  mean  free  surface  are  a  problem  due  to 
the  e**  term  in  potential  equations,  which  would  result  in  very  large  fluid  velocities  and  accelerations 
in  steep  waves.  The  approximation  method  involves  setting  z  equal  to  zero  in  the  ambient  wave 
kinematics  equations  above  the  mean  free  surface  and  using  the  wave  elevation  for  z  in  troughs. 
Figure  6.3  shows  the  distributions  of  acceleration  and  velocity  under  the  waves  as  they  change  with 
depth. 

This  method  of  modeling  nonlinear  forces  is  useful  in  a  time  domain  analysis  of  the  floater.  The 
hydrodynamic  loads  are  determined  by  summing  the  inertial  and  viscous  loads  for  the  strips.  These 
forces  are  calculated  using  the  equation  of  GI  Taylor  and  Morison. 


Figure  6.3  Distribution  of  horizontal  velocities  and  accelerations. 


# 


52 


7  Computer  Model 


The  wind  turbine  /  floater  system  was  best  modeled  in  the  time  domain.  This  was  due  to  the  fact  that 
accurate  representations  of  the  various  components  and  loading  conditions  involve  many  non-linear 
terms.  Several  different  wind  turbine  computer  models  are  used  in  industry  to  simulate  various 
loading  conditions  and  predict  machine  response.  However,  none  of  these  codes  simulate  a  wind 
turbine  mounted  on  a  floating  base.  Modifications  to  existing  codes  and  new  subroutines  were 
written  to  model  the  fully  coupled  response  of  the  combined  floater  wind  turbine  system. 

Several  different  computer  codes  were  used  in  this  research.  The  majority  of  the  work  involved  a 
hydrodynamics  program  that  modeled  the  excitation  forces  on  the  floater,  a  modeling  subroutine  that 
generated  a  computer  model  of  the  floater,  ADAMS  and  FAST.  FAST  was  used  to  generate  the 
ADAMS  model  of  the  wind  turbine  and  to  perform  additional  wind  turbine  calculations.  The 
ADAMS  input  file  was  modified  to  include  the  floating  base.  ADAMS  was  used  to  perform  the  time 
domain  simulation  of  the  system  model.  The  following  chapter  discusses  the  computer  codes  used  in 
the  analysis.  Also,  a  description  of  the  analysis  model  is  given. 

7.1  FAST 

The  FAST  (Fatigue,  Aerodynamics,  Structures,  and  Turbulence)  code  is  a  wind  turbine  simulation 
code  produced  by  the  National  Renewable  Energy  Laboratory  (NREL).  It  is  a  combination  of  three 
different  codes:  FAST2  (2-bladed  horizontal  axis  wind  turbines  (HAWTs),  FAST3  (3-bladed 
HAWTs),  and  AeroDyn  (aerodynamics  subroutines  for  wind  turbines).  A  full  description  of  FAST 
can  be  found  in  the  FAST  User’s  Guide  [4]. 

The  Fast2Adams  subroutine  in  FAST  was  used  to  generate  the  portion  of  the  ADAMS  model  that 
represents  the  wind  turbine.  This  subroutine  was  modified  so  that  it  calls  the  fltr _prop  subroutine, 
which  adds  the  floater  representation  to  the  ADAMS  model  file.  The  Fast2Adams  subroutine  also 
generates  the  ADAMS  control  file  which  governs  the  operation  of  the  solver  analysis. 

7.2  ADAMS 

The  simulation  of  the  total  system  is  carried  out  in  the  time  domain  using  MSC.ADAMS.  ADAMS 
(Automatic  Dynamic  Analysis  of  Mechanical  Systems)  is  a  commercially  available  software  package 
that  is  used  for  dynamic  analysis.  It  is  produced  by  Mechanical  Dynamics,  Inc.  ADAMS  is  a 


program  that  consists  of  many  subprograms  that  can  be  used  to  accomplish  various  engineering 
analysis  tasks.  ADAMS/Solver  was  used  to  combine  the  loading  on  the  floater  with  the  loading  on 
the  wind  turbine  to  create  a  dynamically  coupled  system  . 

ADAMS/Solver™  is  a  numerical  analysis  application  that  automatically  solves  the  equations  of 
motion  for  kinematic,  static,  quasi-static,  and  dynamic  simulations.  It  is  useful  for  analyzing  systems 
that  have  complex  loading  and  flexural  properties.  The  combined  floater  and  wind  turbine  system 
requires  this  kind  of  analysis  tool  in  order  to  achieve  the  proper  coupling.  The  following  section 
describes  how  the  ADAMS/Solver  works  and  how  it  is  used  in  this  research.  The  ADAMS/Solver 
subprogram  will  henceforth  be  referred  to  as  ADAMS.  The  description  of  ADAMS  operation  is 
taken  from  the  ADAMS  Users  Manual  [1]. 

7.2.1  Description  of  an  ADAMS  Model 

The  inputs  to  ADAMS  are  a  model  file  and  a  control  file.  The  model  file  describes  the  system  being 
analyzed.  Mechanical  systems  are  represented  with  interconnected  parts.  Flexible  components  are 
modeled  by  splitting  the  component  up  into  smaller  parts  which  are  then  connected  to  adjacent  parts 
with  springs  and  dampers  that  represent  the  flexural  and  damping  properties  of  the  component  being 
modeled.  For  example  each  blade  consists  of  15  parts  which  range  from  the  blade  root  out  to  the 
blade  tip.  A  system  like  the  wind  turbine  and  floater  combination  consists  of  approximately  60  parts. 

Each  part  is  modeled  as  a  lump  mass  with  a  position  and  inertial  properties  which  reflect  the  region 
being  modeled.  Field  statements  are  used  to  connect  a  part  to  adjacent  parts.  These  statements 
describe  the  translational  and  rotational  restoring  and  damping  properties  of  the  component  in  the 
region  of  the  part.  The  springs  and  dampers  represented  by  the  field  statements  are  attached  to  points 
on  the  parts  called  markets.  Figure  7.1  is  an  example  of  how  individual  parte  are  interconnected  with 
these  field  statements.  When  the  number  of  parts  used  to  describe  a  component  is  chosen  correctly  an 
accurate  prediction  can  be  made  of  the  dynamic  response  for  the  component. 
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Figure  7.1  Example  of  ADAMS  model. 


Once  the  parts  physical  properties  and  interconnections  are  described  the  next  step  in  modeling  is 
application  of  external  forces.  Parts  can  also  be  connected  using  various  kinds  of  joints.  Forces  are 
applied  to  markers  on  parts.  The  force  of  gravity  is  applied  to  the  center  of  gravity  marker.  External 
forces  can  be  applied  to  any  marker  on  a  body. 

Forces  are  calculated  using  functions  or  external  subroutines.  The  forces  exerted  on  the  wind  turbine 
blade  elements  are  computed  by  the  AeroDyn  subroutine.  The  hydrodynamic  forces  on  the 
submerged  portions  of  the  systems  are  computed  by  the  fltrjoads  subroutine.  These  two  subroutines 
are  described  further  in  other  sections. 

7.2.2  ADAMS  Control  File 

The  operation  of  the  ADAMS  solver  is  governed  by  the  control  file.  The  commands  in  this  file 
initiate  motion,  unlock  or  lock  joints,  and  dictate  the  type  of  numerical  integration  used  during  the 
simulation.  The  combined  wind  turbine  /  floaters  system  is  relatively  stiff;  this  fact  leads  to  the 
choice  of  a  stiff  numerical  integration  technique  for  the  simulation. 

The  dynamic  analysis  of  a  mechanical  system  consists  essentially  of  numerically  integrating  the  nonlinear 
differential  equations  of  motion.  The  type  of  system  (soft  or  stiff)  governs  the  choice  of  numerical 
integration  used  by  the  ADAMS/Solver.  For  this  research  the  best  integrator  was  found  to  be  the  GSTIFF 
routine.  This  routine  is  good  for  high  speed  and  accurate  analysis  of  a  wide  range  of  systems.  A  thorough 
description  can  be  found  in  the  ADAMS/Solver  Users  Guide. 

7.3  Floater  Model  Description 

The  input  model  for  the  ADAMS  analysis  is  generated  to  model  the  wind  turbine  and  floater  as  a 
combined  system.  The  wind  turbine  portion  of  the  model  is  generated  by  the  FAST  pre-processor. 
The  floater  portion  is  then  added  to  the  model  file.  The  floater  consists  of  several  different  parts: 
submerged  tower,  spar  cylinder,  spokes,  and  tethers.  All  of  these  parts,  their  associated  properties, 
and  loading  conditions  had  to  be  determined  by  new  subroutines  that  interfaced  with  ADAMS. 

This  program  uses  the  data  in  the  floater  characteristic  input  file  and  generates  the  ADAMS  input 
file.  Figure  7.2  is  an  example  of  a  floater  input  file.  The  jltr _prop  subroutine  serves  two  purposes. 
This  subroutine  calculates  the  physical  parameters  for  the  floater,  and  then  writes  the  description 
statements  to  the  model  file.  The  second  purpose  is  to  perform  additional  calculations  which  are  used 
to  verify  the  output  of  the  time  domain  analysis  against  linear  analysis.  Figure  7.3  is  an  example  of 
the  parameters  calculated  in  the  fltr _prop  subroutine. 


10.0 

!  Floater  Diameter  ALL  UNITS  SI 

10.0 

!  Tower  Diameter 

15.0 

•  Depth  of  tower 

0.0 

!  Length  of  cylinder  (Height) 

0.0 

!  Height  of  spokes  on  cylinder  (0  =  bottom) 

12.0 

!  Length  of  spokes 

2.0 

!  Characteristic  width  of  spokes  (square  cross-section) 

0.0253 

!  Thickness  of  spar 

0.104 

!  Line  Diameter  (m) 

1 

!  Number  of  lines 

743.0 

!  Breaking  Load  in  tonnes 

0.045 

!  Maximum  line  Extension 

200.0 

!  Water  depth 

0.0 

!  Wave  Direction  (degrees) 

0.0 

f  Wave  Direction  Range  (+/-  this  value) 

30 

!  Number  of  Bands  (integer  value) 

.false. 

!  Buoyant  spokes 

2 

!  (1)  PierceMoskowItz  or  (2)  JONSWAP 

5.0 

!  Significant  wave  height  (H1/3) 

8.686 

!  Mean  wave  period,  T1 

Figure  7.2  Example  of  floater  input  file 


Calculated  Floater  Properties 


46675.9 

Spar  Mass 

0.0 

Mass  Cylinder 

15558.6 

Mass  Spar  Plate 

93351.9 

Submerged  Tower  Mass 

79239.6 

Mass  of  Spokes 

219267.4 

Total  Floater  Mass 

201059.2 

Mass  of  Wind  T urbine 

420326.6 

System  Mass 

1212261,1 

Buoyant  Mass 

791934,5 

Reserve  Buoyancy 

188,41 

Percent  Reserve  Buoyancy 

0.0340 

Total  Area  of  Lines 

885.5 

Spring  Constant 

23306.29 

Intial  Line  Stress 

2.19 

Intial  Line  Stretch 

52.5060 

CG  of  Wind  Turbine 

-11.8069 

CG  of  Floater 

18.9565 

CG  of  System 

1225191.9 

A11andA22 

574157.8 

A33 

106552323.0  A44  and  A65 

31229673.8  A66 

721810181.6  Ixx 

721810181.6  lyy 

60009264.4  Izz 

-41993.9 

C11andC22 

-4334783.3 

C33 

-460989441 .0  C44  and  C55 

-713896.9  C66 

ivaiurai  rrequencies  ot  rioaier  oysiem 

0.1598 

Surge  and  Sway  Natural  Frequencies  (rad/sec) 

2.0878 

Heave  Natural  Frequen(^(rad/sec) 

0.7460 

Roll  and  Pitch  Natural  Frequencies  (rad/sec) 

0.0885 

Yaw  Natural  Frequency  (rad/sec) 

Figure  7.3  Output  of fltr _prop  subroutine. 
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The  floater  model  describes  the  physical  characteristics  of  the  floater  so  that  it  realistically  represents 
the  floater  in  the  time  domain  analysis.  The  floater  is  modeled4  as  a  single  rigid  body  which  is 
attached  to  the  tower  base  with  a  spring  and  damper.  Figure  7.4  shows  the  locations  for  the  loading 
applied  to  the  floater.  The  following  loading  characteristics  are  incorporated  into  the  ADAMS  model; 

•  Force  due  to  gravity. 

•  Buoyancy  force  for  the  spar  portion. 

•  Buoyancy  force  for  the  submerged  tower. 

•  Buoyancy  forces  for  the  spokes  (if  buoyant). 

•  Forces  due  to  tethers. 

•  Hydrodynamic  forces  (calculated  in  fltrjoads  subroutine) 


Hydrodynanic 

Figure  7.4  Loading  for  floater  model. 

The  hydrodynamic  forces  and  moments  are  applied  at  different  points  in  the  structure.  As  previously 
discussed  “markers”  are  used  in  the  model  as  application  points  for  forces  and  moments.  The  loading 
markers  for  the  submerged  portion  of  the  tower  and  the  spar  are  located  at  their  bases.  The 
application  points  for  hydrodynamic  forces  on  the  spokes  are  located  at  their  geometric  centers.  The 
solver  can  then  use  these  loads  and  their  application  points  to  determine  the  effect  on  the  rest  of  the 
system.  The  hydrodynamic  loads  are  calculated  by  an  external  subroutine  which  is  discussed  in  the 
next  section. 


7 A  Hydrodynamic  Loading  Module 

Every  time  step  the  ADAMS  solver  calculates  the  loading  at  every  point  in  the  system  several  times 
in  order  to  achieve  the  proper  prediction  of  system  response.  The  forces  on  the  blade  elements  and 
floater  sections  require  external  modules  to  be  called  to  calculate  the  loads  associated  with  them.  The 
floater  consists  of  6  components:  spar,  submerged  portion  of  tower,  and  four  spokes.  Each  of  these 
components  hydrodynamic  forces  are  calculated  in  fltr  loads  subroutine.  The  following  is  a 
description  of  the  calculations  perfonned  in  the  fltrjoads  subroutine. 

7.4.1  Initial  Calculations 

The  first  time  the  fltrjoads  subroutine  is  called  the  floater  input  file  is  read.  The  input  file  consists  of 
parameters  that  describe  the  sea  state  being  analyzed:  number  of  bands,  type  of  spectrum  (Pierson- 
Moskowitz  or  JONSWAP),  H1/3,  and  Tj.  The  number  of  bands  dictates  how  many  different 
frequencies  will  be  used  to  generate  the  sea  state. 

The  type  of  spectrum  is  chosen  based  on  the  geographic  location  where  the  floater  is  being  analyzed. 
The  JONSWAP  spectrum  is  useful  for  modeling  waves  in  the  north  Atlantic.  It  is  the  spectrum  used 
for  most  of  the  analyses  in  this  research.  As  shown  in  Figure  7.5,  the  JONSWAP  spectrum  has  a 
higher  peak  value  but  is  narrower  than  the  Pierson-Moskowitz  spectrum  for  a  given  sea  state.  The 
JONSWAP  spectrum  has  higher  peak  values  of  wave  elevation  due  to  this  difference  in  spectrum. 


Figure  7.5  Graphs  of  Pierson-Moskowitz  and  JONSWAP  spectrums  sea  state  5. 
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The  frequency  bands  are  described  by  a  wave  number,  amplitude,  and  phase.  The  phase  of  each  wave 
is  chosen  randomly  from  a  continuous  distribution  of  -180“  to  +180".  The  random  seed  used  for  this 
process  is  set  at  a  constant  value.  This  ensures  that  identical  frequency  bands  are  generated  every 
time  a  run  is  conducted  with  the  same  number  of  bands.  Changing  the  random  seed  changes  the 
phases  associated  with  these  bands.  The  spectrum  data  is  written  to  an  output  file.  A  time  history  of 
the  wave  elevation  at  the  origin  is  also  written  to  an  output  file.  Figure  7.6  is  an  example  of  a  wave 
elevation  record  for  a  sea  state  5  Pierson-Moskowitz  spectrum. 


-3 


0  100  200  300  400  500  600  700  800 

Time  (seconds) 

Figure  7.6  Record  of  wave  elevation  for  sea  state  5  with  30  bands. 

7.4.2  Hydrodynamic  Loads 

Every  time  step  the  solver  requests  the  forces  and  moments  on  the  six  components  of  the  floater.  The 
hydrodynamic  forces  consist  of  inertial  and  viscous  drag.  Horizontal  forces  on  the  submerged  portion 
of  the  tower  and  the  spar  create  moments  about  their  bases.  The  description  of  the  excitation  forces 
are  given  in  Chapter  6.  The  long  wavelength  approximation  is  assumed  to  be  true  since  the  diameter 
of  the  spar  is  much  less  than  the  expected  wavelengths.  The  structure  is  assumed  to  produce 
negligible  amounts  of  waves  due  to  its  relatively  small  diameter. 

The  spar  and  submerged  tower  calculations  are  very  similar.  There  are  two  differences.  The  first 
difference  is  the  heave  added  mass.  The  heave  added  mass  is  produced  related  to  the  surface  area  of 
the  bottom  of  the  spar.  Since  the  tower  is  connected  to  the  spar,  it  has  no  added  mass  in  heave. 


Therefore,  a  heave  single  added  mass  is  used  for  the  tower  and  spar  combination.  The  second 
difference  is  how  the  viscous  drag  is  calculated.  Drag  in  the  z-direction  is  assumed  to  be  zero  for  the 
submerged  portion  of  the  tower.  The  vertical  component  of  viscous  drag  is  produced  by  the  bottom 
of  the  spar. 

The  solver  sends  the  following  data  to  the  fltrjoads  subroutine:  time,  component  identification, 
component  positions,  component  velocities,  and  component  accelerations.  The  strip  method  is  used 
for  both  the  submerged  tower  and  the  spar.  The  position  of  each  is  used  to  calculate  the  horizontal 
and  vertical  velocity  and  acceleration  of  the  water  displaced  by  the  floater  strip.  These  velocities  and 
accelerations  are  then  used  to  calculate  the  initial  and  drag  forces  and  moments  using  GI  Taylor  (6.4) 
and  Morison  (6.9).  Summing  the  strips  contributions  results  in  the  total  force  (3  components)  and 
moments  (2  components).  No  yaw  moment  is  produced  by  hydrodynamic  loads  on  the  cylinders. 


The  spokes  have  a  veiy  small  contribution  to  the  overall  loads  on  the  system.  The  actual 
configuration  of  the  spokes  is  unknown.  The  spokes  could  be  open  trusses  or  closed  boxes  or 
cylinders.  For  ease  of  rnodeling  the  spokes  are  represented  by  closed  rectangular  boxes  with  a 
characteristic  dimension  (w). 


.  Wavier  Characteristics' 
and  Floater  Dimensions 
:  ■  tFioater  input  File)  f 


ADAMS  Marker 
Number  and  Motions 
Position,  Velocities, 
and  Accelerations 
(from  ADAMS) 


Jrilflal  Calculation  ■ v 

•  Wave.  Spectrum  .V I 

•  Frequenty  Bands. 


Figure  7.7  Flowchart  for  jltr  loads  subroutine. 


7.5  Complete  System  Model  Description 

The  floating  wind  turbine  system  is  modeled  by  coupling  the  dynamic  models  of  a  wind  turbine  and 
floater  using  the  ADAMS  mechanics  dynamic  analysis  program.  The  wind  turbine  model  is  made  up 
of  elements  that  are  connected  by  springs  and  dampers.  The  physical  values  for  the  restoring  and 
damping  coefficients  are  chosen  to  represent  the  mechanical  properties  of  the  wind  turbine  blades, 
drive  train,  and  tower.  Figure  7.8  shows  how  the  model  and  control  files  are  generated.  The 
ADAMS  preprocessor  in  FAST  {FastlAdams  subroutine)  and  the  floater  preprocessor  {fltr jprop 
subroutine)  act  together  to  write  the  files  necessary  to  run  the  fully  coupled  analysis. 


The  control  file  governs  the  operation  of  the  solver.  The  solver  determines  the  motion  and  total 
loading  of  the  elements  that  make  of  the  system.  The  aerodynamic  forces  acting  on  the  individual 
blade  elements  are  calculated  in  the  AeroDyn  subroutine.  AeroDyn  is  an  element-level  wind-turbine 
aerodynamics  analysis  routine.  It  requires  information  on  the  status  of  a  wind  turbine  from  the 
dynamics  analysis  routine  and  returns  the  aerodynamic  loads  for  each  blade  element  to  the  ADAMS 
Solver  [22]. 


The  input  wind  field  to  the  rotor  of  the  wind  turbine  is  a  data  file  that  contains  either  a  time  marching 
set  of  hub  height  wind  velocities  or  a  time  marching  three  dimensional  turbulent  wind  velocity 
profile.  A  program  produced  by  the  NWTC  called  SNWind  is  used  to  generate  the  wind  input  files. 
AeroDyn  uses  the  position  of  the  blade  element  to  determine  the  local  wind  flow  and  then  calculates 
velocity  vectors  for  that  element. 


As  discussed  in  Section  7.4  the  solver  also  uses  an  external  subroutine  to  calculate  the  hydrodynamic 
forces  on  the  floater  components.  The  overall  operation  of  the  solver  and  how  it  interacts  with  the 
external  subroutines  and  data  files  is  shown  in  Figure  7.9.  The  pitch  of  the  wind  turbine  blades  are 
regulated  by  a  pitch  control  subroutine.  Blade  pitch  is  changed  to  control  speed  of  the  rotor  or  power 
of  the  turbine. 


Figure  7.9  Analysis  of  model  flowchart. 


The  output  of  the  solver  consists  of  two  files.  The  first  is  a  graphics  file  that  contains  the  motion  of 
every  element  that  makes  up  the  system.  This  file  can  be  viewed  in  the  postprocessor.  The  second 
output  file  that  is  generated  is  the  request  file.  The  request  file  stores  all  additional  data  requested  by 
the  operator,  such  as  loads,  displacements,  velocities,  and  accelerations. 
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7.6  Model  Verification 


An  important  step  in  the  analysis  of  the  combined  wind  turbine  /  floater  system  is  to  ensure  that  the 
systems  are  being  modeled  correctly.  The  verification  of  the  system  as  a  whole  against  a  physical 
models  response  in  a  simulated  ocean  environment  would  be  ideal.  Unfortunately,  for  this  research, 
that  method  was  cost  prohibitive.  Instead,  the  individual  tools  used  in  the  analyses  are  verified 
independently  of  each  other. 

The  hydrodynamic  modeling  of  the  floater  is  performed  by  calculating  the  excitation  forces  on  the  six 
floater  elements  in  the  Jltr_loads  subroutine  and  integrating  them  into  the  system  dynamic  model  in 
the  solver.  The  verification  of  the  floater  modeling  is  accomplished  by  comparing  the  predicted 
floater  motion  of  a  simulation  with  the  analytical  solution  obtained  from  linear  theory. 

The  operation  of  the  floater  model  and  associated  hydrodynamic  loading  subroutine  was  tested  by 
analyzing  a  tension  leg  spar  buoy  without  a  wind  turbine  attached.  This  allowed  for  the  operation  of 
the  floater  loading  model  to  be  evaluated  without  interference  from  the  wind  turbine.  The  verification 
tests  were  conducted  for  two  different  wave  frequencies.  Both  of  these  test  cases  used  the  same 
floater  design  which  is  shown  below  in  Figure  7.10.  The  floater  input  file  is  included  in  Appendix  A. 
This  file  contains  a  detailed  description  of  the  structures  physical  properties. 

With  the  wind  turbine  removed  from  the  system  the  linear  dynamic  analysis  becomes  easier.  The 
system  shown  in  has  a  center  of  gravity  which  is  below  the  center  of  buoyancy.  The  distance  between 
the  center  of  gravity  and  the  center  of  buoyancy  greatly  affects  the  cross-coupling  between  surge  and 
pitch.  Table  7-1  lists  the  important  physical  characteristics  of  the  floater  used  in  the  verification  tests. 
A  time  domain  simulation  was  conducted  for  the  modified  system  response  in  1  meter  waves  with 
frequencies:  0.75  and  0.9  rads/sec. 
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Figure  7.10  Floater  design  used  in  verification  test. 
Table  7-1  Properties  of  floater  used  in  verification  tests. 


Water  Depth 


System  Mass 


Buoyant  Mass 


Reserve  Buoyancy 


Center  of  Gravity 


Center  of  Buoyancy 


2.53  cm 


50  mm 


Moments  of  Inertia 


451478.0 


The  results  of  the  surge  verification  tests  are  shown  in  Figure  7.1 1.  The  surge  response  of  the  floater 
generated  by  the  time  domain  simulation  using  the  hydrodynamics  loading  module  and  ADAMS 


solver  is  plotted  for  a  100  second  period  of  time  for  the  two  test  frequencies.  The  analytical 
predictions  for  the  system  motion  are  also  plotted.  The  results  of  the  time  domain  simulation  match 
the  expected  response  for  both  frequencies.  The  equation  of  motion  that  governs  the  surge  behavior 
is  given  by  (7.1).  Surge  is  strongly  coupled  with  pitch. 


^15  ^ 

4- 

'c„  c„' 

<  ^51  -^55  *^^55/ 

UJ 

UJ 

UJ 

The  motion  of  the  structure  in  pitch  is  also  given  by  equation  (7.1).  The  cross-coupling  between 
surge  and  pitch  was  very  evident.  When  the  As,  and  Ais  terms  were  set  to  zero  in  the  linear  analysis 
the  difference  between  the  predicted  results  and  simulations  became  very  large.  The  plots  of  the 
predicted  and  simulated  values  of  pitch  are  shown  in  Figure  7.12.  The  predicted  and  simulated  pitch 
values  are  the  same. 

The  agreement  between  the  simulated  and  predicted  response  of  both  modes  of  motion  verifies  the 
accuracy  of  the  method  by  which  the  floater  is  modeled.  Therefore,  the  motion  of  the  complete 
system  can  be  determined  using  the  time  domain  analysis  method  with  the  coupling  between  the 
floater  and  wind  turbine  being  handled  in  the  mechanical  dynamics  solver. 

The  wind  turbine  model  and  its  loading  and  control  subroutines  are  obtained  from  the  National  Wind 
Technology  Center.  The  prediction  of  aerodynamic  lift  and  drag  forces  is  a  difficult  problem. 
Computer  code  results  are  sometimes  very  different  from  wind  tunnel  data.  However,  AeroDyn  is 
seen  as  an  acceptable  method  for  determining  the  loads  on  blade  elements.  Therefore,  this  method  of 
modeling  wind  loading  on  the  rotor  blades  is  used  for  the  analysis  of  the  floating  wind  turbine  system. 

The  ADAMS/Solver  has  been  widely  used  by  industry  and  educational  organizations.  Its  operation 
has  been  fully  verified  for  several  different  mechanical  systems.  The  information  related  to  the 
verification  process  of  the  Solver  can  be  found  in  the  ADAMS  Verification  Guide  in  the  ADAMS 
User’s  Manual  [1]. 


Surge  verification  test  (Wave  Amplitude:  1m.  Wave  Frequency:  0,75  rad/sec) 
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Figure  7.1 1  Results  of  surge  verification  test. 


Figure  7.12  Results  of  pitch  verification  test 


7.7  Description  of  Wind  Turbine 

The  wind  turbine  used  in  this  analysis  was  a  1.5  MW  Horizontal  Axis  Wind  Turbine  (HAWT).  It  was 
a  three  bladed  upwind  design.  The  model  for  this  machine  was  obtained  from  the  NWTC  Design 
Codes  website.  It  is  referred  to  as  the  1 .5  MW  Baseline  wind  turbine.  This  turbine  was  chosen  for 
two  reasons.  First,  it  had  been  used  by  the  NWTC  in  the  past  for  evaluating  the  FAST  analysis  code. 
The  second  reason  for  the  selection  of  the  1.5  MW  Baseline  wind  turbine  was  it  was  the  only  model 
available  which  was  of  sufficient  size  to  warrant  mounting  it  on  a  floating  platform.  Most  currently 
installed  offshore  wind  turbines  are  of  tiiis  size  and  rating  or  larger.  The  principal  characteristics  of 
the  1 .5  MW  Baseline  turbine  are  given  in  Table  7-2. 

The  wind  turbine  model  consists  of  a  description  of  the  systems  mechanical,  aerodynamic,  electrical 
characteristics.  An  ADAMS  model  is  generated  by  the  FAST  preprocessor  based  on  the  input  files 
that  describe  the  wind  turbine  system.  A  FAST  input  file  for  the  1.5  MW  Baseline  wind  turbine  can 
be  found  in  APPENDIX  B.  The  other  information  required  to  describe  the  turbine  are  the  blade, 
tower,  and  AeroDyn  input  files.  The  blade  input  file  details  the  aerodynamic  properties  of  the 
different  airfoils  that  make  up  the  entire  turbine  blade.  The  tower  input  file  lists  the  physical 
characteristics  of  the  tower  section  masses,  elasticity,  and  positions.  The  AeroDyn  input  file  controls 
how  the  subroutine  calculates  aerodynamic  forces. 


Table  7-2  Characteristics  of  the  1.5  MW  Baseline  wind  turbine. 


Power  Output 

1.5  MW 

Nacelle  Mass 

51200  kg 

Hub  Height 

84.29  m 

Hub  Mass 

15100  kg 

Rotor  Diameter 

70  m 

Tower  Mass 

123003  kg 

Number  of  Blades 

3 

Number  of  Tower  Segments  in  Model 

10 

Rotational  Speed 

20  rpm 

Number  of  Blade  Elements  in  Model 

15 /blade 

Cut-in  Wind  Speed 

5  m/s 

Active  Pitch  Control 

yes 

Cutout  Wind  Speed 

25  m/s 

Type  of  Regulation 

speed  Regulate 

The  1.5  MW  Baseline  wind  turbine  model  can  be  modified  to  simulate  different  operating  conditions. 
User  written  controls  subroutines  can  be  substituted  for  the  default  controls  algorithms.  This  model  is 
a  useful  tool  in  analyzing  the  effect  of  placing  a  wind  turbine  on  a  floating  base  and  subjecting  the 
system  to  offshore  loads.  Figure  7.13  is  the  graphical  representation  of  the  wind  turbine. 
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8  Damping 


Damping  is  one  of  the  primary  characteristics  which  must  be  understood  for  a  moored  offshore 
structure.  The  floater  wind  turbine  system  has  three  means  of  damping  motions:  viscous  drag, 
changes  in  aerodynamic  forces,  and  structural  damping.  The  effects  of  these  damping  mechanisms 
are  analyzed  and  comparisons  made  between  them. 

The  damping  of  vertical  cylinders  in  an  ocean  environment  has  been  a  topic  of  concern  for  the 
offshore  industry  for  many  yearn.  The  translational  motion  of  the  system  is  of  the  same  order  of 
magnitude  as  the  diameter  of  the  floater.  The  viscous  drag  forces  on  the  floater  are  discussed  in 
section  6.3.  Viscous  drag  is  a  quadratic  damping  term. 

The  damping  effects  on  the  wind  turbine  arise  due  to  changes  in  inflow  of  air  to  the  wind  turbine 
when  the  system  moves.  Changes  in  airflow  across  the  blades  result  in  changes  in  thrust  and  torque 
on  the  rotor.  For  example;  the  relative  velocity  of  the  wind  into  the  turbine  increases  when  the  system 
surges  into  the  wind.  The  thrust  on  the  rotor  then  increases.  This  increased  thrust  resists  the  motion 
of  the  system  and  thereby  acts  as  a  damping  mechanism.  This  damping  mechanism  will  be  shown  to 
be  a  linear  damping  mechanism. 

Structural  damping  occurs  when  a  material  is  deformed  under  stress.  Some  of  the  energy  exerted  in 
the  deformation  is  absorbed  by  the  material.  The  motion  of  the  stracture  is  therefore  attenuated.  The 
effect  of  structural  damping  is  relatively  small  for  the  floating  wind  turbine.  In  order  to  determine  the 
damping  from  the  structure  all  of  the  other  damping  mechanism  must  be  secured.  A  single  surge  free 
decay  test  was  conducted  with  all  other  damping  mechanisms  secured.  The  results  showed  that  over  a 
4000  second  data  mn  the  amplitude  of  the  oscillation  changed  by  less  than  1  percent.  No  further 
efforts  were  made  to  understand  the  structural  damping. 

The  damping  of  the  system  for  various  modes  of  motion  is  foimd  by  conducting  free  decay  tests  on 
the  model.  These  tests  are  initiated  by  applying  a  loading  to  the  floater  to  cause  a  displacement. 
Then  the  load  is  removed  and  the  response  of  the  system  evaluated.  Thirty  two  damping  runs  where 
conducted  with  and  without  the  wind  turbine  operating.  In  some  tests  viscous  drag  was  set  to  zero  so 
that  the  wind  turbine  could  be  evaluated  separately  from  the  floater.  The  results  of  the  damping  tests 
are  detailed  in  the  following  sections. 


Faltinsen  described  a  method  for  evaluating  the  results  of  free  decays  to  extract  the  linear  and 
quadratic  damping  coefficients  (p.  252,  [13]).  This  method  is  normally  used  to  physical  tests.  It  can 
be  applied  to  the  simulations  that  involve  more  than  one  damping  mechanisms.  Using  this  technique 
can  enable  single  free  decay  data  mns  to  be  used  to  exact  both  damping  coefficients  for  the  system. 
The  equation  of  motion  for  a  body  undergoing  free  decay  is: 


X  +  p^X  +  X  JC  +  PjX  =  0 


(8.1) 


where  pi  is  the  linear  damping  coefficient  and  p2  is  the  quadratic  coefficient.  If  the  damping  is 
constant  with  the  amplitude  of  oscillation  the  damping  coefficients  can  be  found  using  the  following 
relationship: 
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(8.2) 


where  Xn  is  the  amplitude  of  the  nth  oscillation.  The  oscillations  are  separated  by  a  time  equal  to  half 
the  period.  The  left  side  of  equations  (8.2)  is  plotted  as  the  dependant  quantity  of  term  16X„  / 

The  points  are  then  fitted  with  a  straight  line  using  the  least  squares  method.  The  y-intercept  and 
slope  of  the  fitted  line  correspond  to  pt  and  p2  coefficient. 


This  method  allows  for  the  non-dimensional  coefficients  p,  and  p2  to  be  determined  for  a  time  record 
of  a  motion  that  damps  out  with  time.  These  coefficients  are  useful  when  discussing  the  relative 
importance  of  a  damping  mechanism  (viscous  drag  or  wind  turbine)  for  a  given  mode  of  motion.  In 
some  modes  of  motion  viscous  drag  is  found  to  dominate  and  in  others  the  linear  damping  of  the  wind 
turbine.  These  coefficients  are  also  useful  when  trying  to  extract  the  effect  individual  parameter 
changes  (such  as  blade  pitch)  have  on  system  damping. 
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8. 1  Description  of  Fioating  Wind  Turbine  used  in  Damping  Tests 

A  single  floater  design  was  used  for  all  of  the  damping  tests.  It  consisted  of  both  a  spar  cylinder  and 
an  extension  of  the  tower  below  the  surface.  The  floater  supported  the  1 .5  MW  wind  turbine.  Table 
8-1  gives  the  dimensions  and  characteristics  of  the  floating  support. 


Table  8-1  Dimensions  and  properties  of  floater  used  in  damping  tests. 


Spar  Diameter 

10m 

Line  Diameter 

112  mm 

Spar  Height 

12m 

Number  of  Lines 

2 / spoke 

Tower  Diameter 

8m 

Line  Minimum  Break  Load 

371  tonnes 

Tower  Draft 

5  m 

Maximum  Line  Extension 

14% 

Spoke  Length 

20  m 

Line  Length  (un-stretched) 

176.2  m 

Spoke  Width 

1  m 

Water  Depth 

200  m 

Buoyant  Spokes 

No 

Structure  Thickness 

2.53  cm 

The  lines  used  in  the  damping  test  were  chosen  for  their  flexibility.  Having  flexible  lines  allowed  for 
larger  amplitudes  of  oscillation  during  the  free  decay  teste.  Polyester  lines  have  elongation  properties 
similar  to  those  of  the  damping  floater  model.  The  calculated  properties  of  the  system  are  given  in 
Table  8-2.  The  natural  periods  are  included  for  the  sole  purpose  of  comparing  them  to  the  simulation 
results.  A  rendering  of  the  floating  wind  turbine  model  is  shown  in  Figure  8.1.  No  wave  excitation  is 
used  during  the  damping  teste. 


Table  8-2  Calculated  attributes  of  the  floating  wind  turbine  used  in  the  damping  teste. 


Mass  of  Floater 

Center  of  Gravity  (system) 

19.608  m 

Buoyant  Mass 

1228424.6  kg 

Center  of  Buoyancy 

-9.211m 

Total  Mass  of  System 

BG 

28.819  m 

Reserve  Buoyancy 

198.5  % 

Surge  /  Sway  Natural  Period 

39.02  sec 

Initial  Line  Stress 

101.7  MPa 

Heave  Natural  Period 

4.45  sec 

Breaking  Line  Stress 

369.3  MPa 

Pitch  and  Roll  Natural  Period 

10.82  sec 

Center  of  Gravity  (floater) 

-11.82  m 

Yaw  Natural  Period 

10.52  sec 

Figure  8.1  Floating  wind  turbine  system  used  in  damping  tests. 

8.2  Surge  Damping 

The  major  contributor  to  ocean  wave  generation  is  wind  acting  on  the  surface  of  the  water.  Wind 
generated  waves  tend  to  propagate  in  the  direction  of  the  winds.  In  a  steady  wind  enviroiunent,  with 
winds  blowing  in  the  positive  x-direction,  the  wave  motion  of  the  floating  wind  turbine  will  be  in  the 
surge  direction.  Therefore,  surge  damping  is  one  of  the  most  important  properties  to  understand. 

The  free  decay  test  for  the  surge  damping  is  conducted  by  applying  a  load  in  the  x-direction.  This 
load  causes  a  deflection  in  the  floater  away  from  its  rest  position.  Damping  is  determined  with  no 
waves  acting  on  the  structure.  Five  tests  were  run  for  the  surge  damping  analysis.  When  the  wind 
turbine  is  operating  wind  is  assumed  to  be  steady  at  20  m/s.  For  the  first  test  the  wind  turbine  was 
secured  so  that  the  damping  of  the  viscous  drag  could  be  evaluated.  The  next  three  tests  were  run  with 
the  viscous  drag  set  to  zero  for  the  floater.  These  three  tests  evaluated  the  effect  of  blade  pitch  on  the 


system  damping.  The  blade  pitch  of  the  wind  turbine  was  set  at  three  different  values;  18  degrees,  22 
degrees,  and  26  degrees.  The  final  test  consisted  of  the  wind  turbine  operating  with  viscous  drag  on 
the  structure. 


8.2.1  Floater  Damping 

The  damping  properties  of  the  floater  were  analyzed  by  performing  a  surge  free  decay  simulation 
with  the  wind  turbine  secured.  The  results  of  this  damping  test  are  shown  in  Figure  8.2.  There  are 
two  graphs  of  the  response  to  the  surge  impulse.  The  viscous  drag  attenuates  the  amplitude  over  1250 
seconds.  The  bottom  graph  is  the  absolute  value  of  the  surge  amplitude  on  a  logarithmic  scale.  This 
graph  shows  that  viscous  drag  is  a  quadratic  damping  mechanism,  since  the  decay  is  greater  than 
linear. 


The  quadratic  damping  coefficient  for  the  system  is  0.7.  It  should  be  noted  the  frequency  of  the  surge 
response  during  the  free  decay  test  is  39.4  seconds.  Which  corresponds  to  the  natural  fi'equency 
predicted  using  linear  theory.  The  output  of  this  test  was  evaluated  using  the  method  described  by  at 
the  end  of  the  previous  section.  Viscous  damping  results  in  a  rapid  attenuation  of  the  oscillations 
early  on.  This  damping  force  becomes  very  small  once  the  motion  of  the  floater  slows,  since  viscous 
damping  is  proportional  to  the  square  of  the  floater  velocity. 


The  values  for  the  surge  damping  of  the  system  with  the  wind  turbine  secured  were  found  by  fitting  a 
curve  through  the  data  points  in  Figure  8.3.  The  value  for pi  was  found  to  be  very  small  negative 
number.  This  indicates  that  linear  damping  was  negligible  during  the  damping  test.  When  the  wind 
turbine  was  secured  the  system  experienced  small  oscillations  in  modes  other  that  surge  and  pitch 
during  this  simulation.  These  vibrations  resulted  in  data  scatter  which  is  responsible  for  the  negative 
value  for  pi.  Ideally,  this  value  would  be  zero.  The  p2  coefficient  was  determined  to  be  0.0392. 
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Figure  8.3  Plot  used  to  determine  p,  and  p2  for  the  surge  damping  test.  Wind  turbine  secured. 

8.2.2  Damping  due  to  Wind  Turbine 

The  damping  properties  of  the  wind  turbine  were  evaluated  by  setting  the  viscous  drag  of  the  floater 
to  zero.  The  damping  of  the  wind  turbine  was  found  to  be  linear  as  expected.  The  pitch  of  the  blades 
was  found  to  affect  pitch  damping  of  the  wind  turbine.  Free  decay  simulations  were  conducted  for 
three  settings  of  blade  pitch:  18  degrees,  22  degrees,  and  26  degrees.  The  results  of  the  free  decay 
tests  are  shown  in  Figure  8.4.  It  should  be  noted  that  when  in  operation  the  wind  turbine  causes  a 
mean  displacement  of  the  floater  due  to  the  thrast  on  the  rotor. 

The  wind  turbine  was  shown  to  have  a  stabilizing  effect  on  the  system  due  to  its  linear  damping 
characteristics.  Small  oscillations  which  were  present  for  the  surge  free  decay  test  with  the  wind 
turbine  secured  were  not  seen  when  the  wind  turbine  was  operating. 


Figure  8.4  Surge  damping  results  for  wind  turbine  with  blade  pitch  equal  set  at  1 8, 22,  and  26 
degrees.  20  m/s  constant  winds.  No  viscous  drag. 


Changing  the  blade  pitch  changes  the  angle  of  attack  of  the  airfoil.  Lift  and  drag  of  an  airfoil  is 
dependant  on  angle  of  attack.  Therefore,  the  amount  of  thrust  for  a  given  change  in  inflow  to  the 
rotor  changes  with  blade  pitch.  This  is  the  reason  that  the  damping  characteristics  of  the  three  tests 
differ.  In  order  to  extract  the  linear  damping  coefficients  of  the  mean  offsets  are  removed  and  the 
absolute  values  of  the  surge  oscillations  plotted  on  a  logarithmic  scale. 
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Figure  8.5  Absolute  value  of  surge  response  for  free  decay  tests  with  18, 22,  and  26  degree  blade 
pitch.  20  m/s  constant  winds.  No  viscous  drag. 

Evaluation  of  Figure  8.5  yields  the  linear  damping  coefficients  for  the  three  tests.  When  there  is  no 
quadratic  damping  the  p2  in  the  equation  that  governs  a  linear  damping  response  is: 


2  ,  X„, 


Where  X(t)  is  the  amplitude  of  oscillation  at  time  t  and  Xo  is  the  initial  displacement.  A  more  useful 
form  of  this  equation  allows  for  multiple  cycles  to  be  used  in  a  single  calculation: 
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where  is  the  number  of  half  cycles.  The  linear  damping  coefficients  for  the  three  different  blade 
pitches  are  given  in  Table  8-3.  The  damping  coefficient  decreases  with  increasing  blade  pitch. 
Changing  the  blade  pitch  also  affected  the  power  output  of  the  wind  turbine.  This  is  due  to  the  fact 
that  thmst  and  torque  on  the  rotor  are  both  governed  by  inflow  speed  of  the  wind  and  the  blade  pitch 


Table  8-3  Damping  coefficients  for  the  wind  turbine  with  three  different  blade  pitch  angles:  18 
degrees,  22  degrees,  and  26  degrees.  20  m/s  constant  winds.  No  viscous  drag. 


Blade  Pitch 

Damping  Coefficient  {pi) 

18° 

8.37  E-3 

22° 

6.90  E-3 

26° 

5.1  E-3 

8.2.3  Surge  Response  with  both  Linear  and  Quadratic  Effects. 

The  final  free  decay  test  involved  testing  the  floating  wind  turbine  system  as  a  whole.  The  test 
simulated  the  damping  response  of  the  system  operating  near  the  cutout  wind  speed  for  the  turbine. 
The  wind  turbine  was  operating  with  20  m/s  constant  winds  into  the  rotor.  The  pitch  of  the  rotor 
blades  was  fixed  at  22  degrees.  The  second  test  simulated  the  damping  response  of  the  system 
operating  in  12  m/s  winds.  In  normal  operation  the  pitch  of  the  rotor  blades  are  actively  controlled  to 
either  maintain  rotor  speed  or  maximize  generator  power.  For  this  second  test  the  wind  turbine  was 
operating  in  speed  regulate,  with  the  control  system  maintaining  the  rotor  at  20  rpm.  Viscous  drag 
effects  of  the  floater  were  also  included.  These  tests  model  the  response  of  the  system  undergoing  an 
impulse  load  in  the  direction  of  the  positive  x-direction. 


Figure  8.6  is  a  set  of  graphs  of  the  surge  response  to  the  firee  decay  test  for  the  wind  turbine  operating 
in  20  m/s.  By  comparing  these  graphs  to  Figure  8.2  and  Figure  8.4  is  can  be  seen  that  the  combined 
effect  of  the  two  modes  of  damping  greatly  decrease  the  time  it  takes  the  system  motion  to  damp  out. 
The  results  graphed  on  a  log  scale  in  the  following  figure  shows  how  the  damping  is  greater  than 
linear  at  the  beginning  of  the  test  when  quadratic  effects  dominate.  At  the  end  of  the  test  the  linear 
damping  from  the  wind  turbine  dominates. 
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Figure  8.6  Surge  free  decay  test  with  both  linear  and  quadratic  damping  effects.  Wind  turbine 

operating  in  20  m/s  constant  winds. 

The  coefficients  pi  and  p2  are  found  by  evaluating  the  second  graph  using  the  method  previously 
described.  The  terms  on  the  left  and  right  side  of  equation  (8.2)  are  plotted  against  each  other  and  a 
linear  regression  performed  on  the  data  points.  Figure  8.7  is  the  plot  used  to  determine  the  damping 
coefficients.  The  linear  damping  coefficient  {pi)  is  found  to  be  0.0052  and  the  quadratic  damping 
coefficient  {pi)  is  found  to  be  0.0486. 
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Figure  8.7  Determining  coefficients pi  and  p2  for  the  surge  free  decay  test  with  all  mechanisms 

included.  20  m/s  constant  winds. 

The  combined  system  had  linear  and  quadratic  damping  coefficients  which  were  similar  to  the 
individual  tests.  The  differences  were  not  large.  The  difference  in  the  quadratic  damping  term  p2  is 
attributable  to  data  scatter  in  the  results  of  surge  response  with  viscous  drag  only.  Overall,  the 
combination  of  the  two  damping  mechanisms  appears  to  match  the  results  expected  from  the  single 
damping  mechanism  teste.  All  remaining  tests  include  both  viscous  drag  and  linear  damping  from  the 
wind  turbine. 

8.3  Sway  Damping 

In  this  research  the  1.5  MW  Baseline  wind  turbine  is  a  HAWT  with  its  nacelle  aligned  with  the  x- 
axis.  Therefore,  sway  motion  is  the  horizontal  motion  of  the  floater  in  the  plane  of  the  rotor.  This 
type  of  motion  will  have  a  much  less  significant  effect  on  inflow  to  the  turbine  than  surge  motion. 
Therefore,  the  damping  effect  of  the  wind  turbine  is  expected  to  be  lower  in  the  sway  mode  of 
motion. 

A  simulated  free  decay  test  was  conducted  by  loading  the  model  in  the  positive  y~direction  and  then 
removing  the  load.  The  response  of  the  model  was  then  generated.  The  resulte  of  the  sway  free 
decay  test  were  plotted  in  Figure  8.8.  The  mean  offset  in  the  y-direction  was  caused  by  torque  on  the 


rotor. 
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Figure  8.8  Sway  free  decay  test  with  both  linear  and  quadratic  damping  effects.  Wind  turbine 

operating  in  20  m/s  constant  winds. 


Next  the  coefficients  for  the  linear  and  quadratic  damping  terms  were  determined.  The  values  found 
for  Pi  and  p2  are  0.0454  and  0.0015  respectively.  Figure  8.9  is  a  plot  of  the  data  points  used  to 
determine  the  two  coefficients.  The  quadratic  damping  coefficient  pt  is  approximately  the  same  as 
the  value  found  for  the  surge  damping.  This  is  expected  since  the  floater  is  moving  in  sway  is  nearly 
identical  to  the  floater  moving  in  sway. 

The  linear  damping  coefficient  in  sway  was  found  to  be  approximately  30  percent  of  the  surge 
coefficient  for  the  same  initial  conditions.  The  sway  linear  damping  coefficient  is  less  than  the  surge 
coefficient  since  the  motion  of  the  structure  does  not  affect  the  wind  inflow  to  the  wind  turbine  as 
much  as  surge  motion.  This  difference  in  linear  damping  resulted  in  the  sway  oscillations  taking 
longer  to  damp  out  than  those  of  the  surge  test. 


81 


Figure  8.9  Determining  pi  and  p2  for  the  sway  free  decay  test. 

8.4  Heave  Damping 


Vertical  motion  of  the  floating  wind  turbine  is  damped  out  by  a  combination  of  changes  in  wind 
turbine  loading  and  viscous  drag  on  the  floater.  The  natural  frequency  for  the  floating  wind  turbine  in 
heave  is  higher  than  for  the  other  two  modes  of  translational  movement  (surge  and  sway).  The 
response  of  the  system  to  the  heave  free  decay  test  is  therefore  characterized  by  short  periods  with 
small  amplitudes  of  oscillation. 


In  order  to  determine  the  heave  damping  coefficients  a  free  decay  test  was  conducted  with  die  wind 
turbine  operating  in  20  m/s  constant  winds.  The  blade  pitch  was  fixed  at  22  degrees.  The  floater 
model  was  loaded  in  the  positive  z-direction  and  then  released.  The  simulated  results  were  then 
determined.  The  frequency  of  the  damping  response  was  4.44  seconds  which  agreed  with  the  period 
predicted  using  linear  theory  of  4.45  seconds.  The  heave  position  was  plotted  with  respect  to  time  in 
Figure  8.10. 


m 
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Figure  8.10  Heave  free  decay  test  with  both  linear  and  quadratic  damping  effects.  Wind  turbine 
•  operating  in  20  m/s  constant  winds. 

The  linear  and  quadratic  coefficients  were  determined  from  the  free  decay  results.  The  data  points 
used  in  this  determination  were  plotted  in  Figure  8.11.  Linear  regression  of  the  data  resulted  in  a 
value  of  .0039  for  p,  and  .244  for  p2.  The  quadratic  damping  (p^)  due  to  viscous  drag  is  relatively 
large  in  heave  mode  of  motion.  The  linear  drag  coefficient  is  same  order  of  magnitude  as  that  of  the 
sway  mode  of  motion.  Since  heave  and  sway  motion  a  both  moving  the  rotor  nonnal  to  the  wind 
■flow  is  anticipated  that  the  damping  effect  of  the  wind  turbine  be  similar  for  the  two  modes  of  motion. 
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Figure  8,1 1  Determining  pj  and  p2  for  the  heave  free  decay  test. 

8.5  Roll  Damping 

The  three  rotational  modes  of  motion  (roll,  pitch,  and  yaw)  are  coupled  to  other  modes  of  motion  by 
hydrodynamic  effects  or  by  gyroscopic  moments  generated  when  the  rotor  is  moved  from  its  plane  of 
rotation.  These  cross-coupling  effects  make  it  more  difficult  to  extract  the  damping  characteristics  of 
the  floating  wind  turbine  system  than  the  translational  modes  of  motion.  Translational  motions  in  the 
horizontal  directions  (surge  and  sway)  have  relatively  large  magnitudes.  The  cross-coupling  from  the 
rotational  motions  (pitch  and  sway)  do  not  cause  large  changes  in  the  horizontal  motions.  However, 
the  cross-coupling  between  the  translational  modes  to  the  rotational  modes  does  have  a  large  effect  on 
the  rotational  motion  of  the  system  in  pitch  and  sway.  The  result  of  this  behavior  is  that  extracting 
the  damping  coefficiente  for  these  rotational  modes  is  more  difficult 

The  first  rotational  motion  that  was  examined  was  roll.  In  roll  the  structure  pivots  around  an  x-axis. 
This  motion  did  not  change  the  plane  of  rotation  of  the  rotor.  So  gyroscopic  loads  where  not 
generated  by  the  roll  motion.  The  roll  motion  was  cross-coupled  to  the  sway  mode  of  motion. 
Oscillating  displacements  in  roll  caused  loading  hydrodynamic  loading  that  resulted  in  sway  motion. 
The  sway  motion  altered  the  roll  oscillations  of  the  structure  during  flie  free  decay  test.  The  roll 
position  did  not  oscillate  about  a  fixed  point  during  the  free  decay  test  (see  Figure  8.17).  The  period 
for  the  roll  response  was  10.9  seconds  which  is  comparable  to  the  predicted  value  of  10.8  seconds. 
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Figure  8.12  Roll  damping  free  decay  test  Wind  turbine  operating  in  20  tn/s  winds.  Blade  pitch  fixed 

at  22°. 

The  interaction  between  roll  and  sway  can  be  seen  on  the  logarithmic  plot  of  the  data.  The  absolute 
values  of  amplitude  do  not  have  a  smooth  trend.  The  result  of  this  behavior  is  data  scatter  during  the 
determination  of  the  p/  and  p2.  Figure  8.13  shows  the  data  points  extracted  from  the  free  decay  test. 
The  values  for  the  linear  and  quadratic  damping  coefficients  are  0.0285  and  0.0015  respectively. 


Figure  8.13  Determining  pi  and  p^  for  roll  free  decay  test. 

85 


8.6  Pitch  Damping 

The  floater  rotates  about  a  y-axis  when  it  pitches.  This  motion  acts  to  move  the  rotor  directly  into  and 
away  from  the  wind  flow,  assuming  that  the  nacelle  and  wind  both  are  aligned  with  the  x-axis.  It  is 
intuitively  obvious  that  the  damping  affect  of  the  turbine  will  be  very  large  when  compared  to  most  of 
the  other  modes  of  motion.  This  is  due  to  the  fact  that  these  other  modes  don’t  cause  lai^e  changes  in 
the  inflow  to  the  turbine. 

The  damping  effect  of  the  wind  turbine  was  so  large  that  it  made  determining  the  linear  and  quadratic 
damping  coefficients  difficult.  As  shown  in  Figure  8.14,  the  pitch  motion  of  the  floater  due  to  the 
free  decay  test  rapidly  damped  out.  After  only  four  cycles  the  pitch  motion  was  dominated  by  cross¬ 
coupling  with  the  surge  mode  of  motion  (frequency  of  motion  changed  to  natural  frequency  of  surge.) 

Pitch  free  decay  test. 


Time  (sec) 


Figure  8.14  Pitch  free  decay  test.  20  m/s  constant  winds.  Blade  pitch  fixed  at  22  degrees. 

The  pitch  damping  due  to  the  wind  turbine  was  so  strong  that  only  four  data  points  could  be  obtained 
to  perform  the  least  squares  procedure  (see  Figure  8.15).  This  number  is  less  than  optimal  for 
determining  the  damping  coefficients.  The  strong  linear  damping  of  the  system  prevented  an  accurate 
measurement  of  the  viscous  damping  term.  The  values  for  the  damping  coefficients  found  from  the 
linear  regression  technique  were:  pj  =  0.1665  and  =  -0.0195.  The  value  for  pi  is  negative.  We 
know  that  the  viscous  drag  coefficient  should  be  equal  to  that  of  roll  since  the  system  is  symmetrical 
and  should  have  similar  viscous  damping  properties.  Therefore,  it  is  assumed  that  the  value  for  p2  is 
0.001 5  to  match  the  roll  viscous  drag. 


Figure  8.15  Determining  p,  and  p2  for  pitch  free  decay  test. 

8.7  Yaw  Damping 

Yaw  damping  of  the  system  consists  of  both  linear  and  quadratic  effects.  The  linear  damping  occurs 
when  a  moment  is  produced  on  the  rotor  opposing  the  yaw  motion  due  to  changing  aerodynamic 
forces.  As  the  system  yaws,  the  relative  wind  speed  increases  on  one  side  of  the  rotor  and  decrease 
on  the  other.  The  difference  in  relative  wind  speeds  results  in  a  net  moment  in  the  opposite  direction 
of  the  yaw  motion.  The  viscous  damping  comes  from  the  spokes.  As  the  floater  yaws,  drag  forces 
develop  on  the  spokes.  These  forces  act  to  counter  the  system  rotational  motion. 

It  was  difficult  to  conduct  a  yaw  free  decay  simulation  on  the  design  selected  for  the  damping  tests. 
When  the  system  was  displaced  and  then  released,  it  did  not  simply  oscillate  in  the  yaw  mode  (see 
Figure  8.16).  Observation  of  other  modes  of  motion  showed  that  the  system  was  being  excited  into 
other  modes  of  motion  by  the  yaw  motion.  The  frequency  of  the  oscillations  was  also  not  constant. 
The  yaw  damping  loads  are  very  large  and  probably  the  root  cause.  This  behavior  is  exaggerated 
greatly  when  the  wind  turbine  is  secured  and  a  free  decay  test  conducted  (see  Figure  8.17).  The 
cross-coupling  between  roll  and  sway  can  be  seen  by  observing  the  tail  end  of  the  roll  curve.  The 
excited  roll  motion  damps  out  (shorter  period)  but  the  structure  continues  to  roll  due  to  sway  motion 
(longer  period). 
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Figure  8.16  Yaw  free  decay  test  with  system  free  to  move  in  all  modes. 


Yaw  and  roll  response  during  yaw  free  decay  test. 


Time  (sec) 


Figure  8.17  Graph  of  yaw  and  roll  for  free  decay  test  with  wind  turbine  secured. 


The  problems  with  determining  the  yaw  damping  coefficients  were  overcome  by  isolating  the  system 
so  that  it  only  moved  in  the  yaw  mode  of  motion.  This  was  accomplished  by  creating  a  joint  in  the 
bottom  of  the  floater  that  attached  it  directly  to  ground.  This  joint  only  allowed  motion  about  the  z- 
axis.  A  free  decay  test  was  conducted  and  the  damping  coefficients  found  for  the  modified  model 
(Figure  8.18).  The  motion  of  the  system  damped  out  rapidly  due  to  the  large  moments  generated  by 
the  spokes  and  rotor.  The  period  of  the  system  response  was  found  to  be  10.4  seconds  which  is  close 
to  the  period  predicted  by  linear  theory  (10.5  seconds). 
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Figure  8.18  Free  decay  response  with  only  yaw  motion.  Wind  turbine  operating  in  20  m/s  winds. 


The  wind  turbine  dominated  the  yaw  damping  of  the  system.  This  is  evident  by  observing  how  the 
absolute  value  of  yaw  response  decays  with  time  in  Figure  8.18.  The  decay  appears  to  be  nearly 
linear  on  the  log  plot,  which  is  the  expected  behavior  for  a  system  with  a  strong  linear  damping 
mechanism.  The  values  for  pi  and  p2  were  extracted  from  the  response  record  of  the  modified  system 
(Figure  8.19).  The  values  were  found  to  be  .0715  for  p ,  and  .0327  for  p2. 


■log 


X. 


X. 


rt+1 


iii 

3  T 


Figure  8*19  Determining  pi  and  p2  for  yaw  free  decay  test* 

A  mechanism  for  coupling  yaw  motion  of  die  system  to  other  modes  of  motion  exists*  A  gyroscopic 
moment  is  produced  about  the  x-axis  when  the  spinning  rotor  is  forced  to  change  its  orientation.  This 
excites  a  response  in  the  roll  mode*  The  equation  governing  this  behavior  is  discussed  in  3.2*1.  A 
separate  yaw  test  was  run  widi  40  meter  spokes  instead  of  20  meters  to  examine  this  behavior.  Using 
this  modified  floater  design  allowed  for  yaw  tree  decay  tests  to  be  conducted  without  the  need  to 
isolate  the  floater  motion  to  die  yaw  mode  only*  The  effect  of  gyroscopic  cross-coupling  between 
yaw  and  roll  was  simulated.  As  shown  in  Figure  8.20  this  coupling  mechanism  is  relatively  weak. 


Figure  8.20  Yaw  and  roll  cross-coupling  during  yaw  free  decay  test. 
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8.8  Surge  and  Sway  Damping  in  Normai  Operating  Conditions 

For  the  damping  tests  previously  described  the  wind  turbine  was  assumed  to  operate  in  a  wind 
environment  near  the  cutout  speed  for  the  machine.  This  is  the  type  of  environment  that  will  generate 
extreme  wave  events  and  therefore  the  associated  damping  characteristics  of  the  system  are  very 
important.  The  normal  operating  conditions  of  the  wind  turbine  are  winds  in  the  region  of  5-12  m/s. 
The  blade  pitch  at  these  wind  speeds  will  be  much  lower  than  the  22°  at  which  most  of  the  previous 
damping  coefficients  were  determined. 

As  discussed  earlier  blade  pitch  affects  the  damping  of  the  system.  The  surge  damping  coefficient 
was  shown  to  increase  as  blade  pitch  was  reduced  (see  Table  8-3).  However  the  effect  this  change  in 
blade  pitch  had  on  the  sway  damping  coefficients  was  not  discussed.  The  damping  characteristics  of 
the  floating  wind  turbine  in  the  surge  and  sway  modes  of  motion  were  evaluated  in  order  to  determine 
the  system  damping  in  the  normal  operating  region.  Free  decay  tests  were  performed  on  the 
simulation  model  in  both  surge  and  sway  with  hub  height  wind  speeds  of  7  m/s.  The  control  system 
maintained  the  blade  pitch  remained  fixed  at  2.6°  throughout  both  tests  since  this  was  the  optimum 
angle  for  energy  absorption  for  this  wind  speed. 

Due  to  flexibility  and  height  of  the  tower  the  motion  of  the  nacelle  may  not  coincide  with  that  of  the 
floater.  The  nacelle  impulse  responses  for  the  two  tests  are  included  with  the  floater  surge  and  sway 
responses.  Figure  8.21  is  a  graph  of  the  two  free  decay  test  results.  The  motion  of  theiiiacelle  follows 
that  of  the  floater  in  both  tests.  The  initial  amplitudes  of  the  nacelle  surge  and  sway  response  are 
larger  than  that  of  the  floater  due  to  the  effects  of  pitch  and  sway  moving  the  tower  top. 

It  should  be  noted  that  the  offset  between  the  nacelle  and  floater  surge  position  is  due  to  the  system 
having  a  mean  pitch  of  1.5°  caused  by  the  wind  turbine  loading.  A  close  observation  of  the  sway 
response  for  the  nacelle  (lower  graph)  reveals  small  oscillations  superimposed  on  the  larger  sway 
motion.  These  smaller  oscillations  have  a  periodicity  of  10.8  seconds  and  correspond  to  the  roll 
motion  of  the  system.  There  is  no  corresponding  motion  in  the  surge  graphs  due  to  the  fact  that  the 
pitch  damping  of  the  wind  turbine  is  so  strong  that  the  system  only  experienced  negligible  pitch 
motion  (  <.25°)  during  the  free  decay  test.  The  effects  of  pitch  are  only  seen  in  the  first  two 
oscillations  of  the  surge  free  decay  test  and  this  motion  is  1 80°  out  of  phase  with  the  surge  motion 
(offset  between  floater  and  nacelle  became  smaller.)  The  damping  coefficients  for  these  tests  are 
given  in  Table  8-4.  The  damping  coefficients  for  the  nacelle  in  sway  could  not  be  determined  due  to 
interference  from  roll  mode  of  motion. 


Figure  8.21  Surge  and  sway  free  decay  response  for  system  operating  with  7  m/s  winds.  Both  floater 


and  nacelle  positions  are  plotted. 


Table  8-4  Surge  and  sway  damping  coefficients  for  normal  operating  conditions.  Hub  height  wind 

speed  =  7  m/s.  . 


8.9  Summary  of  Single  Mode  Damping  Tests 

The  free  decay  tests  for  the  individual  modes  of  motion  showed  how  the  floating  wind  turbine  system 
damped  out  excited  motions.  Over  fifty  free  decay  tests  were  conducted  in  the  process  of  evaluating 
the  system  model  and  determine  the  damping  coefficients.  Not  all  of  the  modes  could  be  completely 
isolated  from  each  other.  For  example  the  surge/pitch  and  sway/roll  interactions.  When  one  of  these 
modes  was  excited,  the  cross-coupled  mode  was  also  excited.  The  cross-coupling  made  some  of  the 
data  evaluations  difficult  since  there  were  changes  in  the  oscillation  of  the  motion  being  analyzed. 

Two  tests  could  not  be  evaluated  fully  to  determine  either  one  or  both  of  the  associated  quadratic  and 
linear  damping  coefficients.  The  linear  pitch  damping  mechanism  was  very  strong  and  attenuated  the 
pitch  motion  after  only  a  few  cycles.  This  fact  combined  with  the  cross-coupling  with  surge 
hampered  the  data  collection  for  the  pitch  free  decay  test.  While  the  wind  turbine  was  operating,  only 
a  rough  estimate  of  the  pitch  linear  damping  coefficient  could  be  obtained.  The  floater  is  identical 
with  respect  to  roll  and  pitch  and  therefore  the  roll  quadratic  damping  coefficient  was  also  assigned  to 
the  pitch  mode  of  motion. 

For  the  design  chosen  for  the  damping  tests,  it  was  not  possible  to  perform  a  yaw  free  decay  test  on 
the  system  with  the  system  free  to  move  in  other  modes  of  motion.  The  system  was  constrained  to 
oscillate  in  the  yaw  mode  only  and  a  free  decay  test  conducted.  Coefficients  were  then  determined 
for  the  yaw  mode  of  motion.  In  later  tests,  longer  spokes  were  found  to  make  the  system  oscillate  in 
yaw  in  a  stable  manner  which  would  have  allowed  free  decay  data  to  be  gathered. 

It  is  difficult  to  extract  direct  comparisons  between  the  linear  and  quadratic  damping  coefficients. 
However,  some  general  characterizations  can  be  made.  For  a  system  configured  as  in  the  test  case  the 
following  appears  to  be  true.  The  wind  turbine  was  found  to  have  a  large  effect  on  the  systems 
damping  properties.  The  wind  turbine  was  the  major  damping  mechanism  for  both  pitch  and  yaw. 
Viscous  drag  was  the  major  damping  mechanism  in  the  translational  modes  of  motion;  surge,  sway, 
and  heave. 

The  wind  turbine  had  a  damping  effect  even  when  it  was  moved  across  its  plane  of  rotation  (roll  and 
sway).  The  free  decay  tests  where  wind  speed  was  20  m/s  the  sway  linear  damping  coefficient  was 
30  percent  of  the  surge  mode  of  motion.  The  roll  linear  damping  coefficient  was  determined  to  be 
roughly  25  percent  of  the  pitch  coefficient.  These  results  show  that  having  an  operational  wind 
turbine  increases  the  overall  damping  of  the  floating  structure  no  matter  which  mode  is  being  excited. 
The  results  of  the  individual  free  decay  tests  are  given  in  Table  8-5. 


Table  8-5  Linear  and  quadratic  damping  coefficients  determined  from  free  decay  tests. 


Test  Description 

Linear 

Damping 

Pi 

Quadratic 

Damping 

P2 

Surge,  wind  turbine  secured 

-0.0019 1 

0.0392 

Surge,  constant  20  m/s  winds,  blade  pitch  =  22° 

0.0052 

0.0486 

Sway,  constant  20  m/s  winds,  blade  pitch  =  22° 

0.0015 

0.0454 

Heave,  constant  20  m/s  winds,  blade  pitch  =  22° 

0.0039 

0.244 

Roll,  constant  20  m/s,  blade  pitch  =  22° 

0.0285 

.0015 

Pitch,  constant  20  m/s,  blade  pitch  =  22° 

0.1665 

.0015  { 

Yaw,  constant  20  m/s,  blade  pitch  =  22°,  yaw  motion  only 

.0715 

.0327 

Surge  free  decay,  7  m/s  winds,  blade  pitch=  2.6°,  Floater 

.0099 

.048 

Surge  free  decay,  7m/s  winds,  blade  pitch=  2.6°,  Nacelle 

.0101 

.043 

Sway  free  decay,  7  m/s,  blade  pitch=  2.6°,  Floater 

.0007 

.0451 

t  negative  value  attributed  to  data  scatter,  f  substituted  roll  value  of p2  for  pitch. 


The  effect  of  blade  pitch  on  system  damping  was  investigated.  The  pitch  angle  of  the  rotor  blades 
was  shown  to  have  a  large  effect  on  the  linear  damping  coefficient  for  the  system.  The  teste  were 
only  performed  for  surge.  But,  the  results  can  be  applied  to  the  other  modes  that  the  wind  turbine 
affected  also.  This  property  could  be  incorporated  into  a  control  algorithm  that  would  control  blade 
pitch  in  order  to  damp  out  oscillations  of  the  system  before  they  become  dangerous. 

The  surge  and  sway  damping  properties  of  the  system  in  the  region  qf  normal  operations  (7  m/s 
winds)  were  evaluated.  It  was  shown  that  the  linear  damping  coefficient  at  this  wind  speed  is  less 
than  10  percent  the  value  of  the  surge  linear  damping  coefficient.  At  20  m/s  the  sway  linear  damping 
coefficient  was  30  percent  of  the  surge  coefficient.  At  the  lower  wind  speed  the  blade  pitch  angle  is 
only  2.6°  instead  of  22°.  This  change  in  blade  pitch  also  resulted  in  the  surge  linear  damping 
coefficient  being  90  percent  higher  for  the  case  with  the  smaller  blade  pitch. 


8.10  Multiple  Mode  System  Damping 

One  of  the  advantages  of  the  analysis  method  used  in  this  research  is  the  ability  to  simulate  complex 
event  and  gather  meaningful  data  about  the  system  response.  How  different  modes  of  motion  interact 
during  an  event  that  involves  more  than  one  excitation  load.  In  an  effort  to  determine  the  system 
behavior  during  this  type  of  event,  free  decay  tests  involving  three  loading  simultaneously  were 
conducted. 

Cross-coupling  occurs  one  of  three  ways.  First  is  hydrodynamic  cross-coupling.  As  the  floater 
moves  in  one  mode,  the  velocity  and  acceleration  of  the  body  creates  counteracting  forces  and/or 
moments  that  excite  another  mode  into  motion.  During  a  free  decay  test  this  behavior  is  characterized 
by  the  driving  mode  moving  at  its  natural  frequency  and  the  cross-coupled  mode  also  oscillates  at  this 
frequency.  The  second  mechanism  is  associated  with  wind  turbines:  gyroscopic  moments.  With  the 
nacelle  aligned  with  the  x-axis,  moving  the  wind  turbine  in  either  pitch  or  yaw  will  excite  a  cross 
coupled  response  from  another  rotational  motion.  The  equation  governing  this  behavior  is  discussed 
in  3.2.1.  As  discussed  in  the  section  on  yaw  damping,  gyroscopic  moments  weakly  couple  two 
rotational  motions. 

The  final  cross-coupling  mechanism  is  mechanical.  Some  rotation  of  the  structure  cause  the  base  of 
the  floater  to  translate  because  the  rotation  point  is  not  located  at  the  base  of  the  floater,  it  is  above  the 
waterline.  This  mechanism  can  be  seen  by  close  observation  of  how  the  heave  mode  damps  out  with 
surge  or  sway  causing  oscillations  at  their  natural  frequency.  Also,  as  the  structure  surges  or  sways, 
the  lines  pull  the  floater  down  into  the  water.  The  high  frequency  heave  oscillations  continue  as  the 
structure  translates  horizontally  and  the  instantaneous  mean  heave  position  changes  with  time. 

There  are  six  possible  modes  to  excite.  As  previously  discussed  pitch/surge  and  roll/sway  are  cross- 
coupled.  Therefore,  any  excitation  of  one  of  the  pairs  will  excite  the  other.  Several  tests  were  run 
and  two  are  discussed  in  this  section.  The  first  used  a  combined  load  which  included  surge,  heave, 
and  yaw.  In  the  second  test  the  system  was  loaded  in  the  sway,  yaw,  and  negative  heave  modes.  The 
loading  curves  are  for  the  two  tests  are  plotted  in  Figure  8.22  and  Figure  8.23. 

At  the  710  seconds  the  excitation  forces  were  secured  and  the  free  decay  test  commenced.  Figure 
8.24  and  Figure  8.25  are  graphs  of  the  system  response.  The  graphs  show  how  the  individual  modes 
damp  out  with  time  and  interact  with  each  other.  A  simulation  was  conducted  with  the  wind  turbine 
secured  in  order  to  further  determine  the  effect  its  operation  has  on  the  system  (see  Figure  8.26). 


Surge  response 


Figure  8.24  Response  of  six  modes  of  motion  to  multiple  loadings  test  one. 
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Figure  8.25  Response  of  six  modes  of  motion  to  multiple  loadings  test  two 


Figure  8.26  Multiple  loading  free  decay  response  with  no  wind  turbine.  Surge,  heave,  and  yaw. 


The  previous  graphs  show  the  response  of  the  system  to  multiple  system  loadings.  Two  simulations 
with  the  wind  turbine  and  one  simulation  with  the  wind  turbine  secured  were  conducted.  The 
following  observations  were  made.  The  cross-coupling  between  the  translational  motions  surge  and 
sway  were  coupled  with  pitch  and  roll  respectively.  This  cross-coupling  was  by  far  the  strongest 
seen. 

Another  cross-coupling  was  that  between  surge  and  sway  motion  with  heave  motion.  As  described 
previously,  the  system  moves  like  an  inverted  pendulum.  When  the  system  surges  or  sways,  it  also 
sinks  down  into  the  water  a  small  amount.  This  affected  the  heave  record  during  its  damping.  A  slow 
oscillation  approximately  40  seconds  in  length  can  be  seen  on  close  observation  of  the  heave  record. 
So  the  overall  heave  motion  appears  to  be  a  high  frequency  oscillation  (heave  natural  frequency) 
superimposed  on  the  surge  /  sway  response. 

The  cross-coupling  between  the  modes  of  motion  appeared  to  have  no  measurable  effect  on  the 
natural  periods  observed.  Cross-coupled  systems  oscillated  with  their  own  natural  period,  the  period 
of  the  system  they  were  coupled  with,  or  both  natural  periods  simultaneously.  No  unexpected  cross¬ 
coupling  mechanisms  were  discovered. 

When  the  wind  turbine  was  secured  all  modes  took  longer  to  attenuate.  The  largest  effect  of  securing 
the  wind  turbine  was  that  the  amplitudes  of  the  pitoh  and  roll  motions  were  much  larger.  The 
oscillations  for  these  rotational  modes  were  dominated  by  the  cross-coupling  to  the  surge  and  sway 
modes  of  motion.  The  wind  turbine  damped  out  the  high  frequency  rotational  motions  associated 
with  the  oscillation  of  at  the  pitch  and  roll  natural  frequencies  very  rapidly.  The  pitch  and  roll 
motions  then  oscillated  at  the  natural  frequencies  of  surge  and  sway.  The  lower  frequency  motions 
resulted  in  smaller  changes  in  wind  inflow  and  hence  lower  daimping  forces  for  the  translational 
motions. 


9  Operational  Evaluation  of  System 


The  floating  wind  turbine  system  must  be  able  to  operate  in  a  random  wind  and  wave  environment  in 
excess  or  twenty  years  with  a  low  probability  of  failure  for  major  components  like  the  blades  and 
tethers.  The  most  probable  cause  of  failure  would  be  fatigue.  One  of  the  goals  of  this  research  was  to 
understand  how  mounting  the  wind  turbine  on  a  floating  base  affects  its  operation.  This  was 
accomplished  by  comparing  the  wind  turbine  loads  for  the  same  wind  turbine  mounted  on  land  and 
the  floating  base. 

The  floating  support  provides  some  flexibility  to  the  system  that  is  not  present  in  conventional  wind 
turbines.  This  flexibility  can  act  to  reduce  the  magnitude  of  excitation  loads  due  to  wind  events.  The 
base  therefore  can  reduce  the  severity  of  a  wind  induced  event.  However,  in  most  instances  where 
there  are  strong  winds  over  an  ocean  environment  there  are  also  waves.  The  waves  force  the  system 
to  move  in  various  modes  of  motion.  This  forced  motion  resulted  in  increase  wind  magnitudes  into 
the  rotor  which  resulted  in  extra  loading  on  the  wind  turbine.  The  forced  oscillations  also  produced 
additional  inertial  loads.  These  loads  were  highest  when  the  system  was  supported  by  flexible  lines 
that  allow  some  pitch  motion.  When  stiffer  lines  were  used  and  pitch  motion  reduced,  wind  turbine 
loading  was  also  reduced. 

9.1  Description  of  Analysis  Model 

The  floater  used  in  the  damping  numerical  simulations  was  also  used  for  the  operational  simulations. 
Polyester  lines  were  used  in  the  damping  tests  to  increase  the  amplitudes  of  oscillation  to  improve 
data  gathering  and  analysis.  The  flexibility  provided  by  those  lines  resulted  in  some  natural  periods 
for  the  system  which  were  not  ideal  for  a  structure  moored  in  an  ocean  environment.  The  lines  were 
changed  to  have  a  maximum  extension  of  4.5%  instead  of  14%.  The  type  of  line  that  would  have  this 
elongation  property  is  also  a  synthetic  line  made  of  High  Modulus  PolyEthylene  (HMPE),  which  also 
has  a  high  strength  to  weight  ratio. 

Changing  the  elongation  characteristic  of  the  tethers  changed  the  restoring  coefficients  for  the  pitch, 
roll  and  heave  modes  of  motion.  Hence,  the  natural  frequencies  for  these  modes  also  changed.  The 
natural  periods  of  all  modes  (except  yaw)  were  either  above  or  below  the  values  for  which  waves 
were  known  to  have  large  spectral  densities.  The  yaw  natural  frequency  was  not  a  concern  since  a 
majority  of  the  submerged  structure  consists  of  vertical  cylinders  which  do  not  generate  yaw 


moments  in  waves.  The  wave  induced  yaw  moments  on  the  spokes  were  very  small  due  to  the  water 
depth  at  which  they  were  attached.  The  wind  turbine  was  the  major  source  of  excitation  for  the 
system  in  yaw.  Changing  the  yaw  natural  period  was  deemed  unnecessary.  The  following  two  tables 
detail  the  parameters  and  calculated  properties  of  the  system  model  used  in  the  operational  evaluation. 

Table  9-1  Dimensions  and  properties  of  floater  used  in  operational  tests. 


Spar  Diameter 

10m 

Line  Diameter 

lQ4mm 

Spar  Length 

12  m 

Number  of  Lines 

1  / spoke 

Tower  Diameter 

8m 

Line  Minimum  Break  Load 

743  toimes 

Tower  Draft 

5  m 

Maximum  Line  Extension 

4.5% 

1 

Spoke  Length 

20  m 

Line  Length  (un-stretched) 

180.76  m 

Spoke  Width 

Im 

Water  Depth 

200  m 

Buoyant  Spokes 

No 

Structure  Thickness 

2.53  cm 

Table  9-2  Calculated  attributes  of  the  floating  wind  turbine  used  in  the  operational  tests 


Mass  of  Floater 

210435.4  kg 

Center  of  Gravity  (system)  CG 

16.98  m 

Buoyant  Mass 

1228424.6  kg 

Center  of  Buoyancy  CB 

-10.89  m 

Total  Mass  of  System 

411494.5  kg 

CG  -  CB 

27.87  m 

Reserve  Buoyancy 

198.5  % 

Surge  Natural  Period 

39.02  sec 

Initial  Line  Stress 

235.4  MPa 

Heave  Natural  Period 

2.74  sec 

Breaking  Line  Stress 

857.7  MPa 

Pitch  /  Roll  Natural  Period 

5.88  sec 

Center  of  Gravity  (floater) 

Yaw  Natural  Period 

10.51  sec 

9.2  Test  Descriptions 

Wind  turbine  designs  are  proven  to  be  satisfactory  by  evaluating  their  ultimate  and  fatigue  strength 
characteristics.  The  magnitude  and  periodicity  of  loading  events  govern  the  fatigue  life  of  the  system. 
For  some  system  components  fatigue  failure  is  the  dominant  mechanism.  Other  components  are  more 
susceptible  to  failure  due  to  exceeding  the  ultimate  yield  strength  of  the  component.  Both  of  these 
strength  criteria  can  be  evaluated  using  the  appropriate  aero-elastic  time  domain  model.  The 
governing  document  for  wind  turbine  design  and  safety  is  lEC  61400-1  [19]. 

Table  2  in  the  Structural  Design  section  of  lEC  61400-1  listed  the  design  load  cases.  In  order  for  a 
wind  turbine  design  to  become  type  certified  it  must  have  passed  evaluations  for  all  of  the  applicable 
load  cases.  In  the  course  of  the  operational  evaluation  of  the  floating  wind  turbine  system  the  design 
load  cases  were  used  as  guidance  for  selecting  test  cases.  The  goal  of  this  research  was  not  to  achieve 
type  certification;  therefore  not  all  of  the  cases  were  addressed.  Actual  certification  of  the  floating 
design  is  left  for  future  work.  Evaluations  of  the  floating  platform  were  obtained  by  evaluating  a 
specific  load  case  with  the  wind  turbine  mounted  on  the  floater  and  a  wind  turbine  with  a  fixed  base. 
In  accordance  with  section  12.1  of  Guideline  DGOl  [29]  the  following  loads  were  evaluated  during 
the  numerical  simulations: 

•  Blade-root  flapwise  bending  moment 

•  Blade-root  edgewise  bending  moment 

•  Tower-head  torsion 

•  Low-speed  shaft  bending  moments  (2  directions) 

•  Low-speed  shaft  torque  (rotor  torque) 

•  Tower-base  bending  moments  (2  directions) 

The  first  operational  tests  which  were  conducted  involved  simulating  the  operation  of  the  two 
systems  (floater  and  fixed)  at  three  different  wind  speeds.  For  these  tests  the  wind  turbine  was 
operating  in  speed  regulate  (blades  pitch  to  maintain  rotor  speed).  The  wind  profiles  used  in  the 
simulations  were  generated  with  SNWind  (NWTC  design  code).  Wind  profiles  were  generated  using 
the  power  law  equation  (9.1).  For  offshore  wind  turbines  Germanischer  Lloyd  (GL)  recommended 
using  a  power  law  exponent  a  =  0.11  (GL-Regulation  [14]). 


'  z  ' 


\'^r  J 


(9.1) 


where  V{z)  is  the  wind  velocity  at  height  z.  V{z^  is  the  wind  speed  at  the  reference  height  z,. 


The  ultimate  strength  test  which  was  chosen  for  evaluation  was  the  Extreme  Operating  Gust  (EOG). 
This  wind  load  response  was  evaluated  using  the  guidance  of  lEC  61400-1  [19].  For  the  purposes  of 
testing  the  1.5  MW  Baseline  wind  turbine  was  assumed  to  be  a  lEC  Class  II-A  wind  turbine.  Table  1 
of  section  6  in  lEC  61400-1  gives  the  basic  parameters  that  apply  for  the  difference  classes  of  wind 
turbines.  This  class  of  wind  turbine  has  a  reference  velocity  of  42.5  m/s  and  a  turbulence 
intensity  (4^/)  of  0.16.  The  hub  height  gust  has  a  periodicity  of  once  per  year  and  a  gust  magnitude 
given  by: 


l-35(n,-F„);  3.3 


(9.2) 


where 


Standard  deviation  of  turbulence; 

One  year  extreme  wind  speed: 

Longitudinal  turbulence  scale  parameter: 


cT,=/,^(0.75F,„,+b);  b  =  5.6m/s 


F.,=1.12 


ref 


N0.n 
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K^hub  J 
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The  test  wind  record  is  found  using  the  gust  velocity  to  modify  the  wind  velocity  as  follows: 

Jf(z) -  0.37 sin(3;rt /r)(l- cos(2;rt / F))  /or  0  <  t  ^  r 

r  j  —  % 

[  F(z)  otherwise 


(9.3) 


The  extreme  operating  gust  was  generated  for  a  hub  height  wind  velocity  of  7  m/s.  This  gust  was 
used  as  the  wind  input  for  a  10  minute  simulation. 

The  systems  response  to  extreme  wave  evente  was  tested.  Extreme  waves  are  generally  produced  by 
extreme  winds.  For  this  case  the  wind  speed  was  above  the  rated  cutout  speed  for  the  1.5  MW 
baseline  wind  turbine.  The  wind  turbine  model  was  generated  for  a  shutdown  condition.  The  wind 
turbine  was  secured  and  the  rotor  blades  were  set  at  90^  to  minimize  wind  loads.  This  mitigated  the 


wind  loads  on  the  structure.  The  system  was  subjected  to  an  extreme  wave  group  with  two  different 
wave  incident  angles. 

The  amplitudes  of  extreme  waves  can  be  using  a  probabilistic  approach  which  is  described  in 
Sclavounos  [31].  Assumptions  are  made  in  order  to  determine  the  amplitude  for  the  extreme  wave 
simulations.  The  amplitude  of  the  extreme  wave  is  called  the  design  peak  (y ).  For  these  analyses, 
the  design  peak  is  assumed  to  be  the  amplitude  of  largest  wave  passing  the  floating  structure  during  a 
three  hour  time  period.  For  a  given  spectrum  described  by  a  significant  wave  height  (////?)  and 
average  period  {To)  the  design  peak  for  the  time  interval  (7)  is  found  using  the  following  equations: 


Design  Peak:  y  = 


Number  of  Peaks;  n  =  - 


Reliability:  a  s  Pr(y  >  y) 


H 


Standard  Deviation:  cr  s  ■ 


X 


(9.4) 


Where  a  is  the  reliability  parameter  which  is  chosen  to  be  a  small  value  to  indicate  the  low  likelihood 
of  the  design  peak  being  exceeded  during  the  given  time  interval.  The  value  used  in  these  analysis 
was  a  =  0.01.  The  probability  of  the  design  peak  being  exceeded  in  the  three  hour  period  is  therefore 
1  percent. 


The  value  for  the  design  peak  was  determined  using  the  wave  spectmm  which  corresponded  to  the 
reference  wind  speed  which  resulted  in  the  hub  height  wind  speed  exceeding  the  cutout  speed  for  the 
machine.  The  two  wave  incidence  angles  were  assumed  to  be  0°  and  30°.  These  simulations  were 
ten  minutes  in  length.  During  the  test  the  wave  elevation  was  ramped  up  to  the  design  peak,  held  for 
one  cycle,  and  then  ramped  down. 


9.3  Normal  Power  Generation  Cases 

The  purpose  of  the  operational  simulation  tests  were  to  determine  the  effects  mounting  a  wind  turbine 
on  a  tension  leg  floating  platform  has  on  the  wind  turbine.  Also  the  effects  on  the  mooring  lines  had 
to  be  evaluated.  Numerical  simulations  were  conducted  with  different  wind  and  wave  conditions. 
The  loads  at  various  points  in  the  wind  turbine  of  the  floating  system  were  compared  to  the  loads  on 
an  identical  wind  turbine  with  a  fixed  base.  The  wind  environments  are  the  same  for  both  the  floating 
and  fixed  wind  turbine. 

Three  reference  wind  speeds  were  chosen:  7  m/s,  11  m/s,  and  15  m/s  at  a  reference  height  of  10 
meters.  The  ocean  waves  for  the  tests  were  assumed  to  be  generated  by  the  winds.  The  wave 
characteristics  were  chosen  to  match  the  reference  wind  speed.  The  wave  spectrum  was  model  as  a 
Pierson-Moskowitz  type  with  a  significant  wave  height  (Hia)  and  a  mean  wave  period  (T|).  The 
values  for  H1/3  and  Ti  were  extrapolated  from  Table  6.1  of  Newman,  1997  [28]. 


Table  9-3  Environmental  conditions  for  power  generation  cases.  Reference  height  for  wind  -10  m. 

Wave  spectrum  used  -  Pierson-Moskowitz. 


Test  # 

Wind  Inputs 

Wave  Inpute 

Reference 
Wind  Speed 
(m/s) 

Mean  Wind 
Speed  (hub) 
(m/s) 

Sea  State 

mg 

Mean  Wave 
Period  -  T| 
(seconds) 

1 

7 

8.67 

3 

1.2 

4.344 

2 

11 

13.63 

4 

2.4 

6.19 

3 

15 

18.59 

6 

5.0 

8.688 

For  the  first  three  cases  both  the  wind  and  waves  were  directed  along  the  positive  x-axis.  Each  of  the 
first  three  simulations  had  a  length  of  one  hour.  The  random  seeds  used  to  generate  the  wind  and 
wave  records  were  changed  for  each  case.  The  data  collected  during  each  run  was  processed  using  a 
program  called  Crunch  (NWTC  design  code).  This  program  was  used  to  obtain  statistical  and  fatigue 
data  fi’om  the  time  records  of  wind  turbine  and  line  loading.  The  statistical  data  consisted  of  [3]: 

•  Minimum 

•  Mean 

•  Maximum 

•  Maximum  range 

•  Standard  deviation 


# 


106 


During  normal  operations  wind  and  waves  loads  cause  the  floater  to  move.  Some  modes  of  motion 
are  more  affected  than  others.  Figure  9.1  is  an  example  of  the  floater  response  in  a  random  wind  and 
wave  environment.  The  solid  red  lines  represent  the  response  with  both  wind  and  waves.  The  blue 
dotted  lines  represent  the  response  of  the  floater  with  only  wind  acting  on  the  system.  By  comparing 
the  two  sets  of  responses,  it  can  be  seen  that  waves  have  a  great  effect  on  heave  response.  The  thrust 
produced  by  wind  on  the  rotor  results  in  surge  and  pitch  motion.  As  previously  discussed  surge 
motion  is  coupled  with  pitch  motion.  The  surge  motion  is  mechanically  cross-coupled  with  heave 
motion.  Horizontal  motion  causes  set-down.  Y aw  motion  is  primarily  caused  by  wind  acting  on  the 
rotor  blades. 


9.3.1  Wind  Turbine  Results 


The  results  presented  in  this  section  seek  to  contrast  the  differences  between  the  loading  on  a  fixed 
base  wind  turbine  and  a  floating  wind  turbine.  Failure  of  a  component  can  occur  when  maximum 
yield  strength  is  exceeded  or  a  fatigue  limit  is  passed.  First  a  comparison  of  the  maximum  loading  for 
the  floating  and  non-floating  variants  is  made.  The  eight  loadings  listed  in  the  previous  section  were 
collected  for  the  hour  long  simulations  and  their  maximum  (positive  or  negative)  values  charted  for 
the  three  different  wind  and  wave  conditions  (see  Figure  9.2  to  Figure  9.4).  The  following 
abbreviations  were  used  on  the  charts  and  the  following  discussions  for  the  loads: 

TwrBsMxt:  tower  base  moment  (x-axis) 

TwrBsMvt:  tower  base  moment  (y-axis) 

RootMxb2:  blade-root  edgewise  bending 

moment  (blade  2) 

RootMvb2:  blade-root  flapwise  bending 
moment  (blade  2) 

The  charts  show  which  loadings  were  most  affected  by  mounting  the  wind  turbine  on  a  floating  base 
and  subjecting  it  to  excitation  by  waves.  In  general  the  floating  wind  turbine  had  higher  maximum 
values  for  the  eight  loading  values.  However,  two  of  the  parameters  were  higher  for  the  fixed  base 
wind  turbine.  These  two  loads  (RootMxb2  and  LSShtMxa)  were  higher  for  the  fixed  base  wind 
turbine  in  the  first  case  with  reference  wind  equal  to  7  m/s.  All  of  the  maximums  were  greater  for  the 
floating  system  in  the  tests  two  and  three. 

The  bending  moment  at  the  tower  base  showed  the  largest  difference  between  the  two  wind  turbines. 
The  system  was  allowed  some  slight  movement  in  the  roll  and  pitch  modes  when  the  wind  tuibine 
was  mounted  on  the  floating  base.  With  the  wind  and  waves  directed  along  the  x-axis,  large  pitching 
moments  developed  and  pitch  deflection  occurred.  The  pitching  motion  of  the  structure  tilted  the 
tower  and  generated  large  additional  moment  at  the  tower  base. 

A  summary  of  the  three  tests  is  presented  in  Table  9-4.  The  percentage  difference  between  the  floater 
and  fixed  base  values  for  the  eight  loads  are  listed.  These  values  are  the  percentage  of  the  fixed  base 
loading.  Negative  values  reprraent  the  cases  where  the  fixed  base  maxima  were  greater.  Both  the 
low-speed  shaft  tip  bending  and  blade  root  bending  moments  showed  large  increases  their  maxima. 


YawBrMzn:  tower-head  toraion 

LSSTipMvs:  low-speed  shaft  bending  moment 
(y-axis) 

LSSTipMzs:  low-speed  shaft  bending  moment 
(z-axis) 

LSShftMxa:  low-speed  shaft  torque 
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LOADING 

Test  1 

Test  2 

Test  3 

Tower  Base  Bending  (x-axIs) 

140.7% 

164.9% 

145.1% 

Tower  Base  Bending  (y-axis) 

27.6% 

77.1% 

124.0% 

Biade-Root  Edgewise  Bending 

-2.6% 

10.4% 

44.9% 

Blade-Root  Fiapwise  Bending 

0.3% 

15.4% 

34.8% 

Low-Speed  Shaft  Tip  Bending  (y-axis) 

-0.2% 

53.3% 

17.0% 

Low-Speed  Shaft  Tip  Bending  (z-axis) 

29.0% 

43.0% 

37.5% 

Low-Speed  Shaft  Torque 

-0.3% 

6.6% 

3.3% 

Tower-Head  Torsion 

26.9% 

40.2% 

38.6% 

110 


The  largest  percentage  difference  was  between  the  fixed  and  floating  values  for  tower  base  moment 
about  the  x-axis.  The  actual  values  for  the  maxima  for  this  load  are  low  relative  to  the  other  tower 
base  load.  The  percentage  differences  for  the  blade-root  bending  moments  and  the  tower  base 
bending  moment  about  the  y-axis  increased  as  the  environmental  conditions  went  up. 

The  effect  on  wind  turbine  fatigue  due  to  the  floating  support  can  be  evaluated  by  determining  the 
rainflow  cycle  counts  or  raw  cycles  for  the  oscillating  loads.  These  rainflow  cycle  counts  are 
obtained  by  analyzing  the  time  history  for  the  six  simulations  using  the  Crunch  post-processor.  The 
cyclic  loads  are  separated  into  bins.  Each  bin  holds  the  number  of  oscillations  which  have  a  given 
range.  The  width  of  the  bin  defines  the  ranges  of  cyclic  loads  which  are  counted.  The  cycle  counts 
can  be  used  to  determine  the  effect  on  the  fatigue  life  of  the  wind  turbine  components.  Therefore,  the 
cycle  counts  can  be  used  to  compare  the  fatigue  of  the  fixed  base  and  floating  wind  turbines. 

The  results  of  the  rainflow  cycle  counts  for  the  six  simulations  are  shown  in  the  following  charts. 
Each  graph  shows  the  twenty  bins  and  the  number  of  cycles  in  each  for  the  fixed  and  floating  base 
wind  turbines.  The  floating  support  has  more  extreme  oscillations  than  the  fixed  base  support.  From 
a  design  standpoint  the  increased  cyclic  loading  on  the  blade-roots  and  low-speed  shaft  are  a  cause  for 
concern.  A  wind  turbine  which  could  be  used  for  the  fixed  base  machine  would  probably  be 
unacceptable  for  the  floating  system.  Once  again  the  largest  differences  in  loading  are  in  the  tower 
base  bending  moments  which  could  be  compensated  for  by  extra  structural  support  in  the  lower 
portions  of  the  tower. 
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Figure  9.5  Rainflow  cycle  counts  for  test  1  fixed  base  and  floating  wind  turbine.  Reference  wind 

speed  =  7  m/s.  Sea  state  3. 
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Figure  9.6  Rainflow  cycle  counts  for  test  2  fixed  base  and  floating  wind  turbine.  Reference  wind 

speed  =  1 1  m/s.  Sea  state  4. 
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9.3.2  Comparisons  of  Wind  Turbine  Power  Production 

Understanding  how  mounting  a  wind  turbine  on  a  floating  support  structure  affects  the  machines 
power  output  is  important.  If  changing  the  way  it  was  mounted  dramatically  altered  the  power 
output,  the  system  would  be  unfeasible.  This  was  found  not  to  be  the  case  with  the  system  tested  in 
the  three  sets  of  operating  conditions. 

The  power  output  for  the  fixed  base  and  floating  wind  turbines  was  compared  by  determining  the 
average  power  output  for  each  of  the  three  operational  tests.  It  was  unclear  what  the  cause  of  this 
reduction  in  power  was.  Several  parameters  for  each  test  were  reviewed.  The  most  likely  cause  is  the 
flexibility  of  the  system  allowed  some  of  the  power  that  would  normally  be  transferred  to  the 
generator  to  be  passed  on  to  the  floating  support.  The  motion  of  the  floater  absorbs  a  small  portion  of 
the  energy  that  would  have  been  transferred  to  the  generator.  Observation  of  the  sway  and  roll  modes 
of  motions  showed  motion  which  had  to  have  been  excited  by  the  wind  turbine  since  the  waves  were 
directed  along  the  x-axis.  This  motion  was  damped  by  viscous  drag  on  the  floater.  The  damping  due 
to  drag  on  the  structure  is  an  energy  loss.  Another  reason  for  the  power  loss  could  be  the  increased 
pitch  and  yaw  motion  of  the  rotor  due  to  the  motion  of  the  floater.  As  shown  in  Table  9-5,  the 
maximum  reduction  in  power  was  only  1 . 1  percent. 

Table  9-5  Average  power  output  for  fixed  base  and  floating  wind  turbine  for  the  three  operating 

conditions. 


Test# 

Fixed  Base  Average 
Power  Output 
(kW) 

Floater  Average 
Power  Output 
(kW) 

Change  in 
Power  Output 
(%) 

Test1 

740.3 

733.4 

-0.9% 

Test  2 

1373.8 

1365.1 

-0.6% 

Test  3 

1380.0 

1364.8 

-1.1% 
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9.3.3  Tether  Results 

A  majority  of  the  restoring  forces  and  moments  for  the  floating  system  came  from  the  tethers.  As 
previously  discussed  pitch  and  roll  moments  caused  the  largest  changes  in  line  tension.  The  first 
three  tests  cases  used  a  single  104  mm  synthetic  line  on  each  spoke.  The  reserve  buoyancy  of  the 
system  resulted  in  a  static  line  tension  of  approximately  30  percent  of  the  ultimate  stress  for  the 
chosen  lines.  When  the  wind  and  waves  aligned  with  the  x-axis  the  highest  stresses  were  experienced 
on  line  number  3  (on  negative  x-axis).  Lines  1  and  3  absorbed  the  pitching  moments  acting  on  the 
system.  , 


Figure  9.8  is  a  graph  showing  a  100  second  record  of  line  stress  for  the  three  sets  of  environmental 
conditions.  This  graph  shows  that  the  oscillations  in  the  line  stress  increase  in  severity  as  the  wind 
speed  and  sea  state  increase.  The  graphs  for  line  number  1  and  3  look  similar  except  that  they  are 
opposite  in  sign.  Line  3  stress  increases  and  line  1  stress  decreases  for  positive  pitching  momente. 


Line  stress  response  for  three  operational  conditions. 
Line  number  3. 


The  time  records  of  the  line  stress  were  analyzed  in  order  to  evaluate  the  effect  of  the  three  test 
conditions  on  the  tethers.  The  statistical  results  were  recorded  in  Table  9-6.  As  expected  the 
maximum  line  stresses  increased  as  the  environmental  conditions  worsened.  The  maximum  line 
stress  for  the  three  operational  cases  was  433,5  MPa,  which  was  approximately  50  %  of  the  breaking 
stress. 
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Table  9-6  Line  stress  statistics  for  operational  tests  1,  2,  and  3. 


Test  1 

Test  2 

Test  3 

Minimum 

Stress 

(MPA) 

Mean 

Stress 

(MPA) 

Minimum 

Stress 

(MPA) 

Minimum 

Stress 

(MPA) 

Mean 

Stress( 

MPA) 

Maximum 

Stress 

(MPA) 

Minimum 

Stress 

(MPA) 

Mean 

Stress 

(MPA) 

Maximum 

Stress 

(MPA) 

Uriel 

139.2 

204.9 

270.2 

79.2 

204.1 

328.1 

30.8 

214.0 

386.3 

Line2 

223.7 

237.3 

250.7 

215.1 

238.4 

259.7 

200.6 

239.5 

279.9 

Line3 

200.0 

267.3 

332.4 

132.0 

268.1 

392.2 

79.7 

258.1 

433.5 

Line4 

222.1 

234.9 

247.3 

208.1 

233.8 

254.3 

192.7 

232.6 

272.5 

The  tethers  were  made  of  a  high  modulus  polyethylene  material  which  is  reported  to  have  better 
fatigue  properties  that  steel  wire  ropes  [8].  Some  test  data  was  found  that  backed  up  this  claim.  But 
no  comprehensive  tension  -  tension  fatigue  life  graphs  were  found  for  the  high  modulus  polyethylene 
It  was  decided  that  a  conservative  estimation  of  the  fatigue  loading  of  the  tethers  could  be  obtained  by 
using  the  fatigue  graphs  associated  with  steel  wire  ropes.  Using  the  dimensions  and  fatigue  properties 
allowed  the  equivalent  damage  due  to  fatigue  loading  to  be  calculated  for  the  three  one  hour 
simulations.  Damage  was  assumed  to  accumulate  linearly  for  each  oscillating  load  using  an  equation 
take  from  Annex  G  of  61400-1  [19]: 


Damage  = 


(9.5) 


where  5/  is  the  load  range  for  cycle  i,  and  N(Si)  is  the  number  of  cycles  to  failure  for  the  load  range  St. 
The  values  for  S,  are  obtained  from  the  load  range  values  of  the  rainflow  cycle  bins;  The  number  of 
cycles  to  failure  is  derived  from  the  fatigue  life  cycle  graph.  The  equivalent  damage  for  each  of  the 
three  tests  is  found  by  combining  the  damage  from  every  cycle  in  the  bins.  A  value  of  damage  of  1.0 
or  higher  means  failure  is  predicted.  The  damage  values  for  lines  1  and  3  during  the  second  and  third 
tests  are  the  highest  (see  Table  9-7).  The  life  of  the  chosen  lines  could  be  extremely  low  (<  25  hours 
for  lines  1  and  3  in  test  3)  in  high  wind  and  sea  conditions  and  would  therefore  be  unacceptable. 


Table  9-7  Total  damage  due  to  fatigue  for  lines  1-4  during  tests  1 ,  2,  and  3. 


Test# 

Line  1 

Line  2 

Line  3 

Line  4 

Test  1 

1.94E-05 

1.77E-09 

2.15E-05 

1.45E-09 

Test  2 

2.81  E-03 

3.27E-08 

3.84E-03 

2.65E-08 

Test  3 

4.04E-02 

1.34E-06 

4.26E-02 

8.85E-07 

9.4  Design  Modification 

The  results  of  the  firet  three  operational  tests  showed  that  the  floating  design  resulted  in  increased 
loading  on  the  wind  turbine.  The  fatigue  life  of  the  wind  turbine  was  also  shown  to  be  adversely 
affected.  These  problems  could  be  mitigated  by  design  changes.  Although  design  optimization  was 
not  a  goal  for  this  research,  it  was  decided  that  a  single  change  would  be  made  to  the  system  to 
demonstrate  how  the  loading  properties  of  the  system  could  be  adjusted  to  overcome  problems. 

The  first  step  in  this  design  iteration  was  to  determine  the  most  likely  cause  for  a  majority  of  the 
increased  load  on  the  wind  turbine  components.  It  was  found  that  the  most  of  the  increase  in  load  on 
the  tower  base  was  due  to  the  structure  pitching,  by  comparing  the  pitch  angle  response  of  the  floater 
to  the  tower  base  bending  moment  (y-axis).  A  positive  pitch  angle  resulted  in  a  positive  bending 
moment  (see  Figure  9.9).  Other  wind  turbine  loads  were  also  found  to  be  largely  affected  by  either 
pitch  angle  or  rate  of  pitch. 


Figure  9.9  Graph  of  pitch  angle  and  tower  base  bending  moment  (y-axis)  for  operational  test  3. 


The  next  step  was  to  decide  which  design  parameter  to  change  in  order  to  reduce  the  pitch  motion  of 
the  floater.  The  three  parameters  which  had  the  largest  effect  on  pitch  motion  were  the  type  of  lines 
selected,  number  of  lines,  and  spoke  length.  Either  of  these  parameters  could  be  changed  to  achieve 
the  desired  result  of  reduced  wind  turbine  loading.  The  number  of  lines  was  increased  to  three  per 
spoke  and  lines  with  a  higher  stiffness  were  used.  The  maximum  elongation  was  changed  from  4.5% 


to  2.0%,  which  is  consistent  with  another  type  of  high  modulus  polyethylene  line  [9].  This  change 
resulted  in  a  system  which  was  much  stiffer  in  the  pitch  mode  of  motion  with  natural  periods  of  pitch 
of  2.3  seconds.  Operational  test  3  was  repeated  for  the  new  design.  The  surge,  pitch,  and  heave 
responses  for  both  designs  were  plotted  in  Figure  9.10.  Three  of  the  wind  turbine  loads  were  also 
plotted  for  both  designs.  The  largest  load  reduction  was  the  y-axis  tower  bending  moment.  Limiting 
the  tower  angle  reduced  the  gravity  induced  moment  at  the  base  of  the  tower.  The  reduced  pitch 
motion  improved  the  systems  behavior  for  all  eight  wind  turbine  loads  (see  Table  9-8). 


Figure  9.10  Surge,  heave,  and  pitch  response  and  wind  turbine  loads  for  1  and  3  lines  per  spoke  for 

operational  test  3  conditions. 


Table  9-8  Comparison  of  fixed  base  and  floating  wind  turbines  for  1  and  3  lines  per  spoke.  The 
reference  wind  speed  was  15  m/s.  Sea  state  6.  The  values  represent  the  percentage  difference 
between  the  two  floating  variants  and  the  fixed  base  wind  turbine  maximum  loading  values. 


LOADING 

Fixed  Base 
Loading 
(kNm) 

1  1  Line  /  Spoke 

3  Lines  /  Spoke  | 

Loading 

(kNm) 

Delta 

(%) 

Delta 

(%) 

Tower  Base  Bending 
(x-axis) 

5956.0 

145.1% 

6604  ■ 

10.9% 

Tower  Base  Bending 
(y-axis) 

35060.0 

78550.0 

124.0% 

45390 

29.5% 

Blade-Root  Edgewise 
Bending 

700.5 

1015.0 

44.9% 

741.9 

5.9% 

Blade-Root  Flapwise 
Bending 

1797.0 

2422.0 

34.8% 

1866 

3.8% 

iTfnfifffifiiil 

1213.0 

1419.0 

17.0% 

1325 

9.2% 

486.6 

669.1 

37.5% 

467.9 

-3.8% 

884.0 

913.4 

3.3% 

843 

-4.6% 

741.2 

1027.0 

752.5 

As  shown  in  Figure  9.1 1,  the  fatigue  loading  of  the  wind  turbine  was  greatly  reduced.  Increasing  the 
stiffness  in  pitch  resulted  in  smaller  pitch  motions  and  angular  velocities.  The  total  effect  of  these 
changes  reduced  the  range  of  cyclic  load  for  the  wind  turbine.  The  results  of  the  numerical 
simulation  indicated  that  some  of  the  loadings  had  lower  cyclic  loads  than  the  fixed  base  wind  turbine 
even  with  sea  state  6  seas  acting  on  the  floater.  The  flexibility  of  the  system  in  surge  was  not 
significantly  affected  by  the  change  of  lines.  The  system'  surged  when  wind  loading  events  occurred 
and  thereby  relieved  some  of  the  loads  on  the  wind  turbine.  The  reduction  in  wave  induced  pitch 
motion  (since  system  was  more  stiff  in  the  pitch  mode)  dramatically  reduced  the  loads  on  the  wind 
turbine. 

The  increased  number  of  lines  improved  the  fatigue  results  for  the  tethers.  The  maximum  value  for 
total  damage  was  1 .9  E-7  instead  of  4.3E-2.  This  change  would  greatly  increase  the  fatigue  life  of  the 
lines.  The  maximum  line  stress  for  the  three  line  design  was  126.6  MPa  while  the  single  line  design 
had  a  maximum  stress  of  433.5  MPa.  The  results  of  this  test  showed  that  design  changes  could  be 
used  to  overcome  adverse  loading  conditions  on  the  wind  turbine  and  lines.  The  improved  design 
was  used  in  all  subsequent  tests. 


9.4.1  Design  Verification 

While  the  results  of  a  single  simulation  may  be  enough  to  determine  the  unsuitability  of  a  design,  one 
simulation  is  not  adequate  to  prove  the  acceptability  of  a  design.  The  environmental  inputs  can  vaiy 
greatly  between  different  simulations.  The  random  number  generators  used  to  produce  the  wind  and 
wave  records  determine  the  environmental  inputs  for  the  floating  wind  turbine  system.  A  single  test 
is  insufficient  and  several  simulations  must  be  conducted  before  a  design  can  truly  be  proven 
adequate. 

The  response  of  the  modified  floating  wind  with  three  tethers  per  spoke  was  further  tested  using 
operational  test  3  conditions  (sea  state  6  and  1 5  m/s  reference  winds).  Five  floating  and  five  fixed 
base  simulations  were  conducted  using  different  random  seeds  to  generate  the  input  wind  records  and 
wave  spectrums.  Each  simulation  was  ten  minutes  in  length.  Table  9-9  is  a  summary  of  the  ten 
simulations.  The  maximum  loading  values  for  all  but  the  tower  base  bending  moments  vvere  higher 
for  the  fixed  base  than  the  floating  support  results.  This  general  behavior  is  not  shown  in  the  single 
one  simulation.  During  the  one  hour  simulation  wave  loading  was  the  driving  factor  behind  the 
maximum  loading  on  the  floating  wind  turbine  since  most  of  floating  loadings  were  highest  for  the 
floating  support.  The  value  of  ruiming  multiple  tests  was  therefore  demonstrated  since  the  results  of 
the  single  test  showed  maximum  loading  values  being  higher  for  the  floater  being  higher  than  the 
fixed  base  for  four  additional  loads. 


Table  9-9  Comparison  of  mean  and  maximum  loading  values  for  5  fixed  base  and  5  floating  support 

10  minute  simulations  using  operational  test  3  conditions. 


Loading 

Mean  Loading 

Maximum  Loading 

Delta 

(%) 

Tower  Base  Bending 
(x-axis) 

1223.4 

1260.5 

3.0% 

5790.0 

6075.0 

4.9% 

730219 

7393.7 

1.2% 

41720.0 

53180.0 

27.5% 

Blade-Root  Edgewise 

Bending 

36.9 

38.0 

2.8% 

1254.0 

1009.0 

-19.5% 

Blade-Root  Flapwlse 
Bending 

503.7 

504.1 

0.1% 

2060.0 

2022.0 

-1.8% 

253.5 

250.8 

-1.0% 

1318.0 

Bi 

-6.0% 

Low-Speed  Shaft  Tip 

Bending  (z-axis) 

89.9 

92.4 

2.8% 

846.2 

-27.4% 

Low-Speed  Shaft 

Torque 

677.6 

■^1 

-0.3% 

1010.0 

-7.7% 

Tower-Head  Torsion 

128.1 

130.8 

2.1% 

862.0 

636.3 

-26.2% 

Running  multiple  tests  using  different  random  seeds  for  the  environmental  inputs  showed  how  a  stiff 
support  system  on  the  floater  can  result  in  lower  maximum  operating  loads  than  a  fixed  base  wind 
turbine.  This  effect  of  wave  loading  on  the  maximum  loading  values  for  all  but  the  tower  base 
bending  moments  was  mitigated  by  the  fact  that  the  stiffer  tethers  did  not  allow  the  waves  to  induce 
excessive  motion  of  the  nacelle  due  to  pitch  /  roll.  The  flexibility  of  the  floating  system  in  surge  and 
sway  relieved  some  of  the  wind  stresses  placed  on  the  system. 

The  fatigue  response  of  the  system  was  also  determined  for  the  ten  simulations.  The  rainflow  cycle 
counts  were  determine  for  each  of  the  tests.  Then  the  results  of  the  five  floating  and  five  fixed  base 
test  were  combined  into  two  rainflow  cycle  records  for  the  two  wind  turbine  base  conditions  (see). 
Higher  numbers  of  cycles  at  the  higher  loading  values  will  result  in  decreased  component  life.  As 
shown  in  Figure  9.12,  both  of  the  tower  base  bending  moments  will  have  higher  fatigue  for  die 
floating  wind  turbine.  The  loading  values  for  the  other  six  components  were  shown  to  be  very  similar 
for  both  the  fixed  and  floating  cases.  The  additional  loading  on  the  system  had  very  little  effect  on 
the  cyclic  loading  on  the  parts  at  the  tower  top. 


Cy«l«$  /  Hour  |  Cyclot  /  Hour  I  I  Cycles  /  Hour 


Figure  9.12  Combined  rainflow  cycle  counts  for  five  floating  (3  lines  per  spoke)  and  five  fixed  base 
simulations.  Simulation  times  =  1 0  minutes.  Sea  state  6  and  1 5  m/s  reference  winds. 
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9.5  Extreme  Wind  Loading  Case 

The  response  of  the  system  to  an  abnormal  wind  event  can  be  evaluated  using  the  Extreme  Operating 
Gust  (EOG)  described  in  section  9.2.  During  the  EOG  event  the  wind  inflow  to  the  rotor  initially  dips 
down,  then  rises  sharply  to  the  peak  wind  velocity,  followed  by  another  dip,  and  finally  comes  back 
to  the  original  wind  speed.  This  event  occurs  over  a  10  second  period.  The  hub  height  wind  record 
used  in  this  testing  is  shown  in 


Figure  9.13.  This  wind  record  was  generated  using  the  lECWind  program  from  the  National  Wind 
Technology  Center. 


Figure  9.13  Wind  velocity  record  during  Extreme  Operating  Gust  for  a  class  II-A  wind  turbine  with  a 

rated  speed  of  1 0  m/s  and  cutout  speed  of  25  m/s. 

Numerical  simulations  were  conducted  on  the  system  for  four  different  cases.  First  the  fixed  based 
wind  turbine  was  evaluated.  Next,  a  floating  wind  turbine  with  no  incident  waves  was  subjected  to 
the  EOG.  Finally,  two  tests  with  sea  state  3  waves  and  the  EOG  were  conducted.  The  second  wave 
test  used  the  same  wave  spectra  as  the  first  sea  state  3  test  with  the  phase  angles  for  the  spectrum 
bands  shifted  by  180  degrees.  This  phase  shift  resulted  in  the  wave  excitation  forces  on  the  floater 
being  likewise  shifted  by  180  degrees. 

As  shown  in  Figure  9.14,  mounting  the  wind  turbine  on  a  floating  base  changes  the  maximum  values 
of  loading  during  the  extreme  operating  gust  event.  Most  of  the  values  for  loading  were  higher  for 
the  floater  that  the  fixed  base  wind  turbine.  But  in  most  cases  the  differences  were  small.  The  waves 
acting  on  the  floating  base  were  seen  to  affect  the  values  for  maximum  loading. 
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Figure  9.14  Load  comparisons  for  extreme  operating  gust.  Four  cases;  fixed  base,  floating  with  no 

waves,  floating  with  sea  state  3-1"  case,  and  floating  with  sea  state  3-2"**  case  (phase  shifted  180°). 

In  general,  loading  values  for  the  floating  wind  turbines  were  higher  than  the  fixed  base  wind  tarbine. 
However,  some  of  the  loading  values  were  lower  for  the  floating  wind  turbine.  All  of  the  blade-root 
moments  were  slightly  lower  for  the  floating  turbine  than  the  fixed  turbine.  The  largest  differences 
were  between  the  values  of  low-speed  shaft  tip  bending  moments  (see  Table  9-10).  But,  the  values 
are  low  compared  to  the  moments  of  the  shaft  tip  about  the  y-axis  and  therefore  are  not  very 
troubling. 


Waves  acting  on  the  floater  produced  forces  and  moments  on  the  system.  In  the  firet  test  the  wave 
loads  acted  in  the  opposite  direction  as  the  wind  loads  during  the  10  second  period  of  the  EOG. 
Wave  loads  acted  in  the  same  directions  as  wind  loads  during  the  second  test.  This  factor  resulted  in 
six  of  the  eight  wind  turbine  loads  being  lower  for  the  second  test  being  lower  that  the  first  sea  state  3 
test.  Comparison  of  the  EOG  results  with  those  of  operational  test  3  (see  Table  9-8)  shows  that  that 
the  EOG  event  is  not  a  limiting  case  for  the  floating  wind  turbine  system. 
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Table  9-10  Summary  of  maximum  loading  for  three  operational  tests.  Percentages  represent  the 
difference  between  fixed  base  wind  turbine  and  three  floating  wind  turbine  cases:  no  waves,  sea  state 
3,  and  sea  state  3  (shifted  by  180°).  Negative  delta  values  indicate  fixed  base  loading  was  higher. 


LOADING 

Fixed  Base 
Loading 
(kNm) 

No  Waves 

SS3  -1*‘ 

sss-a""  1 

Loading 

(kN-m) 

Delta 

(%) 

Loading 

(kNm) 

Delta  (%) 

Loading 

(kNm) 

Delta 

(%) 

Tower  Base  Bending 
(x-axis) 

1532 

1909 

24.6% 

2063 

34.7% 

1889 

23.3% 

Tower  Base  Bending 
(y-axis) 

27820 

31100 

11.8% 

31180 

12.1% 

25640 

-7.8% 

Blade-Root  Edgewise 
Bending 

668 

653 

-2.3% 

643 

-3.7% 

665 

-0.4% 

Blade-Root  Flapwise 
Bending 

2120 

2116 

-0.2% 

2112 

-0.4% 

2061 

-2.8% 

Low-Speed  Shaft  Tip 
Bending  (y-axis) 

239 

266 

11.3% 

260 

8.9% 

241 

0.8% 

Low-Speed  Shaft  Tip 
Bending  (z-axis) 

-30 

-75 

145.6% 

-68 

122.5% 

-46 

50.7% 

Low-Speed  Shaft 
Torque 

814 

827 

1.5% 

810 

-0.6% 

817 

0.3% 

Tower-Head  Torsion 

-104 

-140 

34.5% 

-132 

27.3% 

-109 

4.5% 

The  load  imparted  by  the  extreme  operating  gust  was  transferred  to  the  tethers.  The  wind  transients 
were  directed  along  the  x-axis  and  therefore  generated  surge  forces  and  pitching  moments.  Lines 
attached  to  the  upwind  spoke  (line  3)  had  the  highest  line  stresses  during  the  EOG  tests.  None  of  the 
line  stresses  exceeded  15  percent  of  the  breaking  stress  for  the  lines.  The  lines  all  remained  in  tension 
throughout  all  of  the  tests. 

9.6  Extreme  Wave  Loading 

The  final  simulations  dealt  with  how  the  system  responded  to  large  amplitude  waves  while  the  wind 
turbine  was  shutdown.  This  state  would  be  associated  with  extreme  wind  Conditions  which  prevent 
normal  wind  turbine  operations.  The  floating  wind  turbine  was  assumed  to  have  an  operational 
requirement  that  necessitated  securing  the  wind  turbine  once  the  sea  state  reached  a  certain  level.  It 
was  important  to  understand  how  the  system  responds  to  extreme  waves  when  in  this  condition.  The 
response  of  the  system  to  an  extreme  wave  event  was  evaluated  with  the  wind  turbine  secured  with 
the  hub  height  wind  speed  greater  than  the  cutout  speed  of  the  wind  turbine.  The  reference  wind 


speed  (at  10  meters  above  the  water)  which  corresponded  to  this  hub  wind  speed  was  20  m/s. 
According  to  Table  6.1  of  Newman,  the  sea  state  would  be  approximately  7  for  winds  of  20  m/s.  The 
Pierson-Moskowitz  sea  state  7  spectrum  is  chamcterized  by  a  significant  wave  height  Hj/s  =  8.9 
meteiB  and  a  avemge  period  To  =  10.7  seconds.  These  values  were  used  in  (9.4)  to  determine  the 
design  peak  y  =  10.68  meters.  This  value  was  used  to  generate  the  wave  group  for  the  two  test  cases 
(see  Figure  9.15). 


Figure  9.15  Wave  elevation  record  for  extreme  wave  tests. 


The  extreme  wave  is  the  largest  wave  expected  to  pass  the  floating  wind  turbine  during  a  three  hour 
period.  The  value  used  for  y  had  a  reliability  factor  of  a  =  0.01,  meaning  that  fliere  was  only  a  one 
percent  chance  of  the  actual  highest  wave  exceeding  the  value  of  y .  The  maximum  wave  amplitude 
during  the  test  was  10.68  meters.  When  die  trough  of  the  wave  passed  the  floater,  the  entire 
submerged  portion  of  the  tower  and  approximately  5  meters  of  the  spar  were  temporarily  above  the 
surface  of  the  water.  This  led  to  die  generation  of  large  nonlinear  loads  on  the  floater.  To  show  this 
behavior,  the  combined  horizontal  force  on  the  submerged  portion  of  the  tower  and  the  spar  was 
plotted  in  Figure  9.16.  The  horizontal  forces  in  the  positive  x-direction  were  larger  than  those  in  the 
negative  x-direction  due  to  the  non-linear  effecte  of  the  wave  passing  the  tethered  floater. 


Figure  9.16  Horizontal  wave  loading  of  submerged  tower  and  spar  during  extreme  wave  test  1. 

Incident  wave  angle  =  0°. 

The  effect  of  the  wave  loading  on  the  system  was  determined  for  two  wave  incident  angles:  0°  and 
30°.  The  other  loading  on  the  system  was  due  to  random  winds  which  had  an  average  hub  height 
velocity  of  25.3  m/s.  The  winds  were  directed  along  the  positive  x-axis.  The  wind  record  also 
reflected  the  fact  that  y  and  z  wind  velocity  components  were  present.  The  combination  of  wind  and 
wave  loading  on  the  system  produced  loads  on  the  various  wind  turbine  components.  A  baseline  test 
was  also  run  with  the  floater  fixed  in  position  for  comparison  purposes!  The  baseline  wind  turbine 
loads  were  then  compared  to  the  loads  from  the  two  floating  cases.  The  baseline  system  was  only 
acted  on  by  wind.  For  all  cases  the  rotor  was  fixed  in  position  with  the  blades  pitch  =  90°  to 
minimized  wind  loads.  The  maximum  values  for  the  three  tests  are  charted  in  Figure  9.17. 
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Figure  9.17  Maximum  wind  turbine  loads  during  extreme  wave  tests.  Wind  turbine  secured.  Blade 
pitch  =  90°.  Reference  speed  of  20  m/s.  Incident  wave  amplitude  =  10.68  m.  Incident  wave  period  = 
1 1 ,62  seconds.  1)  Fixed  base  (only  wind  loading),  2)  floating  (incident  wave  angle  =  0°),  and  3) 

floating  (incident  wave  angle  =  30°) 

The  extreme  wave  events  caused  translational  and  rotational  motion  of  the  system.  In  some  cases, 
these  motions  increased  the  loading  on  the  wind  turbine  components.  In  a  few  instances  the  ability  of 
the  floater  to  move  when  the  blade  elements  were  acted  on  by  the  wind  reduced  the  component 
loading  values.  Once  again  the  largest  differences  in  loading  were  between  the  fixed  and  floating 
support  tower  base  bending  moments.  The  tower  moment  about  the  y-axis  when  the  waves  were 
directed  along  the  x-axis  (0°)  was  over  nine  times  the  value  of  the  fixed  base  case  with  only  wind 
loading  on  tiie  feathered  rotor  blades. 
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The  wind  turbine  in  these  tests  is  in  a  different  configuration  during  the  extreme  wave  tests,  therefore 
the  wind  loading  oh  the  system  is  very  different  from  the  other  tests.  However,  by  comparing  the 
values  in  Table  9-1 1  with  those  of  Table  9-8  for  the  three  lines  per  spoke  design,  it  can  be  seen  that 
only  the  tower  base  bending  moments  for  the  extreme  wave  tests  are  greater  than  the  associated 
baseline  and  the  sea  state  6  operational  loading  tests.  Therefore,  with  exception  of  the  tower  base, 
exceeding  the  maximum  allowed  loading  on  wind  turbine  would  not  be  a  limiting  factor  for  the 
chosen  extreme  wave  event.  If  the  system  could  withstand  normal  operational,  it  could  also  handle 
the  design  peak  extreme  wave. 


Table  9-1 1  Extreme  wave  loading  results.  The  two  extreme  wave  cases  are  compared  with  the 
loading  on  the  fixed  base  wind  turbine  (baseline).  Reference  wind  speed  =  20  m/s.  Blade  pitch  = 
90°.  Maximum  wave  amplitude  =-10.68  m.  Wave  period  =  10.7  seconds. 


Component 

Baseline 

Loading 
(kN  m) 

Extreme  Wave  Tests 

Wave  Angle  =  0* 

Wave  Angle  =  30“ 

Loading 
(kN  m) 

Change 

{%) 

Loading 
(kN  m) 

Change 

(%) 

Tower  Base  Bending  (x-axis) 

19370 

18170 

-6.2% 

30700 

58.5% 

Tower  Base  Bending  (y-axis) 

-5247 

-53010 

910.3% 

-46530 

786.8% 

Blade-Root  Edgewise 
Bending 

333 

349 

4.8% 

382 

14.8% 

Blade-Root  Flapwise 

Bending 

1307 

1414 

8.2% 

1411 

8.0% 

Low-Speed  Shaft  Tip 

Bending  (y-axis) 

-402 

-425 

5.9% 

-419 

4.4% 

Low-Speed  Shaft  Tip 

Bending  (z-axis) 

361 

352 

-2.6% 

361 

-0.2% 

Low-Speed  Shaft  Torque 

1130 

1105 

-2.2% 

1070 

-5.3% 

Tower-Head  Torsion 

458 

454 

-0.9% 

620 

35.3% 

The  line  configuration  for  the  design  which  was  used  in  the  numerical  simulations  was  relatively  stiff. 
The  lines  and  spokes  resisted  the  pitching  moments  caused  by  the  wind  and  waves.  As  shown  in 
Figure  9.18,  even  when  the  system  was  rapidly  translating  the  mooring  system  maintained  the  vertical 
alignment  of  floater  and  tower. 


Figure  9.18  Surge  response  during  extreme  wave  event.  Left:  graph  of  surge  response.  Right  Still 

image  of  animation  of  simulation  at  time  =  279.7  seconds  (time  of  maximum  deflection). 

The  other  factor  which  was  evaluated  during  the  extreme  wave  cases  was  the  line  stress.  In  order  for 
the  system  design  to  be  acceptable  for  the  chosen  environment  the  lines  must  not  exceed  their 
breaking  stress.  A  safety  factor  would  also  be  incorporated  into  the  analyses.  None  of  the  line 
stresses  exceeded  15  percent  of  the  breaking  stress  (857.6  x  10*  Pa).  The  minimum  line  stress  during 
the  extreme  wave  simulations  was  found  to  be  76  x  10*  Pa.  If  any  of  the  line  stresses  had  gone  to 
zero  during  the  tests  the  pitch  angle  would  have  increased.  Opposite  lines  absorbed  moments  placed 
on  the  system.  Lines  1  and  3  (on  the  x-axis)  resist  pitch  motion  and  lines  2  and  4  (on  the  y-axis) 
resist  roll  motion.  As  shown  in  Figure  9.19,  during  the  first  test  the  lines  stresses  for  lines  1  and  3 
oscillated  with  line  1  stress  increasing  as  line  3  stress  decreased  and  vice  versa.  During  the  first  test 
the  waves  propagated  along  the  x-axis  and  therefore  created  no  roll  moments.  For  this  reason  lines  2 
and  4  stresses  only  had  to  absorb  roll  moments  produced  by  wind  effects  on  the  blades  and  inertial 
loads.  During  the  second  test  (incident  wave  angle  =  30°),  the  waves  produced  both  large  roll  and 
pitching  moments.  Lines  1  and  3  absorbed  the  pitch  moments  and  lines  2  and  four  absorbed  the  roll 
moments,  and  their  line  stresses  changed  accordingly. 


Figure  9.19  Line  stress  response  during  extreme  wave  events.  Wave  angles  =  0°  and  30°. 

During  the  extreme  wave  event  which  was  tested  all  of  the  wind  turbine  components  loading  values 
were  below  those  determined  during  the  normal  operations  test  cases.  The  line  stresses  caused  by  the 
extreme  wave  events  were  well  below  the  limitations  for  the  tethers.  Also,  the  lines  remained  in 
tension  throughout  which  prevented  excessive  rotational  motion  of  the  system.  Therefore,  it  was 
proven  that  the  tested  extreme  wave  event  would  not  be  a  limiting  case  for  this  configuration  of  the 
floating  wind  turbine  system. 


10  Conclusions  and  Future  Work 


The  thesis  focused  on  floating  offshore  wind  turbines.  The  first  goal  of  the  woric  was  to  develop  a 
method  to  analyze  a  floating  wind  turbine  system.  The  second  goal  was  to  determine  a  promising 
design  for  the  floating  support  which  could  allow  for  the  safe  Operation  of  a  wind  turbine  in  realistic 
wind  and  wave  environments.  Determining  the  dynamic  properties  of  the  system  was  the  third  goal 
of  this  research.  The  final  goal  was  evaluate  the  effects  of  normal  operations  and  extreme  wind  and 
wave  events  on  both  the  wind  turbine  and  tethers. 

In  order  to  support  a  wind  turbine  and  its  own  weight  the  floater  must  have  a  minimum  displacement. 
The  center  of  gravity  of  the  wind  turbine  is  well  above  the  free  surface  this  leads  to  a  problem  with 
the  center  of  gravity  being  above  the  center  of  buoyancy  for  floater.  This  arrangement  would  be 
unacceptable  and  could  be  overcome  by  either  increasing  the  size  and  weight  of  the  floater  to  bring 
die  system  center  of  gravity  down  below  the  center  of  buoyancy  or  using  tethere  to  maintain  the 
system  in  an  upright  position.  The  former  design  would  be  much  larger  than  the  later  design  and 
subsequently  more  expensive.  It  was  assumed  that  the  support  system  should  be  relatively  stiff  in  the 
pitch  and  roll  modes  of  motion  in  order  not  to  put  undo  stresses  on  the  wind  turbine  portion  of  the 
system.  Other  possible  designs  which  had  multiple  wind  turbines  supported  on  a  single  floater  were 
deemed  to  be  cost  prohibitive  and  would  be  subjected  to  very  large  wave  loads. 

The  design  chosen  for  analysis  in  this  work  was  a  tension  leg  spar  buoy,  which  consisted  of  a  vertical 
cylinder  most  of  the  displacement  volume  and  spokes  with  tethers  under  tension  attaching  the  floater 
to  the  seafloor.  This  design  required  smaller  displacement  volumes  and  hence  lower  material  costs 
for  the  floater  than  for  a  design  which  had  the  system  center  of  gravity  below  the  center  of  buoyancy. 
The  use  of  spokes  and  tethers  under  tension  creates  a  support  system  that  has  the  potential  to  be  very 
stiff  in  the  pitch  and  roll  modes  of  motion  depending  on  the  type  of  lines  used  for  the  tethere.  A 
design  parameter  analysis  of  the  system  was  conducted  for  the  system  and  is  included  in  Appendix  A. 

The  analysis  method  used  for  the  floating  wind  turbine  incorporated  models  of  a  wind  turbine  and 
floater.  These  two  systems  were  analyzed  in  a  fully  coupled  manner  in  the  time  domain.  The  wind 
turbine  was  modeled  using  a  blade  element  method.  Tbe  loads  on  the  individual  blade  elements  were 
calculated  using  the  AeroDyn  subroutine.  The  loads  on  the  floater  were  calculated  using  the 
JltrJoads  subroutine.  These  loadings  were  determined  for  the  system  model  using  the  ADAMS 
Solver  dynamics  program.  The  elastic  and  damping  properties  of  the  wind  turbine  and  tower  were 


incorporated  into  the  ADAMS  dynamic  model.  The  dynamic  model  of  the  system  also  included  the 
behavior  of  the  tethers. 

The  wave  loads  on  the  floater  were  calculated  using  a  combination  of  equations  of  G1  Taylor  and 
Morrison.  The  inertial  loads  on  the  floater  were  found  using  GI  Taylor.  Viscous  drag  calculations  on 
the  cylinders  and  spokes  on  the  system  were  determined  by  Morrison’s  equation.  These  calculations 
were  carried  out  in  the  time  domain  so  no  linear  approximations  of  drag  were  required.  The  operation 
of  the  floater  loading  subroutine  was  verified  for  small  amplitude  waves  with  linear  calculations  using 
a  floater  which  did  not  include  a  wind  turbine. 

Time  varying  wave  elevations  produce  nonlinear  wave  forces.  These  forces  result  in  a  mean 
displacement  of  the  system  in  the  direction  of  wave  propagation.  These  nonlinear  loads  are  modeled 
using  a  linear  approximation  method  described  in  Faltinsen,  1 990  [13]. 

Numerical  simulations  were  conducted  for  free  decay  tests  of  the  system  in  the  different  modes  of 
motion.  The  non-dimensional  linear  (wind  turbine)  and  quadratic  (viscous  drag)  damping  coefficients 
for  the  floater  and  wind  turbine  were  evaluated.  Effects  of  wind  speed  and  rotor  blade  pitch  angle 
were  investigated.  The  natural  periods  of  the  system  during  the  free  decay  tests  corresponded  well 
with  those  predicted  using  linear  theory. 

The  wind  turbine  damping  was  also  shown  to  be  linear.  The  wind  turbine  was  found  to  be  a  very 
large  source  of  damping  for  pitch  motion.  The  quadratic  damping  of  the  system  was  the  result  of 
viscous  drag.  Viscous  drag  dominated  the  initial  damping  of  the  system  for  most  of  the  free  decay 
simulations.  Blade  pitch  had  a  large  effect  on  the  system  damping.  This  fact  could  possibly  be 
incorporated  into  future  control  systems  to  change  damping  behavior  of  the  system. 

The  effects  of  normal  operating  conditions  on  various  components  of  the  wind  turbine  were 
evaluated.  Three  different  combinations  of  reference  wind  speed  and  sea  state  were  used.  The 
maximum  loading  and  fatigue  response  of  the  floating  wind  turbine  system  was  compared  to  that  of 
the  same  wind  turbine  with  a  fixed  base.  The  initial  floating  design  was  shown  to  have  increased 
maximum  loading  values  and  markedly  reduced  fatigue  life  when  compared  to  the  fixed  base  design. 
The  number  and  types  of  lines  used  to  tether  the  structure  were  changed  to  demonstrate  how  the 
initial  problems  could  be  overcome.  The  stiffness  of  the  system  in  pitch  and  roll  was  increased  by 
marking  the  tethers  stiffer.  The  altered  design  showed  marked  improvements  in  both  maximum 
loading  and  fatigue  response.  This  process  demonstrated  how  design  optimization  could  be  used  to 
improve  system  performance. 


The  final  tests  involved  determining  the  response  of  the  floating  wind  turbine  systems  to  extreme 
wind  and  wave  events.  The  extreme  operating  gust  case  described  in  lEC  61400-1  was  performed  for 
both  floating  and  fixed  base  wind  turbines.  The  maximum  loading  values  for  the  wind  turbine  were 
compared.  The  maximum  loading  was  determined  for  an  extreme  wave  event  with  the  wind  speed  at 
above  the  cutout  speed  for  the  wind  turbine.  With  one  exception,  the  values  for  maximum  load 
during  the  extreme  wind  and  wave  tests  were  found  to  be  less  than  those  determined  during  the 
normal  operations  cases  or  less  than  that  of  the  fixed  base  wind  turbine  results.  The  loading  of  the 
tower  base  was  identified  as  the  load  which  had  flie  largest  differences  between  the  fixed  base  and 
floating  support  loads. 

The  some  of  the  behaviors  of  the  tension  leg  spar  buoy  predicted  by  this  research  could  be  verified  by 
model  testing.  Unfortunately,  accurately  it  is  not  possible  to  build  a  scaled  model  of  an  operational 
wind  turbine  and  attach  it  to  a  scaled  floater.  Uie  behavior  of  the  wind  turbine  would  not  realistically 
mimic  the  response  of  a  large  scale  machine.  However,  some  useful  data  could  be  gathered  using 
scaled  models.  The  best  way  to  prove  the  viability  of  the  chosen  floating  support  would  be  to  build  a 
prototype  and  operate  it  in  an  offshore  environment.  Yet,  the  fully  coupled  dynamic  method 
developed  in  this  thesis  is  expected  to  be  accurate  since  the  independent  validation  of  the  wind 
turbine  modeling  methods  (AeroDyn,  FAST,  and  ADAMS)  as  well  as  the  fundamental  wave  loading 
algorithms  used  in  the  present  research  has  been  performed  extensively  in  earlier  studies. 

The  use  of  active  controls  is  important  for  large  wind  turbine  systems.  These  systems  are  designed  to 
improve  power  output  and  mitigate  adverse  loading  conditions  on  machine  components.  When 
mounted  on  a  floating  platform  the  control  systems  will  have  to  be  adapted  to  also  respond  to 
undesirable  loading  situations  of  the  tethers.  Further  research  is  required  to  better  understand  how 
theses  systems  might  best  be  utilized  for  a  floating  wind  turbine  system. 

In  conclusion,  the  fully  coupled  dynamic  analysis  method  proved  to  be  a  useful  tool  in  the  evaluation 
of  a  floating  offshore  wind  turbine  system.  The  tension  leg  spar  buoy  was  shown  to  be  a  promising 
design  for  supporting  wind  turbines  in  water  depths  that  prohibit  the  use  of  towers  fixed  to  the  sea 
floor.  In  the  future  the  design  should  be  further  analyzed  to  include  a  detailed  structural  model  for  the 
floater.  An  optimization  of  the  floater  configuration  would  be  required  in  onier  to  minimize  the  cost 
and  maximize  the  system  performance.  The  optimized  floater  design  should  then  be  evaluated  using 
the  techniques  developed  in  this  research.  Alternate  tether  /  mooring  line  configurations  should  also 
be  investigated. 
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Appendix  A:  Parametric  Design  Issues 


It  is  necessary  to  understand  how  changing  floater  design  parameters  affect  the  system  response,  in 
order  for  the  floater  design  evaluated  in  this  research  to  be  adapted  for  specific  wind  turbines  and 
environmental  conditions.  Figure  A.0.1  is  a  diagram  of  the  tension  leg  spar  buoy.  The  design 
parameters  that  are  used  to  describe  the  floater  are: 

•  spar  diameter  •  spoke  length 

•  spar  length  •  spoke  width 

•  submerged  tower  diameter  •  line  stiffness 


Figure  A.0.1  Diagram  of  Floating  Support 

The  line  stiffness  is  determined  by  the  type  of  line  used  (steel,  HPME,  polyester,  etc)  and  the 
diameter  of  the  line.  The  design  parameters  affect  the  restoring,  damping,  and  inertia  terms  in  the 
system  equations  of  motion.  Most  of  the  parameters  all  affect  the  wave  excitation  forces  acting  on 
the  submerged  portions  of  the  structure.  The  water  depth  has  an  effect  on  the  restoring  coefficients 
since  it  drives  the  length  of  the  tethers. 

It  is  necessary  to  understand  how  the  floater  design  can  be  modified  to  achieve  a  desired  change  in 
response.  The  wave  environment  in  which  the  system  must  operate  is  a  key  design  consideration  for 


the  system.  Changing  one  floater  design  parameter  can  affect  the  system  response  in  many  ways.  For 
example,  changing  the  floater  diameter  changes  the  excess  buoyancy  of  the  system  which  changes  the 
natural  periods  of  the  system  for  the  surge,  sway,  and  yaw  modes  of  motion.  This  change  has  small 
affects  on  the  roll  and  pitch  restoring  coefficients.  The  amount  of  wave  excitation  changes  when  the 
floater  diameter  is  changed. 

Natural  Periods 


Frequently^  designers  are  concerned  with  the  natural  periods  of  a  system.  For  waterborne  shiictures  it 
is  normal  to  design  the  system  such  that  its  natural  periods  lie  either  about  or  below  the  range  where 
wave  action  might  excite  resonance.  This  section  examines  the  affects  the  design  parameter  of  the 
floater  have  on  the  systems  natuml  periods.  The  general  equations  used  to  determine  the  natural 
periods  are: 


^natural 

I  C 

Jm+a 


where  C  is  the  restoring  coefficient.  Mis  the  inertia  term,  and  A  is  the  added  mass.  Due  to  symmetry 
of  the  floater  the  natural  periods  for  surge  and  sway  are  identical  and  roll  and  pitch  are  identical.  For 
the  following  analysis  the  design  parameters  are  varied  about  a  reference  design.  This  reference 
design  is  the  same  as  that  used  in  the  extreme  event  cases  discussed  in  the  previous  chapter.  Table 
A-1  and  Table  A-2  list  the  properties  of  the  reference  system. 


Table  A-1  Dimensions  and  properties  of  floater  used  in  operational  tests. 


Spar  Diameter 


Spar  Length 


Tower  Diameter 


Tower  Dra 


Spoke  Length 


Spoke  Width 
Buoyant  Spokes 


Line  Diameter 


12  m 

Number  of  Lines 

8m 

Line  Minimum  Break  Load 

5  m 

Maximum  Line  Extension 

Line  Length  (un-stretched) 
Water  Depth 
Stracture  Thickness 
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104mm 


3 / spoke 


743  tonnes 


180.76  m 


200  m 
2.53  cm 


Table  A-2  Calculated  attributes  of  the  floating  wind  turbine  used  in  the  operational  tests 


Mass  of  Floater 

210435.4  kg 

Center  of  Gravity  (system)  CG 

16.98  m 

Buoyant  Mass 

1228424.6  kg 

Center  of  Buoyancy  CB 

-10.89  m 

Total  Mass  of  System 

411494.5  kg 

CG-CB 

27.87  m 

Reserve  Buoyancy 

198.5  % 

Surge  Natural  Period 

39.0186  sec 

Initial  Line  Stress 

235.4  MPa 

Heave  Natural  Period 

1.6591  sec 

Breaking  Line  Stress 

857.7  MPa 

Pitch  /  Roll  Natural  Period 

3.3561  sec 

Center  of  Gravity  (floater) 

-13.82  m 

Yaw  Natural  Period 

10.51  sec 

The  important  parameters  used  in  determining  the  natural  periods  of  systems  are  the  restoring 
coefficients,  the  inertia  terms  of  the  structure,  and  the  added  mass.  The  following  equations  are 
relevant  to  the  discussions  of  natural  periods: 
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Spar  Diameter 


The  first  design  parameter  which  is  evaluated  is  spar  diameter.  Changing  the  spar  diameter  changes 
the  cross-sectional  area  of  the  cylinder,  which  increases  the  added  mass  in  every  mode  except  yaw. 
The  restoring  coefficients  are  also  affected  since  the  buoyancy  of  the  system  changes.  The  spar 
diameter  is  varied  from  12  meters  to  8  meters  and  the  natural  periods  determined.  The  spoke  length  is 
changed  to  keep  the  lines  at  the  same  location  relative  to  the  center  of  the  structure. 

Figure  A.0.2  shows  how  increasing  the  diameter  causes  the  surge,  sway,  and  yaw  natural  periods  to 
decrease.  This  decrease  is  due  to  the  fact  that  the  as  the  diameter  increases  the  amount  of  reserve 
buoyancy  increases  (even  though  the  weight  went  up  also).  The  restoring  coefficients  for  all  three  of 
these  modes  of  motion  are  dependant  on  the  amount  of  reserve  buoyancy.  As  the  restoring 
coefficients  increase  the  natural  periods  decrease. 

The  pitch,  roll,  and  heave  natural  periods  increase  when  the  spar  diameter  increases.  This  is  due  to 
the  fact  that  the  heave  added  mass  associated  with  these  modes  of  motion  also  increases  with  the  spar 
diameter.  And  since  the  restoring  coefficients  are  not  greatly  affected  by  the  change  in  diameter 
(heave  restoring  increases  slightly  with  increasing  diameter),  the  natural  periods  increase  with  the 
spar  diameter. 


Figure  A.0.2  Effect  of  spar  diameter  on  natural  periods. 


Spar  Height 

Spar  has  a  similar  effect  on  the  systems  natural  periods  as  spar  diameter  on  all  of  the  modes  of  motion 
except  for  heave  (see  Figure  A.0.3).  Changing  the  spar  length  changes  the  buoyancy  and  mass  of  the 
system.  The  buoyancy  changes  up  linearly  unlike  spar  diameter.  The  added  mass  terms  for  surge  and 
sway  are  affected  while  yaw  added  mass  is  unaffected.  The  net  result  of  changing  the  spar  length  are 
natural  periods  of  surge,  sway,  and  yaw  which  decreases  less  than  linearly  as  length  increases. 


Figure  A.0.3  Effects  of  spar  height  on  natural  periods. 


Heave  is  relatively  unaffected  by  changing  the  length  of  the  cylinder  in  water  depths  of  200  meters. 
The  mass  of  the  system  increases  with  spar  height,  therefore  the  inertia  term  gets  bigger.  Small 
changes  in  line  length  affect  the  spring  constant  of  the  lines  and  hence  the  heave  restoring  coefficient. 
However,  the  increase  in  the  inertia  term  dominates  and  heave  natural  period  goes  up  with  increasing 
spar  length.  The  roll  and  pitch  natural  periods  are  more  affected  by  changes  in  length  since  the  added 
mass  terms  go  up  proportional  to  the  square  of  the  length.  The  moments  of  inertia  also  change.  The 
combined  effect  of  increased  inertia  and  added  mass  in  roll  and  pitch  cause  the  natural  periods  to 
increase  for  longer  spar  lengths. 


Spoke  Length 

The  use  of  spokes  for  the  tension  leg  spar  buoy  increases  the  stiffness  of  the  system  in  all  three  of  the 
rotational  modes  of  motion.  As  was  shown  in  the  previous  chapter,  the  increased  stiffness  reduced 
maximum  and  fatigue  loading  on  the  floating  wind  turbine.  Unfortunately,  as  the  spoke  length 
increases  challenges  arise  in  mounting  very  long  spokes  on  a  mlatively  small  spar.  Design 
optimization  should  be  used  to  balance  the  additional  cost  of  longer  spokes  to  increasing  the  stiffiiess 
of  the  lines.  This  decision  will  be  based  on  the  water  depth  and  the  size  of  the  machine  being 
supported  by  die  floater. 

Since  the  spokes  are  modeled  as  closed  box  structures,  the  inertia  terms  for  every  mode  of  motion  are 
affected  by  changing  the  spoke  length.  The  added  mass  terms  for  all  of  the  modes  increase  with 
spoke  length.  These  increases  in  the  inertia  and  added  mass  terms  cause  all  of  die  modes  except  pitch 
and  roll  to  increase  with  spoke  length.  The  yaw  restoring  coefficient  increases  with  spoke  length,  but 
not  as  much  as  the  added  mass  and  inertia  terms  increase.  A  different  model  of  spoke  which  had  less 
impact  on  the  inertia  and  added  mass  terms  would  result  in  the  yaw  natural  period  decreasing  with 
increasing  spoke  length. 

As  shown  in  Figure  A.0.4,  the  pitch  and  roll  natural  periods  decrease  less  than  linearly  with  the 
increases  in  spoke  lengdi.  The  spoke  length  greatly  affects  die  pitch  restoring  coefficient,  which  is 
die  dominant  term  for  these  modes  natural  period. 


Figure  A.0.4  Effect  of  spoke  length  on  natural  periods. 


Spoke  Width 


Spoke  width  affects  all  of  the  modes  in  a  similar  manner  (see  Figure  A.0.5).  An  increase  in  spoke 
width  increases  the  mass  and  moments  of  inertia  of  the  system.  Larger  spoke  cross-sections  have 
more  added  mass.  The  restoring  coefficients  in  surge,  sway,  and  yaw  go  down  as  the  reserve 
buoyancy  increases.  Therefore,  these  effects  combine  to  increase  the  natural  period  with  increases  in 
spoke  width. 
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Figure  A.0.5  Effect  of  spoke  width  on  natural  periods. 

Number  of  Lines 

Changing  the  number  of  lines  parameter  is  equivalent  to  increasing  the  spring  constant  acting  at  the 
end  of  each  spoke.  Changing  the  line  diameter  or  material  type  would  have  similar  effects.  The 
surge,  sway,  and  yaw  restoring  coefficients  are  constant  when  lines  are  changed  since  they  depend  on 
the  amount  of  reserved  buoyancy  (generates  line  tension)  and  the  length  of  the  lines.  The  number  of 
lines  has  a  negligible  effect  on  the  inertia  properties  of  the  system.  Surge,  sway,  and  yaw  natural 
periods  are  unaffected  by  the  changes  in  the  lines. 

Pitch  and  roll  natural  periods  are  greatly  affected  by  the  number  of  lines.  For  most  floater  designs  the 
line  tension  and  relaxation  components  of  the  pitch  and  roll  restoring  coefficients  determine  a 
majority  of  the  stiffiiess  (see  equation  (4.6)).  Changing  the  number  of  lines  affects  the  restoring 
coefficients.  Therefore,  pitch  and  roll  natural  frequencies  change.  Figure  A.0.6  is  a  graph  showing 
how  increasing  the  number  of  lines  reduces  the  roll  and  pitch  natural  periods.  Similarly,  the  heave 
restoring  coefficient  changes  when  the  number  of  lines  is  changed,  this  causes  the  natural  period  to 
change. 


Appendix  B:  Input  file  for  floater  verification  test 


4.0 

4.0 

0.0 

10.0 

0.0 

5.0 

1.0 

0.0253 

0.050 

1 

200.0 

0.045 

200.0 

0.0 

0.0 

1 

.true. 

1.0 

0.75 


!  Floater  Diameter  ALL  UNITS  SI 
!  Tower  Diameter 
!  Depth  of  tower 
!  Length  of  cylinder  (Height) 

!  Height  of  spokes  on  cylinder  (0  =  bottom) 

!  Length  of  spokes 

!  Characteristic  width  of  spokes  (square  cross-section) 
!  Thickness  of  spar 
!  Line  Diameter  (m) 

!  Number  of  lines 
!  Breaking  Load  in  tonnes  743 
!  Maximum  line  Extension 
!  Water  depth 
!  Wave  Direction  (degrees) 

!  Wave  Direction  Range  (+/-  this  value) 

!  Number  of  Bands  (integer  value) 

!  Buoyant  spokes 
!  Wave  Amplitude  (meter) 

!  Wave  Frequency  (rad/sec) 


========Line  Length  Calculations  ~  ~  — - ======= 

229.40  190.00  1.0785  Initial  Values  of  KofLine  kN/m,  Line  Length,  and  L 
230.70  188.92  1.0725  Final  Values  of  KofLine  kN/m,  Line  Length,  and  Lin 

=========Calculated  Floater  Properties — ^ - — 

32362.0  Spar  Mass 

24893.8  Mass  Cylinder 
2489.4  Mass  Spar  Plate 

0.0  Submerged  Tower  Mass 

16640.3  Mass  of  Spokes 

49002.3  Total  Floater  Mass 
201059.2  Mass  of  Wind  Turbine 

49002.3  System  Mass 

149887.8  Buoyant  Mass 

100885.5  Reserve  Buoyancy 
205.88  Percent  Reserve  Buoyancy 
0.0079  Total  Area  of  Lines 
230.7  Spring  Constant 

126.01  Intial  Line  Stress  (MPa) 

1 .07  Intial  Line  Stretch 

========  System  Properties  ===^^============= 

0.0000  CG  of  Wind  Turbine 

-6.6979  CG  of  Floater 

-6.6979  CG  of  System 

-1.0114  BG  of  System 

141537.5  AllandA22 
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46183.0  A33 

14281945.8  A44andA55 
546258.4  A66 

1140080.9  Ixx 

1140080.9  lyy 
451478.0  Izz 

5208.9  CllandC22 

1049673.8  C33 

27364352.8  C44andC55 

36462.2  C66 

- -  - - ==  Natural  Frequencies  of  Floater  System 

0. 1653  Surge  and  Sway  Natural  Frequencies  (rad/sec) 

3.3208  Heave  Natural  Frequency(rad/sec) 

1.3321  Roll  Natural  Frequency(rad/sec) 

1.3321  Pitch  Natural  Frequency  (rad/sec) 

0.1912  Yaw  Natural  Frequency  (rad/sec) 

-  -  — -  ^=  Natural  Periods  of  Floater  System  === 

38.0015  Surge  and  Sway  Natural  Period  (sec) 

1.8921  Heave  Natural  Period  (sec) 

4.7169  Roll  Natural  Period  (sec) 

4.7 1 69  Pitch  Natural  Period  (sec) 

32.8675  Yaw  Natural  Period  (sec) 


Appendix  C:  1.5  MW  Baseline  FAST  Input  File 


- FAST  INPUT  FILE - — _ 

FAST  certification  Test  #13:  WindPACT  1 .5  MW  Baseline  with  many  DOFs  with  VS  and  VP  and  FF  turbulence. 

Model  properties  from  'TnputDatal  .5A08V07adm.xls"  (from  C.  Hansen)  with  bugs  removed.  Compatible  with  FAST  v4  3 

- SIMULATION  CONTROL - 

False  Echo  -  Echo  input  data  to  "echo. out”  (switch) 

2  ADAMSPrep  -  ADAMS  preprocessor  mode  {1:  Run  FAST,  2:  use  FAST  as  a  preprocessor  to  create  an  ADAMS  model,  3:  do 
both}  (switch) 

2  AnalMode  -  Analysis  mode  { 1 :  Run  a  time-marching  simulation,  2:  create  a  periodic  linearized  model}  (switch) 

3  NumBl  -  Number  of  blades  (-) 

1000.0  TMax  -  Total  run  time  (s) 

0.05  DT  -  Integration  time  step  (s) 

- TURBINE  CONTROL - 

2  PCMode  -  Pitch  control  mode  (0:  none,  1:  power  control,  2:  speed  control}  (switch) 

5.0  TPCOn  -  Time  to  enable  active  pitch  control  (s) 

1  VSContrl  -  Variable-speed  control  {0:  none,  1:  simple  VS,  2:  user-defined  VS}  (switch) 

1 800.0  RatGenSp  -  Rated  generator  speed  for  simple  variable-speed  generator  control  (HSS  side)  (rpm)  [used  only  when  VSContrl=  1  ] 

0.002456  Reg2TCon  -  Torque  constant  for  simple  variable-speed  generator  control  in  Region  2  (HSS  side)  (N-m/rpm"^2)  fused  onlv 
whenVSContrl^l] 

1  GenModel  -  Generator  model  { 1 :  Simple,  2:  Thevenin,  3:  User  Defined)  [used  only  when  VSContrl=0]  (-) 

True  GenTiStr  -  Method  to  start  the  generator  {T:  timed  using  TimGenOn,  F:  generator  speed  using  SpdGenOn}  (switch) 

True  GenTiStp  -  Method  to  stop  the  generator  {T:  timed  using  TimGenOf,  F:  when  generator  power  =  0}  (switch) 

9999.9  SpdGenOn  -  Generator  speed  to  turn  on  the  generator  for  a  startup  (HSS  speed)  (rpm) 

0.0  TimGenOn  -  Time  to  turn  on  the  generator  for  a  startup  (s) 

9999.9  TimGenOf  -  Time  to  turn  off  the  generator  (s) 

9999.9  THSSBrDp  -  Time  to  initiate  deployment  of  the  HSS  brake  (s) 

9999.9  TiDynBrk  -  Time  to  initiate  deployment  of  the  dynamic  generator  brake  [CURRENTLY  IGNORED]  (s) 

9999.9  TTpBrDp(l)  -  Time  to  initiate  deployment  of  tip  brake  1  (s) 

9999.9  TTpBrDp(2)  -  Time  to  initiate  deployment  of  tip  brake  2  (s) 

9999.9  TTpBrDp(3)  -  Time  to  initiate  deployment  of  tip  brake  3  (s)  [unused  for  2  blades] 

9999.9  TBDepISp(l)  -  Deployment-initiation  speed  for  the  tip  brake  on  blade  1  (rpm) 

9999.9  TBDepISp(2)  -  Deployment-initiation  speed  for  the  tip  brake  on  blade  2  (rpm) 

9999.9  TBDepISp(3)  -  Deployment-initiation  speed  for  the  tip  brake  on  blade  3  (rpm)  [unused  for  2  blades] 

9999.9  TPitManS(l)  -  Time  to  start  override  pitch  maneuver  for  blade  1  and  end  standard  pitch  control  (s) 

9999.9  TPitManS(2)  -  Time  to  start  override  pitch  maneuver  for  blade  2  and  end  standard  pitch  control  (s) 

9999.9  TPitManS(3)  -  Time  to  start  override  pitch  maneuver  for  blade  3  and  end  standard  pitch  control  (s)  [unused  for  2  blades] 

9999.9  TPitManE(l )  -  Time  at  which  override  pitch  maneuver  for  blade  1  reaches  final  pitch  (s) 

9999.9  TPitManE(2)  -  Time  at  which  override  pitch  maneuver  for  blade  2  reaches  final  pitch  (s) 

9999.9  TPitManE(3)  -  Time  at  which  override  pitch  maneuver  for  blade  3  reaches  final  pitch  (s)  [unused  for  2  blades] 

23.0  BlPitch(i)  -  Blade  1  initial  pitch  (degrees) 

23.0  BIPitch(2)  -  Blade  2  initial  pitch  (degrees) 

23,0  BlPitch(3)  -  Blade  3  initial  pitch  (degrees)  [unused  for  2  blades] 

23.0  BlPitchF(l)  -  Blade  1  final  pitch  for  pitch  maneuvers  (degrees) 

23.0  BlPitchF(2)  -  Blade  2  final  pitch  for  pitch  maneuvers  (degrees) 

23.0  BlPitchF(3)  -  Blade  3  final  pitch  for  pitch  maneuvers  (degrees)  [unused  for  2  blades] 

- ENVIRONMENTAL  CONDITIONS - 

9.80665  Gravity  -  Gravitational  acceleration  (m/s^2) 

- FEATURE  SWITCHES - 

True  FlapDOFl  -  First  flapwise  blade  mode  DOF  (switch) 

True  FlapDOF2  -  Second  flapwise  blade  mode  DOF  (switch) 

True  EdgeDOF  -  First  edgewise  blade  mode  DOF  (switch) 

False  TeetDOF  -  Rotor-teeter  DOF  (switch)  [unused  for  3  blades] 

True  DrTrDOF  -  Drivetrain  rotational-flexibility  DOF  (switch) 

True  GenDOF  -  Generator  DOF  (switch) 

False  TiltDOF  -  Nacelle-tilt  DOF  (switch) 

False  YawDOF  -  Yaw  DOF  (switch) 

True  TwFADOiFl  -  First  fore-aft  tower  bending-mode  DOF  (switch) 

True  TwFADOF2  -  Second  fore-aft  tower  bending-mode  DOF  (switch) 

True  TwSSDOFl  -  First  side-to-side  tower  bending-mode  DOF  (switch) 

True  TwSSDOF2  -  Second  side-to-side  tower  bending-mode  DOF  (switch) 

True  CompAero  -  Compute  aerodynamic  forces  (switch) 

False  CompNoise  -  Compute  aerodynamic  noise  (switch) 
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- - -—INITIAL  CONDITIONS----^ - — - — — - 

0.0  OoPDefl  -  Initial  out-of-p!ane  blade-tip  displacement,  (meters) 

0.0  IPDefl  -  Initial  in-plane  blade-tip  deflection,  (meters) 

0,0  TeetDefl  -  Initial  or  fixed  teeter  angle  (degrees)  [unused  for  3  blades] 

0.0  Azimuth  -  Initial  azimuth  angle  for  blade  1  (degrees) 

20.0  RotSpeed  -  Initial  or  fixed  rotor  speed  (rpm) 

-5,0  NacTilt  -  Initial  or  fixed  nacelle-tilt  angle  (degrees) 

0.0  NacYaw  -  Initial  or  fixed  nacelle-yaw  angle  (degrees)  **********************^ 

0.0  TTDspFA  -  Initial  fore-aft  tower-top  displacement  (meters) 

0.0  TTDspSS  -  Initial  side-to-side  tower-top  displacement  (meters) 

- - —TURBINE  CONFIGURATION—. - — — — ~ 

35.0  TipRad  -  The  distance  from  the  rotor  apex  to  the  blade  tip  (meters)  . 

1.75  HubRad  -  The  distance  from  the  rotor  apex  to  the  blade  root  (meters) 

1  PSpnEIN  -  Number  of  the  innermost  blade  element  which  is  still  part  of  the  pitchable  portion  of  the  blade  for  partial-span  pitch 

control  [I  to  BldNodes]  [CURRENTLY  IGNORED]  (-) 

0.0  UndSling  -  Undersling  length  [distance  from  teeter  pin  to  the  rotor  apex]  (meters)  [unused  for  3  blades] 

0.0  HubCM  -  Distance  from  rotor  apex  to  hub  mass  [positive  downwind]  (meters) 

-3.3  OverHang  -  Distance  from  yaw  axis  to  rotor  apex  [3  blades]  or  teeter  pin  [2  blades]  (metei^) 

^.1251  ParaDNM  -  Distance  parallel  to  shaft  from  yaw  axis  to  nacelle  CM  (meters) 

-0.2328  PeipDNM  -  Peipendicular  distance  from  shaft  to  nacelle  CM  (metem) 

S2.39  TowerHt  -  Height  of  tower  above  ground  level  (meters) 

1 .61  Twi2Shft  -  Vertical  distance  from  the  tower  top  to  the  yaw/shaft  intersection  (metera) 

0.0  TwrRBHt  -  Tower  rigid  base  height  (meters) 

0.0  Delta3  -  Delta-3  angle  for  teetering  rotors  (degrees)  [unused  for  3  blades] 

0.0  PreCon€(l) -Blade  1  cone  angle  (degrees) 

0.0  PreCone(2)  -  Blade  2  cone  angle  (degrees) 

0.0  PreCone(3)  -  Blade  3  cone  angle  (degrees)  [unused  for  2  blades] 

0.0  AzimBlUp  -  Asmutfi  value  to  use  for  I/O  when  blade  1  points  up  (degrees) 

- mass  and  INERTIA  - - — - - 

51.170E3  NacMass  -Nacelle  mass  (kg) 

15.148E3  HubMass  -Hub mass  (kg) 

0.0  .TipMass(l) -Tip-brake mass, blade  I  (kg) 

0.0  TjpMass(2)  -  Tip-brake  mass,  blade  2  (kg) 

0.0  TipMass(3)  -  Tip-brake  mass,  blade  3  (kg)  [unused  for  2  blades] 

49.130E3  NacYIner  -  Nacelle  inertia  about  yaw  axis  (kg  m^2)  ' 

58.720E3  NacTIner  -  Nacelle  inertia  about  tilt  axis  ^cg  m^2) 

53.036  Geniner  -  Generator  inertia  about  HSS  (kg  m^2) 

34.6W)E3  Hublner  -  Hub  inertia  about  teeter  axis  (kg  m^2)  [unused  for  3  blades] 

- - - —  DRIVETRAIN — — - - - - 

100.0  GBoxEff  -  Gearbox  efficiency  (%) 

95.0  GenEff  -  Generator  efficiency  [ignored  by  the  Thevenin  and  user-defln^  generator  models]  (%) 

87.965  GBRatio  -  Gearbox  ratio  (-) 

False  GBRevera  -  Gearbox  reversal  (T:  if  rotor  and  generator  rotate  in  opposite  directions}  (switch) 

9999.9  HSSBfTqF  -  Fully  deployed  HSS-brake  torque  (N-m) 

9999.9  HSSBrDT  -  Time  for  HSS-brake  to  reach  full  d^loyment  once  initiated  (sec) 

D>mBrkFi  -  File  confining  a  mech-gen-torque  vs  HSS-speed  curve  for  a  dynamic  brake  [CURRENTLY  IGNORED]  (quot^ 

string) 

5.6E9  DTTorSpr  -  Drivetrain  torsional  spring  (N-m/rad) 

1.0E7  DTTorDmp  -  Drivetrain  torsional  damper  (N-m/s) 

- - - — - - -  SIMPLE  INDUCTION  GENERATOR 

9999.9  SIG_SlPc  -  Rated  generator  slip  percentage  [>0]  (%)  Now  HSS  side! 

9999.9  SIG_SySp  -  Synchronous  (zero-torque)  generator  speed  [>0J  (rpm)  Now  HSS  side! 

9999.9  SIG^RtTq  -  Rated  torque  [>0]  (N-m)  Now  HSS  side! 

9999.9  SIG_PORt  -Pull-outratio(Tpullout/Trated)  [>l]  (-) 

— - -  THEVENIN-EQUIVALENT  INDUCTION  GENERATOR  — — — - 

9999.9  Xrc^Freq  -  Line  frequency  [50  or  60]  (Hz) 

9998  TEC_NPol  -  Number  of  poles  [even  integer  >  0]  (-) 

9999.9  TEC_SRes  -  Stator  resistance  [>0]  (ohms) 

9999.9  TEC_RRes  -  Rotor  resistance  [>0]  (ohms) 

9999.9  TEC_VLL  -  Line-to-!ine  RMS  voltage  (vote) 

9999.9  TEC_SLR  -  Stator  leakage  reactance  (ohms) 

9999.9  TEC_RLR  -  Rotor  leakage  reactance  (ohms) 

9999.9  TEC_MR  -  Magnetizing  reactance  (ohms) 

* - - — - TOWER  — — - - - 

10  TwrNodes  -  Number  of  tower  nodes  i^ed  for  analysis  (-) 

**BaseIinejrower.dat"  TwrFile  -  Name  of  file  containing  tower  properties  (quoted  string) 

- — NACELLE- YAW - - — - — - — — - — — 

0.0  YawSpr  -  Nacelle-yaw  spring  constant  (N-m/rad) 

0.0  YawDamp  -  Nacelle-yaw  constant  (N-m/rad/s) 

0.0  YawNeut  -  Neutral  yaw  position— yaw  spring  force  is  zero  at  this  yaw  (degrees) 
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- - NACELLE-TILT - 

0.0  TiltSpr  -  Nacelle-tilt  linear-spring  constant  (N-m/rad) 

0.0  TiltDamp  -  Nacelle-tilt  damping  constant  (N-m/rad/s) 

0.0  Tilts StP  “  Nacelle-tilt  soft-stop  position  (degrees) 

0.0  TiltHStP  -  Nacelle-tilt  hard-stop  position  (degrees) 

0.0  TiltSSSp  -  Nacelle-tilt  soft-stop  linear-spring  constant  (N-m/rad) 

0.0  TiltHSSp  -  Nacelle-tilt  hard-stop  linear-spring  constant  (N-m/rad) 

- ROTOR-TEETER - 

0  TeetMod  -  Rotor-teeter  spring/damper  model  (0:  none,  1 :  standard,  2:  user-defined)  (switch)  [unused  for  3  blades] 

0.0  TeetDmpP  -  Rotor-teeter  damper  position  (degrees)  [unused  for  3  blades] 

0.0  TeetDmp  -  Rotor-teeter  damping  constant  (N-m/rad/s)  [unused  for  3  blades] 

0.0  TeetCDmp  -  Rotor-teeter  rate-independent  Coulomb-damping  moment  (N-m)  [unused  for  3  blades] 

0.0  TeetSStP  -  Rotor-teeter  soft-stop  position  (degrees)  [unused  for  3  blades] 

0.0  TeetHStP  -  Rotor-teeter  hard-stop  position  (degrees)  [unused  for  3  blades] 

0.0  TeetSSSp  -  Rotor-teeter  soft-stop  linear-spring  constant  (N-m/rad)  [unused  for  3  blades] 

0.0  TeetHSSp  -  Rotor-teeter  hard-stop  linear-spring  constant  (N-m/rad)  [unused  for  3  blades] 

- TIP-BRAKE  — - - - - 

0.0  TBDrConN  -  Tip-brake  drag  constant  during  normal  operation,  Cd*Area  (m'^2) 

0.0  TBDrConD  -  Tip-brake  drag  constant  during  fully-deployed  operation,  Cd*Area  (m*2) 

0.0  TpBrDT  -  Time  for  tip-brake  to  reach  full  deployment  once  released  (sec) 

- : - BLADE - - - 

'’Baseline_Blade.dat'’  BldFile(l)  -  Name  of  file  containing  properties  for  blade  1  (quoted  string) 

"Baseline_Blade.dat"  BldFile(2)  -  Name  of  file  containing  properties  for  blade  2  (quoted  string) 

’’Baseline_Blade.daf'  BldFile(3)  -  Name  of  file  containing  properties  for  blade  3  (quoted  string)  [unused  for  2  blades] 

- AERODYN - 

"Testl  3_AD.ipt"  ADFile  -  Name  of  file  containing  AeroDyn  input  parameters  (quoted  string) 

- NOISE - 

"Noise.dat”  NoiseFile  -  Name  of  file  containing  aerodynamic  noise  input  parameters  (quoted  string) 

- ADAMS - 

"BaseHne_ADAMS.dat”  ADAMSFile  -  Name  of  file  containing  ADAMS-specific  input  parameters  (quoted  string) 

- LINEARIZATION  CONTROL - 

''Baseline_Linear.dat"  LinFile  -  Name  of  file  containing  FAST  linearazation  parameters  (quoted  string) 

- OUTPUT - 

True  SumPrint  -  Print  summary  data  to  ''<RootName>.fsm''  (switch) 

True  TabDelim  -  Generate  a  tab-delimited  tabular  output  file,  (switch) 

"ES10.3E2"  OutFmt  -  Format  used  for  tabular  output  except  time.  Resulting  field  should  be  10  characters,  (quoted  string)  [not  checked 
for  validity!] 

10.0  TStart  -  Time  to  begin  tabular  output  (s) 

1 0  DecFact  -  Decimation  factor  for  tabular  output  [  1 :  output  every  time  step]  (-) 

1 .0  Sttstime  -  Amount  of  time  between  screen  status  messages  (sec) 

0.99  ShftGagL  -  Distance  from  rotor  apex  [3  blades]  or  teeter  pin  [2  blades]  to  shaft  strain  gages  [positive  for  upwind  rotors]  (meters) 
0  NBlGages  -  Number  of  blade  nodes  that  have  strain  gages  for  output  [0  to  5]  (-) 

0  BldGagNd  -  List  of  blade  nodes  that  have  strain  gages  [  1  to  BldNodes]  (-) 

OutList  -  The  next  line(s)  contains  a  list  of  output  parameters.  See  OutList.txt  for  a  listing  ofavailable  output  channels  {-) 
”Windyxt,WindVyt,WindVzt”  -  Wind-speed  components 
”HorWndDir,VerWndDir''  -  Wind  directions 

'’BldPitch2”  -  Blade  2  pitch  angle 

"GenPwr,LSShftFxa"  -  Generator  Power,  Rotor  Thrust 

"IPDefll,  IPDefl2"  -  IP  blade  1,2  tip  deflections 

'TipAycl,TipAyc2”  -  Blade  IP  tip  acceleration 

”RootMxb2,  RootMyb2,  RootMzb2”  -  Blade  2  root  moments 
”LSShftFys,  LSShftFzs”  -  Non-rotating  LSS  shear  forces 
"LSSTipMys,  LSSTipMzs”  -  Non-rotating  LSS  bending  moments  at  the  shaft  tip 
"YawBrTDxp,  YawBrTDyp"  -  FA  and  SS  tower-top  deflections  wrt  tilted  tower 

"YawBrMxn,  YawBrMyn,  YawBrMzn"  -  Tower-top  /  yaw  bearing  roll,  pitch,  and  yaw  moments 
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